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STATISTICAL ASPECTS OF COLLISION WARNING SYSTEM DESIGN* t 
UNDER THE ASSUMPTION OF CONSTANT VELOCITY COURSES 


Judah Rosenblatt 


University of New Mexico 


INTRODUCTION 

The increase in aircraft speed, due to the advent of high-speed jet airplanes, has made 
it necessary to devise superior ways of avoiding mid-air collisions. Because so much distance 
is covered in the time it takes an aircraft to react to its controls, any maneuvers taken to avoid 
a collision must be executed well in advance of such a possibility. For this reason, a visual 
determination of the likelihood of collision is no longer practical. In a previous article, Pina 
and Landau |1], have proposed using radar equipment, and devised a scheme for taking and 
processing the observations to achieve a useful collision warning system. 

The purpose of this article is to find the best procedure for taking and processing these 
observations in a manner consistent with certain restraints, namely a limit on the number of 
available observations, which results from the limited detection range, and the existence of a 
required warning time. 

This article attempts to refine both the method of taking observations and the method of 
processing observations which have previously been proposed. 

Under the idealized assumptions made, these results are an extension and refinement of 
the work of Pina and Landau |1]. Even if the system proposed here is not used, possibly due 
to obvious weight restrictions, it will serve as a yardstick by which the performance of any 


proposed system can be measured. 


ASSUMPTIONS 
1. The two airplanes of interest each fly constant velocity (straight line, constant 
speed). (Thus, the important case of flight in the vicinity of aircraft terminals is 
not handled in this article.) 
The first airplane, P,, is equipped with the complete collision warning system, and 
the second airplane, Po, has at least a responder. 
P, can take measurements on range. Measurements on range are considered to have 


maximum error of (9g feet. (In the Boeing system o9= 50, for reasons to be given.) 


* 
Manuscript received May 8, 1961. 


t This study was initiated during the summer of 1957 while the author was a summer faculty 
associate in what is presently the Analysis Unit, Engineering Computing and Analysis Staff, 
Transport Division, Boeing Airplane Company, and supported by Boeing until the major 
results were obtained. 
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P, can take measurements on bearing angle. Bearing angle measurements are 
assumed to be normally distributed, with mean value equal to the true bearing angle, 
and standard deviation of o radians. (In the Boeing system o = 0.00087. This is 
an attempt to realize mathematically the condition that maximum error in bearing 
angle measurement is supposed to be 0.15 degrees = 0.00261 radians.) 

All measurement errors are (mutually) statistically independent. 


6. If the antenna of Pj faces Po, and if the range does not exceed some given range S, 


the probability of P; detecting Po is unity. (The results of this assumption can 
easily be modified to handle the case where "Probability of Detection" versus "Range" 
is known.) 

There is one complete antenna revolution every t, seconds (in the Boeing system 

tp = 4), with one measurement taken on each antenna revolution. 

In the event of a warning, a maneuver to shift the miss distance by 5 feet is to be 
executed (in the Boeing system 5 = 1000), at least ty, seconds before the impending 
collision, ty, chosen so this maneuver can be executed, (in the Boeing system 

ty = 12.5). ; 

Somé discussion of these assumptions is in order. While it is true that even if it is 
attempted, airplanes cannot fly a perfect constant velocity course, it is felt that the deviations 
from constant velocity can be held well below the error in measurement. 

The reasons for not using a simple reflection radar (i.e., one based on the reflection of 
transmitted waves, much like a flashlight) are: 

Since the transmitted wave need only be strong enough to be detected by the oncoming 
airplane (rather than so strong that its reflection from the oncoming airplane be 
detectable at the transmitter), the range in which detection is likely is increased. 
Measurement error is greatly reduced, since now all measurements are from antenna 
to antenna. 

The maximum error in distance measurement is derived from the error in determining 
the time when the signal from the responder is first received. It takes a known finite time for 
a pulse to rise to its maximum value. It is assumed that the pulse will surely be detected by 
the time it has reached this maximum value, and will not be detected before it starts its rise. 
For this reason, the error in measuring time of arrival cannot exceed this known "rise time." 
From the maximum error in measurement of time of arrival, and knowledge of the speed of 
radio waves, it is simple to determine the maximum error in distance measurement. 

The normality assumption, concerning the bearing-angle measurements has some basis 
in theory (Central Limit Theorem), but is open to serious question. However, results based 
on "normal theory" often hold up well when this assumption is violated. Further, this theory 
is mathematically tractable. In any case, before the proposed scheme should be adopted, it 
would be worthwhile to investigate its performance under deviations from the "normality" 
assumption, possibly using a Monte Carlo approach. 
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DISCUSSION OF RESULTS 

It is desired that an alarm be sounded in P, if Py will approach within some predeter- 
mined distance, | do| , of Pj, and that an alarm be sounded at least ty, seconds before such an 
occurrence. If Po will never approach within | dg] of P; it is desired that no alarm be sounded. 
The distance | dp| is chosen so that if the minimum miss distance exceeds it, a collision is 
impossible, taking into account the finite size of the airplanes. The technique suggested here 
allows (to within high accuracy) the probability of sounding an alarm when minimum miss 
distance is | dp| to be set at any desired level 1 - a. 

Because of errors in measurement, it is often impossible to achieve the desired results. 
The achievable results are computable as a one parameter family of graphs (functions) of 
"Probability of Sounding an Alarm" versus "Minimum Miss Distance."" The parameter is | v|, 
the magnitude of the relative velocity. 

If the assumption of normality of bearing-angle measurements is reasonable, the theory 
developed here is a refinement of that of Pina and Landau [1]. It yields improvement where 
most needed, at the expense of requiring a more complicated computer and greater storage 
capability. Compared to the Pina-Landau results, this procedure improves the ability to 
discriminate between collision and non collision courses. (Note that maneuvers, in the case 
of a course far from collision, can be quite dangerous.) Such improvement in the performance 
curve is to be expected, since this procedure assumes normality, while the Pina-Landau 
procedure leaves the distribution unspecified, and puts only a restriction on its range. 


NOTATIC’ 

Capi.al letters with a “ on top of them denote random variables. 

The corresponding capital letters (without a “ ) denote the values taken on by the 
random variables. 

The corresponding small letters denote "true" values. 

Thus, if R is a random variable corresponding to range measurement, and R is one 
corresponding to bearing angle measurement, particular observed values of R and g are 
denoted by R, X. The true range and bearing angle are denoted by r, x. From numbers 3 
and 4 in the section on Assumptions, we conclude that 


(1) P {|R-rj< og}=1 
and 


(2) E(R) = x. 


As noted in number 4 in the section on Assumptions, g is assumed to be normally distributed 
with standard deviation 0. 


THEORY NEGLECTING ERRORS IN RANGE MEASUREMENT 

All of this analysis will be done relative to P,, that is with P,; considered as a fixed 
point, and it will be a two-dimensional analysis. This is reasonable, since though it is azimuth 
and elevation that may be measured, these can be converted to a "measurement" of bearing 
angle. For details, see Pina-Landau [1], pp. 796-798. 
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We shall postpone until a later section (The Sampling Plan — Which Measurement to 
Store and When to Cease Taking Samples) the questions of when to cease taking observations 


e di 
and which observations are to be stored for use in making the final decision, and make : 
Assumption I _— 

We shall assume that we are given the measurements Rj, ...,R, on range, and the with | 
corresponding measurements X;,...,X, on bearing angle, which are to be used in making the for tl 
final decision, and that these measurements are ordered in time. We shall further assume that 
for given ranges Pyoceeoly and first bearing angle X4> (i-e., given the ouuer in wates specify Assu: 
completely the true geometric situation with respect to the manent frame), a... wer R,) and 
K,» we yk ,) are statistically independent of each other, and K,, — &) are mutually statistic- Pye: 
ally ciliates. be de 

We shall indicate in the section dealing with the sampling plan how to satisfy Assumption secti 
I. Assumption I simply means that if rj,... Ths Xy ave given, the knowledge of Rj,...,Rn meas 
does not change the statistical description of X,,... »Xn; and knowledge of any of the Xj does on th 
not alter the statistical description of any of the other Xj: Clearly, we must be given x, for 
the above assumption to be reasonable, due to the arbitrariness of the reference frame, for indep 
without Xx given, clearly the R, could not be independent. As is seen, the performance of the Thus 
test to be developed in no way depends on xj. 

It is intuitively seen that the only relevant information provided by the X; is given by (6) 

A A 
(3) ¥,=%,1-%) i=1,..., nel 
has @ 


(Y; being the measured bearing angle changes), since bearing angle is measured relative to an 
arbitrary reference frame. This intuitive reasoning may be further justified to believers in 
Likelihood Ratio Tests, by showing the test about to be derived, to be the Likelihood Ratio Test (7) 
of the hypothesis d = 0 against the hypothesis d # 0, based on X1,-++,Xp, where 
(4) |d| = minimum | miss distance|. 
Assuming rj >rj, 4, an assumption justified in the section dealing with the Sampling Plan, it and ¢ 
can easily be seen from paragraph 1 of the section on Assumptions, that if three different 
ranges, together with their corresponding times, are known perfectly, then | d| can be 
computed perfectly. Thus, we might consider a collision warning system based on range and (8) 
time measurements only. However, as shown in the following illustration, a change in | d| 
from 0 to 1000 feet, only alters those ranges in excess of 15,000 feet by at most 9 feet, for 
given 

P; ; .. at t=0 














first range of 60,000 feet. Thus, an error in range measurement of 9 feet could alter the Ther 
computed value of | d| by 1000 feet, making a collision warning system based on range measure- 


ments alone quite useless under number 3 of the section on Assumptions. 
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Thus, it appears that measured range values alone furnish us with practically none of 





the desired information. However, from 








(5) Y; = Xj - Xj+1 = arc sin d - em eS a 1 «ae, 
7 Ti+1 rj Tisd Tj 


we see that for r; - rj, 4 Sufficiently large, a knowledge of the bearing angle change together 
with knowledge of the corresponding ranges should enable accurate computation of d. Thus 


for the purposes of this section, we shall make Assumption II. 
iat 


> ify Assumption II 

d Tyo+++sTp are known perfectly. A collision warning system based on Xj,... ,X,, using 

ic- Tyoeeeeo lp only to find d for the measured Xj,---,Xp,, using the approximation in Eq. (5) will 
be derived. We shall later determine the effect of using the approximation (5). In the next 

on section we Shall relax Assumption II, i.e., in the test procedure replace T'y,+++,Tp by the 
measured ranges R;,..-,Rnp, and we shall show that such a modification has very little effect 
on the results derived in this section. 

Under Assumption I, for given|v|, rj,--- Tp» X,, the % , i=1,...,n are statistically 

independent normal (Gaussian) random variables with mean values x; and common variance o2, 


Thus, the random vector is 


Yy 


> 
iT] 


(6) 
Yn- 1 


has a multivariate normal distribution with mean vector, assuming the approximation (5), 











5 ve | 
r 3 
st (1) E(M=al 7: 1 Jade 
- _ 1 
rh Pn-l 
and covariance matrix = ” 
20 -0 0 
-0 20 -o” ‘ 
(8) a(n-1) x (n-1) - 0 a 902 . ; 
0 0 -o2 202 
i.e., 
rt=0 
— 202 if i=j, 
(A); = -o? if j=i+1, j =i-1, 
0 otherwise. 
Therefore, the density £$( ~ 3@) of Y is given by 
re- 
1 -1 
*< (A-de)' A “(A-de) | 
(9) g(r 3d) = He 
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and the constant ~ can be chosen independent of d so that 
2) 
ous {5 (A;d) dA=l. 
is 


With a slight bit of manipulation, we see that 


-3 wana + [2a rate aerate] = g(A) + h[z(A);d], 
-e 


(10) {(A; d) = ve 


where 
(11) 2(a) = vAntc. 


Making use of the Neyman factorization theorem, we see that the random variable 2(Y) is 
sufficient for d(Lehmann | 2].) What this means is that all of the information present in 9 


concerning d is contained in the random variable 2(9). 


Since 
i(n-j) if i<j, 


i(n-i) if i>, 


we see that 
n-1l i 


1% = 2, 


i=1 


(13) 


Note that ay > 0 if Ti41 < rj (see the section on the sampling plan, where this is justified). 


Since any continuous one-one function of a sufficient statistic for d is sufficient, it is 


preferable to work with 





with 
agains 


asser' 


(15) 


(Note t 
procedi 
most px 
Will be 
approxi 
4much 


To 
there aj 


vith ¢ 

















STATISTICS OF COLLISION WARNING SYSTEM DESIGN 











which is a normal random variable with mean value proportional to d, (assuming the approxi- 


mation in Eq. 5), and unit variance. This follows, since for any Cc; 


(2 n-2 
2 | 2 
Var c:%, = 20 », Cc. - 7. C.Ci.y >0 
sg | \ai i ii i~i+ 
if not all cj are 0, since 
n-1 n-2 n-1 
2 1 2 
i=] i= i=0 


with cg =c, = 0. Thus, aj > 0. A reasonable test of the hypothesis «that |d|< |d,| 
against the hypothesis #, that|d| > |d,|, of level @ is given by the procedure: 


assert that i, is true if and only if 


n-l 
(15) », aj Y; Se. 
i=l 


where cy is chosen so that 


_ 


n-e- 


Prdj=Ido} (2, 4% 


i=1 


=l-@. 


—K> 
[A 
fe) 


(Note that cg is dependent on |do|, though this is not explicit in the notation.) This test 
procedure is optimal, under the approximations in Eq. (5), in the sense of being a uniformly 
most powerful unbiased test of x , against Xj of level a (Lehmann) | 2], but this optimality 

will be of little concern to us here, and we will try to justify the whole procedure (including our 
approximations) by comparing this test's performance with that of certain tests which are optimal in 


amuch more desirable sense to be mentioned soon. 


To accomplish what is desired, let us examine the simpler situation in which we assume that 


there are only two possible miss distances |4,| and |4,| > |do|, 


i A 
vith Y having density i. when d= la.|, (s = 0,1). 
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To test the hypothesis 1 5 : d=|d,| against the hypothesis ¥ | : d= |d,|, we make use 


of a theorem of Neyman and Pearson (Lehmann, | 2]) which states that under certain mild 
mathematical restrictions (which are here satisfied), there is a class of tests of a simple 


hypothesis W > (i.e., a hypothesis which uniquely specifies a density) against the simple 


hypothesis 4 *, for which the one satisfying 
OO 
(l') P , {Reject Xo}=a 
Ho 


is the unique test minimizing 


(21) P x { Reject x F} . 


of all tests satisfying (1'). 


No test of size a can hope to do better in deciding between ue and . . We note that if 


for each a e€ (0,1) there is a Neyman-Pearson test, then since any test T of ue versus MT 


must be of some size, there is always a Neyman- Pearson test which does better than T unless 


T is a Neyman-Pearson test. 


’ 


The Neyman- Pearson Lemma is applicable to the problem of testing xO against N 1 


because each of these hypotheses specifies one unique density for y ° 


The Neyman- Pearson tests of 14 against 1 are those for which we accept ‘4 if and only if 


Se 
<k, 
fo(¥y,--+»¥p) = 





where k is some nonnegative real number. Under the assumption of normality of bearing 


measurements the density of ? under s{ is normal, with mean vector 


Yi,1 


Vi, n-1 
and covariance matrix A as given in Eq. (8). Omitting the rather detailed computations 


involved, we can show that in this case the Neyman- Pearson tests of {i versus i 1 are those 


for which we accept Ko if and only if 


n-1 i n-1 n-1 

(16) * | rit) Do ily, 5-¥0,j) + 1 2, -i04,5-%,;) | Y= Dd. vpy Sk, 
i=1 j=l j=i+1 i=1 

where k' is some real number. 












we ¢ 


whic 
exac 


| 4, 
Lett: 


(17) 


and « 
Eq. 


Eq. 
(18) 


wher 


(this 
dy, s 


comp 
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Again, rather than work with 


i=l 


we choose to use the random variable, 


a 


Re n- 


b, 2 i ‘, 
i=l n-1 9 n-2 
Y20 >. b’ 4 7 b* bay 
j=1 j=l 


}= 


io” 
* 








_ 


i= 


which is a normal random variable with unit variance. 
We now show that under the exact part of Eq. (5) the test given in Eq. (15) does almost 


exactly as well as the Neyman-Pearson test in distinguishing between certain crucial values of 








if | d, | and |4, | in those cases where the Neyman-Pearson test does well in such distinguishing. 
‘ Letting 
SS 1 na i ld | ld | 
ai? = z (-1) (n-i) z j {| are sin S' -arcsin —8 
1 Tia] Yr; 
- ‘ i+ d/ 
s=0 j=l 
n-l d d \ 
(17) +i Z. (n-j) { are sin | s| - arc sin |4s | 
*j+1 rj , 
nly if j=it+l 


* 
and defining ai in the same manner as b; was defined from b; after Eq. (16), then from 
Eq. (16) the Neyman- Pearson test of size a of N 4 versus ? assuming the exact formulas of 


Eq. (5) is to accept HG if and only if 


(18) 2. a ¥, < cy, 


where cy is chosen so that 


l-a 


ae) 
ae 
ry) 
eat 
> 
me 
1A 
re) 
R- 
i} 


yse i=1 


(this latter being computed under the exact part of Eq. (5). Note that cy depends on d, and 
qi, since ai depends on them. The desired result will be inncediaiamel in two sihien 
complicated steps, all computations being done under the exact formulas of Eq. (5). Inthe 
first step, we show for the crucial values of Iq, | and |4, | that the test given in Eq. (15) does 


almost as well in distinguishing between | d,| and |4,| as the test: 
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n-1 va 
accept KG if and only if ?. a Y; $ Cy? 
i=1 
n-1 A at 
(19) me? bi a, ¥, Se, ¢ =1-a. 
i=1 


In the second step, we show for the crucial values of |d,| and | 4, | that this latter test 


does almost as well in distinguishing between |4,| and | 4, | as the Neyman-Pearson test Eq. (18), 


In order to find the crucial values of \d,| and | dy], we assume that the sum of the wingspans 


of P, and Pz is 2w feet. Since the maneuver in case of impending collision involves a shift in 


miss distance of 1000 feet, it is vital that a useful test distinguish well between 
|d,| =© and |dy| = 1000 - ». 


For modern high-speed aircraft we may assume 100 < »< 300. Thus, the crucial values 


satisfy 
(20) | 4,| <9|q,| 
and 
7 
(21) 3 do] < [41] - 


Since it is clear that a miss distance greater than d > 0 has a higher probability of being 


interpreted as a noncollision than d, we will restrict our attention to the cases 
(22) 0 < d < 1000. 


(Completely similar results hold if we examine d = - | 4o| versus d= - |4,| <- |d 
and restrict to - 1000 < d < 0.) 


ol 


Also, since we are primarily interested in high relative speed cases, due to the required 
12.5-second warning time, we restrict consideration to the cases 


(23) r; > 11,000. 


We note that for 0<x<1 


xe 
(24) arc sin x = x ++ ?: o, xt, 0< eo <1. 


(See Courant, [3], Vol. 1, p. 408.) Thus, if 1>x>y>0 













(25 


Nov 


(26, 


Now 


we | 


(27) 


with 


chan 


and | 


(28) 


[dy], 


isa 
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(x3+ 2k _ y3+2k ) 


y 





‘ , 1 
x-y < arc sin x - arc Ssiny<x-y+¢% 


i 
° 


x2 + Ay + y2 - x2y2 


(1-x) (1-y2) 





1 
= (x-y) |1 + j 


Now making use of Eqs. (22), (23), and (25), we find 

















18). (26) ( a + < are sin - arc sin 4 < ( S. .S.) s.008. 
oS id i+1 a | 3) ed 
_ Now, letting 
Y= ° ’ 
a 
we let 
n-1 % . 
(27) g(y ;d) = Ey Y, 
aT Va n-1 9 n-2 
20 2 Y; mie Yj Yj41 


j=l j=1 


with the expectation taken under the assumption of the exact formula for bearing angle 
change, i.e., 


- 
A — — 
Eq(Y;) = arc sin —— - arcsin — , 
Ti+] rj 
and let h(y;d) be the same expectation under the assumption, 
A 
1 1 
(28) BRAY, = d{[—— - —}. 
ai Ti+] rj 
d 


In order for the Neyman- Pearson test Eq. (18) to do well for the crucial values of | do | and 


|dy|, since 


i=1 
is a normal random variable with unit variance, it is necessary that 
(29) (0's | 41) . g(a’; [4g |) > 5 


(where a' has the obvious meaning), as a glance at the standard normal distribution will 
show. 
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Rather indirectly (due to the extreme difficulty of working directly with the ay and aj), we 
shall prove a similar inequality for g(a; |4,| ) - g(a; |do|)-. 


Noting that the values of g and h are independent of magnification of y, we see from 
the stated optimality of the Neyman- Pearson test, and the fact that both 
n-1 mn n-1 A 
z. a, Y; and > a Y; 
i=1 i=1 


are normal random variables with unit variance, that 


(30) ; maximizes g(y; |di]) - gy; |do| ). 


maximizes h(y; |dy|) - h(y; |do|)- 


g(a; |dy]) - g(a; |do]) < efa’; |dy]) - efa's | dol) 
< 1.0042 h (a'; |dy|) - h(a’; |dg|) < 1.0074 [h (a'; |dy|) - h(a’; |do])] 


< 1.0074 [h (a; |dy|) - h(a; |d,])]< 1.0074 [ g(a; |dy]) - g(a; |dol)], 


g(a; |dy|) - g(a; |do]) < g(a’; |da|) - gat; |do]) < 1.0074 [ g(a; |dy]|) - g(a; |do|)] - 
Here, the first inequality follows from (30), the second inequality from (26), the third 


inequality from (28) (linearity of h in d) and (21), the fourth inequality from (30), and the 
fifth inequality from (27), (28), and (24), because the derivative w.r.t. d of h(a; d) - h(a; | do} ) 


does not exceed that of g(a; d) - g(a; |dp|) for d > | do| while these two expressions are equal 
for d= |do|. 


Hence, we have from (29) and (31): 
(32) g(a; | dy|) - ea; |do|) > 5/1.0074. 


Now, 


(33) g(a; |di| ) < 1.0042 h(a; |dy|) <9 + 1.0042 h(a; |d|) < 9.038 g(a; |do]), 


the first inequality following from (26), the second from (20) and the third from (26), all using 
(27) and (28). Now (33) and (32) imply 


(34) g(a; |dg|) > 5/8.1. 
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But Eq. (34) proves that the tests given by (15) and (19) are just about identical in distinguishing 
the crucial |d,| and [41 | values, as yet another glance at the standard normal distribution shows. 


From (31), we can also carry through the final portion of the "near optimality" proof, since 


the significance of (31) is that for the important |d, | the statistic 


n-1l 


>» a,%,. 


i=1 
can discriminate near |d, | almost as well as can the optimum statistic 
n-l 
rH as 2, . 
i=1 
To be precise, the size a test of x 4, (19), based on 


n-1 
yay Yat d=|do] + 1.0074 [|dy | - |do]], 
i=1 


has at least as great a probability of choosing 4 i as does the Neyman- Pearson test, (18), 


based on 
n-1 A 
>, aj ¥,at d= | d,| * Ss, 
i=1 
at worst there is a "loss of vision" of 0.0074 [ |d3| - | d| ] feet. This follows because from (31) 


g(a’; |dy|) - gia; [do] ) < 1.0074[ h(a; |dy|) - h(a; |dg|)], 
But, since h(a;d) - h(a; |do|) is linear in d - | d)|, at d = | do| + 1.0074[ |43| - | do], 
h(a;d) - h(a; |d,|) > e(a'; | dy| ) - g(a"; | do])- 
Hence, since from (24) the derivative of g(a;d) - g(a;|d,|) w.r.t. d is at least as great as 


that of h(a;d) - h(a; |d)|), these two expressions being equal for d= |d |, we have at 
d=d, + 1.0074 [ |dy| - }do| ] 
g(a;d) - g(a; |dg|)> h(asd) - h(a; | do| ) 


from which the desired result follows. It is shown in the following illustration 
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a 


4 
g(a;d) od g(a; |4,]) 
LZ. 





1.0074m > g(a’; |dy|) - g(a’; lag)4-7 h(a;d) = h(a} |do| ) 


m = h(a; |d,|) = h(a; |d,|) 








UL 


| do| |4)| d axis 


L | 


1.0074 [ [dy] - |dol | 














We can now prove that for all d satisfying |d)| <d<|d,|, with |d,| satisfying (21), 
g(a;d) - g(a; |d,|) is close to g(a';d) - g(a’; |d,|). (Note that a and a' are functions of 
| dol; |4| satisfying (21). This computation is done in three stages: . 

(a) We show g(a';d) - g(a’; |do |) close to h(a';d) - h(a’; | do| ). 


From (24) we see that the derivative of h(a';d) - h(a'; | do| )w.r.t. d does not exceed that 
of g(a';d) - g(a’; |d,]). Since these two expressions equal 0 at d= | do| » we conclude that 


(35) h(a';d) - h(a"; |do]) < efa'sd) - g(a’; [dg|) for d> [do]. 


From (24) we see also that the derivative of g(a';d) - g(a’; | do|) is increasing. 
Thus, since g(a';d) - g(a’; \a,[) is less than or equal to 1.0074{ h(a';d)-h(a'; | do])] 
for d= | dy | (by (31)) and these are equal for d = | do| , with the derivative of the latter being 
constant 
1.0074{ h(a';d) - h 


(35 do D1. g(a’;d) - g(a’; |do|) 






1 





| 4, | |44| 
it follows that 
(36) e(a';d) - g(a’; |d,|) < 1.0074 [h(a';d) - h(a’; |d,|)], for d > |do|. 
(b) We show h(a';d) - h(a’; |d,| ) close to h(a;d) - h(a; |d)|). 
From (31), g(a; | d |) - g(a;|d|) < 1.0074 | h(a’; | dy]) - h(a’; |do|)] 


< 1.0074 [ h(a; | dy|) - h(a; |d,])] < 1.0074 [ g(a; |dy|) - g(a; |do|)]. 
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But from (28), we see that we can write 
h(y; [dy] ) - hy; [dol )= [| 4y| - [dol] 7 (y)- 

Hence, 

(37) h(a';d) - h(a’; |d,] ) < h(a;d) - h(a; |do| ) < 1.0074[h(a';d) - h(a’; |do|)], 


for d2>|d,| . 


(c) We show h(a;d) - h(a; | do|) close to g(a;d) - g(a; |do| ). 
By the same reasoning as in (a), 
(38) h(a;d) - h(a; [Jd |) < g(asd) - g(a; |d,|) < 1.0074 [ h(a;d) - h(a; |d,| )]- 
Combining (35), (36), (37), (38), we conclude that 
(39) g(a;d) - g(a; | d,|) < (1.0074)?| e(a';d) - g(a"; |dg|)] < (1.0074)9| g(asd) - g(a; |do| )]. 


In a fashion similar to (31), (39) asserts that the properties of the tests based on 


are Similar for all d> |4,| : 


We summarize the results of this section. Based on certain approximations, we derived 


atestof x): |d| < |d,| versus 1): |a| > | do | which was shown to be almost as good in 
distinguishing ‘ _ «= |d, | versus i 1: d= |d) | as the best that can be done using no approxi- 


mations. Thus, the characteristics of this test can be used as a yardstick against which the 


: performance of any test can be compared, when the model is valid. 


In practice, of course, this test cannot be performed, since it assumed perfect knowledge 
of the ranges r;- In the next section, we shall show that the performance is not greatly altered 
when the r; are replaced by the available measurements Rj. That is, we present a performable 


test whose properties are close to the unachievable optimum. 


REPLACING THE SUGGESTED TEST BY A PERFORMABLE TEST 

The test procedure given by (15) is not a workable one, because the r; are not known 
perfectly. As was pointed out, a test based on the rj; alone would be extremely sensitive to 
errors in the knowledge of them. In this section, we shall investigate the effect of replacing 


the rj by the measured Rj in the test given by (15). Let 







(40) Resim) 5 i(q>-shi yp ool 
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: L R. R. 4 R. 
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A A 
A* A* 
i i+l 


if and only if 


(42) Ae 8 = © 


A A 
where A, is the value assumed by Aj when R. assumes the value Ri» and where, to be ex- 


plicit, for the given R;> C. is chosen from the relation 


n-l 
A 
(43) Pp f+ DAY. < Cc cine, 
a 1 ol a 
l= 


Aw 
with 6 a standard normal random variable and 


1 1 
alles“) 


(Note that in (43) we cannot replace = ’ ; by 2 ? as would be desired, since ¥. v the true 
bearing-angle changes when d = |4, 7 » are a Note also that C. is the value ae a 


veniam variable because it depends on the values A, and Y = of the neiees variables A. 


and q oi’) As an alternative to (42) if we rE. rejecting Ko when d = | 4 | a ieee 
as serious as accepting Ky when d = | 4, | we could accept x, if and only if 
n-l 
z: A. c, 
i 
i=1 


where for given Ri C is chosen from the relation 


n-1l 
< C>, with 








A n-1 
5 + >: Av. 
1 ol 

i=1 

















The test procedure actually proposed to handle the collision warning problem is to accept x ‘ 





Ino 


(44) 


and 


and 


(45) 


(46) 


whil 
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In order to compare the test of (42) with that of (15), we must compare the two test statistics 


MD-l pp A_A 
(44) (a) 7 A, Y,= Q, 
i=1 
and 
n-l A — 
(b) a A” 
t 2 
i=1 
and the two significance points 
(45) (a) Cy (or C) and 


(b) C (or c 


if (15) is modified as was suggested above for (42)). We begin by noting that for the same 


reason as in (30). ‘ 











1 
y=] =a maximizes 
an-1 
1 1 
n-1 Yily. yy, 
(46) h(y)= i+ 1 i 
i=1 n-1l 9 E n-2 
j=1 j=l 
while A 


=A maximizes 
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We now assume that 


(48) ee ae se ee 


Recall that o% is the maximum possible error in range measurement (number 3 in section 
on Assumptions). We shall see in the next section that this assumption is satisfied. For ex- 
ample, if .* 50 we may, in storing observations, only make use of those observations for 
which re- Tia > 3700. This already must occur in high-speed cases where relative speed 
is 1900 fps and one antenna revolution occurs every 4 seconds. Should antenna revolution 
rate be increased, in order for our results to be of much use, the error in range measure- 
ment would have to be decreased accordingly. The discarding of certain measurements is 
justified in the lower-speed cases because the suggested test gives such excellent results 


there even with the information being discarded. 


We also assume that 
(49) <<, 


which naturally follows for a, s 50 due to the 12.5 seconds required warning time. 


In order to compare the test statistics Q and 4, we take note of the following facts: 
4 is a normal random variable with unit variance; given A; the value assumed by AQ is 
a normal random y vm with unit variance. We shall show that no matter what value A; is 
assumed by A, E(QA, = A, oe ere n-1) is "close" to E(4), and also that Cc, and c 
are "close" to each other. This will prove that the test given by (42) and that given by (15) 
behave in almost the same way. In fact, the technique to be given will permit a quantitative 


comparison of the behavior of these two tests. In the development, we shall be guided by the 
approximations in (5). 





Using these approximations, for given Tyrese Th 
A.A 
(50) E,(QIA, = A, i= 1,---,n-1)+ nia) a? 
20 
and 
ae 

(51) E (a) = h(a) —S 

V2 o 


Our first task is therefore to show that h(a)/h(A) is close to unity. This will be 


accomplished by first demonstrating that 





(52) a is close to unity, 










and tl 


(53) 


conce 


samp 


(54) 


It fol 


(55) 


and 


(56) 


Let 


(57) 


Then 


(58) 


and 


(59) 
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and then that 


H(A) 





(53) h(A) is close to unity. 







In proving (52) and (53), we shall make use of (48) and (49) as well as the statements 


concerning (46) and (47). 


First we note, assuming rs > Vind (an assumption justified in the section on the 


sampling plan), from number 3 in the section on Assumptions, (48), (14), (41) 























(54) ” Te - i > 0, _1 - 1 > 0, a. * 0, A, > 0. 
r, ¥. R, R. i i 
i+ 1 i i+1 i 
It follows, using number 3 in the section on Assumptions, assuming r, > rae that 
ae | 
R. R. 20 
(55) i+ 1 i eis Ce) 1 
r. ~ ¥, 20 
.f i” iva) (,. =| 
“oa 5 Ti+ 
and 
- eee 
—. 20 
(56) 5 Tn i a foe a awe 
R. - R. 20 
Ff i i+1 (1 fe) 
Riva m, Riv 
Let 
5 = mi -r = mil - 5 
(57) f me (Ti-1 ri); F a (Ry R, ) 


Then, using (54), we see that 


oa) 
| 
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Now, from (46), we know that 








1 1 1 1 
60 A. --—|\< a. -— 
wes ), Fae )* 2, a ae a 


and from (47) 


n-1 n-1l 
1 1 1 1 
(61) ) ale -t)s , } A(z -e): 
- i+1 i+1 i 











20a 
i ~ ae 
r n-1 n-l 
n 1 ae ; ae | 
(62) : erro A, tar . R | : ai ts r, 








\ ~ 
n-1l i+ —— n-l 
(63) mo Nee Semen 2 

i Y 


— 
bar 
~ 
+ 
_ 
~ 
™— 
Pes 
! 
nN 
Q 
So a 
~~ 
| 
_ 
> 
eke, 
wy 
>; | 
_ 
' 
| 
~ 
ae 


Equations (62) and (63) together with (48), (49) and (57) prove (52). 


Using (56) and (57), we have 








20 
n-l 1 + 2 %0 n-l 
1 1 | F ae. 1 
6 ) | es = : ) ) ; 2 yaa 
(64) a, r. y. ° 9 A, | R. 1 R. 
cl i+1 i (, 0 7 \ “i+ i 


Equations (64) and (58) prove (53). 


We have thus shown, combining (62) with (64) and (63) with (58), that 






















(65) 


Note 


is wit 


is wit 


sults 


(66) 


Thus, 


upper 
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The term, 


n-1 \ 


d >. A. | 3 - ; | dh(A), 
é ta ad |, 


A A 
is within 0. 42 percent of the expected value E(QIA = A). while 


is within 0.42 percent of E(4), under (23) for |d| < 1000, as seen in (26). 


Rather crude inequalities were used to obtain (65), and for any particular case the re- 


sults could be sharpened greatly. We just sketch how this might be accomplished. 


We note that 


(00) n-1 : 1 1 \ ? aid ‘ % : n-1 *h4 ~* 
) | Et oe x. | R R R. 
‘ i+] i/ n 1 fom i 
i=l i= 
Thus, for given Ty veey PL (and therefore given a.) we can compute, using (1), sharper 


upper and lower bounds for the left number. Thus, we can greatly improve (59) and thus (63). 


Writing 


n 
A* = - (n-i) ze. + i 

or on de 

é itt 


we can get upper and lower bounds for A* , and thus as above obtain upper and lower 


od La 


— 
it 
—_ 
oa) 
_ 
~ 


bounds for 
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“i 
i+1 + 


1 “e | 
i+ 1 R; 





Thus, (58) can be improved, and therefore (62). Since this also improves (64), we have shown 


that (62), (63), and (64), (58) are improvable. Thus, (65) can be sharpened. 


As an example of the results of the crude (65), if * 15,000, 1" 50, f = 3750, 


while if r = 24,000, o 
n oO 


in co: 





perce 


distay 


ee ae 
J “os 0) Or i+] j 


When d = |4,|> we see from (30) that with very little error we can take the expected value of minin 


obser 


measi 


not be 


lv] by 


Ri an 
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for given range measurements R to be 


is proportional to |, , the effect of using 


rather than 


in computing . is to change the miss distance at which P {acc Ho } = 1 - @ by at most the 
percent indicated in (68). Thus, if | 4,| = 200, *. = 50, R, = 10,000, F = 5000, the miss 


distance at which P {acc x ” } = 1 - @ is achieved is between 194 and 206 feet. 


To summarize, we have shown the effect of replacing the test of (15) by a performable 
one to be small, giving quantitative results which estimate the "loss of vision.'' Equations 
(65) and (68) give the cruder general results. In the discussion following (68), a sharper 
method for estimating the "loss of vision" in any particular case, was outlined. In any case, 
the loss of vision can be made small by keeping error-in-range measurements suitably small 
(how small, to be found from (65), (68)). 


THE SAMPLING PLAN-WHICH MEASUREMENTS TO STORE AND 
WHEN TO CEASE TAKING OBSERVATIONS 
From the original assumption that all measurement errors, prior to discarding any 
observations, are mutually independent it can be shown that any sampling scheme for deter- 
mining which measurements to use which is dependent on only the original observed range 
measurements, R1, ... Ri will satisfy Assumption I. Due to space limitations, this will 


1 
not be proved in this article. 


Here is one reasonable sampling plan. Store R!, X! only if * > oY Estimate 


1 


lv|by|v'| =(R! - R'\/7_, where R' - R! > 13,000, and 7_ is the time between measuring 
s/s Ss s 


1 1 
R and R’: (Note that since required warning range is at most 24,000 feet and detection 
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range at least 50,000 feet, this is feasible.) From |v"| we can obtain a better estimate of the 
required warning range (error in|v'|is approximately at most 2 o,/ ‘. for small |d| , which 
is fairly small) and thus obtain near the maximum number of usable observations. Any sub- 


set of Ry» cons x. for which successive observations are F feet apart will do for 


po Ri This restriction is so that the results of the preceding section are applicable. 


As can be seen from an examination of test performance, it is desirable to let R, = Ri 


‘ A . 
= R'.. We let n be the random variable whose value is the number of used observations, 


A m 


A . lee ; ; 
-y 1a be the random variables whose values denote the indices of those observations 


which are used. 


COMPUTATION OF THE OVER-ALL PERFORMANCE OF THE TEST 


We have seen that for any given r, Xy {v| » we can compute upper and lower bounds for 
the probability of sounding an alarm. What is desired is to compute the probability of sound- 
ing an alarm for given | v| and | a| (or good upper and lower bounds). In theory, this is a simple 


matter, and is done as follows: 


Give ry . ~~ xX} and the probability distribution of the R: from the sampling 


A A : 
plan we can compute the probability distribution of A, i 1A. From this and the upper 


po 
and lower bounds derived in an earlier section (replacing the Suggested Test by a Performable 
Test), we can compute upper and lower bounds for the probability of sounding an alarm for 


true ranges r', and first bearing angle xy) in the obvious way, i.e., we use the relation 


-> F°., XL {sound alarm 
m’ 1 { } 


3 


ver’ 


1A A ° 4 
P ,_, {sound alarm |n =n, i ae <= 
rsx) 1 


Now, given |v|and|d|, we can assume some plausible probability distribution for the 


random variables m, rj r' , x}, and proceeding just as above compute upper and lower 


Ae a 


bounds for P {sound an alarm}. 


Iv|> [al 
This being rather laborious, we can obtain quick estimates of the performance of this 
test by looking at the "'best"’ and "'worst"' cases. 
For |v| = 1900 fps, assuming a detection range of 46,000 feet, antenna rotation time 


of ty = 4 seconds, required warning time of 12.5 seconds, maximum error in range 





test's 
secol 
seric 
This 

wher 


two 
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measurement ~* 50 feet, and standard deviation of bearing-angle measurements 


g = 0.00087 radians, roughly we have 


Worst Case: ry = 46,500 = 38,900 = $1,300 = ‘ 


Best Case : r= 46, 800 Ty = 39, 200 r, = 31, 600 r= 24,000 = r 


Setting P{alarm}= 0.9995 for ld] = 200 we find approximately: 





P{alarm, worst case} P{alarm, best case} 








El 





0.9995 0.9995 








0.9454 0. 5848 


0.4616 


Ls 


0.0001 








0.0000+... 























REMARKS 


In low-speed cases, where there is good reason to believe the constant velocity 
assumption seriously violated, it seems reasonable to require the system to warn the pilot 


when the second airplane is close, say closer than 5000 feet. 


It is seen that an increase in bearing-angle measurement accuracy would improve the 
test's performance. Performance could also be improved if the number of observations per 
second were increased, while not permitting bearing-angle measurement accuracy to deteriorate 
seriously, and increasing the accuracy of range measurement correspondingly (see (65), (68)). 
This latter suggestion would seem the only feasible one in case of very high-speed aircraft, 
where with the given number of observations per second it might not even be possible to obtain 


two observations. 
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LOGNORMAL DISTRIBUTION AND MAINTAINABILITY 
IN SUPPORT SYSTEMS RESEARCH* 


R. L. Bovaird,! and H. I. Zagor! 


Technical Military Planning Operation - TEMPO 
General klectric Company 


Santa Barbara, California 





This article presents a method for predicting the downtime 
distribution for a new complex electronic equipment. A theoretical 
discussion is given which indicates that the downtime distribution for 
such an equipment can be represented by a lognormal function. Ex- 
perimental data are shown which establish that the two-parameter 
lognormal distribution function adequately describes the downtimes of 
a diverse group of electronic equipments. It is concluded that this 
analytical expression for the downtime distribution of complex elec- 
tronic equipments provides the basis for a quantitative and scientific 
approach to the measurement and prediction of maintainability. This 
approach takes account of equipment maintenance characteristics and 
maintenance personnel capabilities. 











INTRODUCTION 


Of increasing importance is the recognition of maintainability as a major character- 
istic of military equipment. If maintainability is to have practical meaning in, for example, 
such areas as equipment specifications, it must have the qualities of being capable of definition, 
measurement and prediction. It is to this problem of measurement and prediction that we 


address this article, starting first with a discussion of downtime as a measurement of main- 


tainability.! 


DOWNTIME AS A MEASURE OF MAINTAINABILITY 


A review of the literature dealing with the subject of maintainability shows that many 


different equipment characteristics and operational conditions are included under this heading. 
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There are also many different viewpoints toward maintainability factors. Asa result, a num- 
ber of different measures and indices of maintainability have been found useful in one applica- 
tion or another, depending on which output functions and managerial decisions are of primary 
importance in each application. 

For example, it can happen that the user of a system is interested in the fraction of 
demanded operation times that the system fails to perform satisfactorily. This is termed 
operational unavailability in this article and depends on the frequency of system failure (relia- 
bility), as well as on the duration of individual failures. In this case, an appropriate measure 
of maintainability is the downtime per failure which depends on both environmental parameters 
(e.g., the type of maintenance personnel, facilities, and spares available) and system param- 
eters (e.g., fault detection and isolation features and ease of replacement of failed items). 

To the user, maintainability is a function of the dollars spent for support operations. 
In the case just mentioned, at a given support cost the system with the least downtime per 
failure is the most maintainable. Also, the system with the lowest failure rate is the most 
reliable, and the system with the least total downtime in a given operational period is judged 
to be most supportable. 

The system designer may not find these measures to be the most useful, for he is pri- 
marily interested in design parameters which he can influence, such as thoroughness of test, 
fault isolation capability, replacability of units in the system, and the like. For example, one 
design for a given system may permit 90 percent of its failures to be repaired in the field by 
isolating and replacing the faulty subassembly or part whereas an alternate design permits 
this for only 50 percent of its faults. The rest of system faults require unit replacement. To 
many designers, the second system is the more maintainable of the two. Others may judge the 
first to be more maintainable. 

There are great difficulties in resolving conflicting views (such as those just men- 
tioned), regarding an appropriate measure of maintainability acceptable to all. However, 
these difficulties must be overcome if satisfactory contractual agreements, regarding main- 
tainability specifications, are to be written. In particular, an unequivocal measure of main- 
tainability which can be measured and predicted will have to be found if this goal is to be 
accomplished. 

One approach to this problem which appears promising is outlined in what follows. 

The suggested avenue of attack rests on the experimental observation that the distribution of 
downtimes per failure for a variety of electronic systems can be satisfactorily approximated 
by a lognormal distribution with zero mean. Further, the two parameters of this distribution 
appear to be correlated with the maintainability parameters of the system and of its environment. 

The lognormal distribution provides a vehicle for quantitatively relating the various 
significant maintainability parameters in a form which is operationally meaningful. At the 


same time, guidelines are provided for structuring maintainability studies and discriminating 


between maintainability parameters, measures, and indices. Before indicating how this might 


be accomplished, the basic rationale and some supporting experimental evidence are presented. 
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CATEGORIES OF DOWNTIME 
When it is observed that an equipment fails to perform its functions satisfactorily, all or 
part of it is taken out of operation to locate and correct the fault. The fault may be corrected by 
an adjustment or a part may be replaced with a spare. When it has been verified by appropriate 
test that the fault is corrected, the equipment is returned to service. It may be placed back in 
operation, or it may be placed in standby, depending on the operational conditions at that time. 
Downtime can thus be divided into the following categories of time: 
Recognition 
Location 
Correction 


Verification. 


In addition, part of downtime may be occupied with getting a spare or in waiting for tools or test 
equipment. The additional time required to do this, over and above the actual time required 
for the repair activities, is supply and administrative time. Further, the time between the 
actual occurrence of a fault in the equipment and its recognition is usually not known. 
Downtime may, for most practical purposes, thus be divided as follows: 
Locate time (time from recognition that a fault exists to its location) 
Correction time (time from location of a fault to its correction, including post 
replacement adjustments, less supply time) 
Supply time (time required to get a spare, over and above time required for actual 
repair activities) 
Verify time (time from supposed correction of fault to verification by appropriate 
test). 
It is to be noted that these categories of time refer primarily to a single functional group of 
equipment which is capable of only one mode of operation (unimoded). Whenever a failure of 
such equipment occurs, the equipment must be taken out of operation. If the equipment is 
multimoded it may not necessarily have to be taken out of operation; however, it can be ana- 


lyzed by reducing it to a set of unimoded functional groups. 3 


DOWNTIME DISTRIBUTION 
The total length of time an equipment is down for maintenance when a failure occurs 
varies statistically from one failure instance to another. This is also true for the previously 


noted individual time categories into which total downtime per failure is divided. 


In some systems this category of downtime contributes significantly to system unavailability, 
although it is extremely difficult to obtain quantitative data on the actual magnitude of these 
times. In such cases, one solution often adopted is to strive for maximum thoroughness of 
test in order to minimize this contribution to downtime. 


4 System which can operate at several performance levels is multimoded. It can be divided 
into groups of equipment according to the specific functions performed during operation. The 
failure of some of its functional groups causes the system to pass into one of its alternate 
operational modes. (Cf. [ 4] ). 
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Analysis of experimental data on downtimes shows that the observed range of variables 
can be fitted by a lognormal distribution. Figure 1 shows a typical fit of the downtime data by 
a lognormal distribution. The particular curve shown is for the APS-20E radar [ 2] but is 
representative of much of the data obtained[3]. A rationale which leads one to expect this 


result is discussed subsequently. 


NOTE: 
DATA GATHERED IN TWENTY-FOUR 
SYSTEMS IN TWO TEST SQUADRONS 


OBSERVED REPAIRABILITY FUNCTION 
m——-—- THEORETICAL REPAIRABILITY FUNCTION 


(@ = 0.5891, m = 1.0008 HRS) 


PROBABILITY THAT REPAIR ACTION 
1S COMPLETE AT TIME t¢ 





A 1 1 1 1 


1.0 2.0 5.0 10.0 
ACTIVE REPAIR TIME, t, - HOURS 





Figure 1 - Observed and theoretical repairability functions for 
AN/ APS-20E radar systems 


The lognormal distribution is one in which the logarithm of the distributed variable is 
normally distributed. The parameters which determine the function are the median value of 
the variable and the standard deviation of the logarithm of the variable. The frequency function 


can be written in the form 


(1) f{(t/m,o) d(t/m) = yy lt/m) exp| - 3 Gr int (t/m)} at, 


where = distributed variable, 
m = median value oft, 
and o = standard deviation of ln t(dispersion parameter). . 

The variable of the frequency function is seen in Eq. (1) to be t/m while 02 appears in 
the exponential factor. The parameter o thus determines the shape of the function while the 
parameter m determines the scale on its variable. The multiplying constant 1/(moV2z) normal- 
izes the frequency function so that its integral over the entire range of t/m equals unity. This 
4 


The density function of a variable x, 0 < x < ©, where ln x is normally distributed with mean [J 
and standard deviation 0, is 


f(x, uu, 0) dx = (1/0xV2m) exp [-(an x-1)*/207] in 
If 4 = 0 and x = t/m, this reduces to Eq. (1). Note that m may be any constant or parameter 
of proper dimension, It is convenient here to let m be the median value of t. 
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frequency function is plotted in Figure 2 for several values of 0. The corresponding distribu- 
tion functions are plotted in Figure 3. The lognormal frequency function (Figure 2) is seen to 

be asymmetrical with a long tail to the right and a sharp cutoff to the left. As the dispersion 
parameter o becomes smaller, the distribution becomes more symmetrical. It is seen in Figure 
3 that the distribution is very nearly symmetrical and normal for values of og less than about 0.2. 


The limit as o——-0 is a normal distribution function with a standard deviation of zero. 
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Figure 2 - Lognormal frequency function 


As o (Figure 3) increases, the distribution of observed values of the variable is dis- 
persed over a larger range although the median or half way value remains constant. The 
small values to the left of the median are relatively evenly distributed from 0 to m, while the 
large values (i.e. , those which lie to the right of the median), are distributed relatively evenly 
out to very large values (i.e., the tail of the frequency function becomes very long and uniform). 

It can be shown| 5 | that the arithmetic mean "'t" (expected value of the variable) is given 
by 
(2) me/2, 
where median, 

e base of natural logarithms (e = 2.718), 

and o0 = dispersion parameter. 
Inspection of Eq. (2) shows that for a given othe expected or average value of the variable is 
directly proportional to the median. Equation (2) also shows that, for a fixed median, the 
average value increases exponentially with the square of the dispersion parameter. For this 


reason it is important in a practical case to keep o as small as possible. 
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Figure 3 - Lognormal distribution function 


Lognormal Theory 

The lognormal distribution function is frequently observed in two general areas of 
phenomena. One of these is in connection with growing organisms; the other is in connection 
with partitioning processes (e.g., placing objects in logically derived categories). There is 
as yet no general theory which gives the necessary and sufficient conditions for the occurrence 
of a lognormal distribution. However, it has been shown that whenever the "law of proportion- 
ate effect'' rules a statistical variable, it is lognormally distributed [5|. The law of propor- 
tionate effect holds whenever the change in a variable at any step of a process is a random 
proportion of the previous value of the variable. 

It has been known for a long time that response times of animals and humans to simple 
physical stimuli or perceptions are frequently lognormal (the so-called Weber-Fechner law). 
More recently, experimental studies of the reaction time of humans to more complicated per- 
ceptual patterns involving some degree of learning indicate that these reaction times are also 
frequently lognormal. 

It has also been observed that ''when a large number of items is classified on some 
homogeneity principle, the variate defined as the number of items in a class is often approxi- 
mately lognormal. Examples of this phenomena. ..are the number of persons in a census 
population class, the number of Sino-Japanese characters in a lexicographical group, and the 
outlay by households on classes of commodities,''[5 ]. Another reported observation is that the 
number of spores per leaf of a tree is lognormally distributed [5]. According to the above 


statement regarding classification on a homogeneity principle, this fact is a consequence of 
the manner of counting spores rather than an attribute of the spores on a tree. 
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Horvath has noted that the service times at atoolcrib and ina library can be fit very wellbya 
lognormal distribution | 6]. The process of locating ademanded toolor book is a partition process 
characterized by logic, for a coordinate system is used to find the tool or book (aisle, shelf, binor 
placeon shelf). Inboth cases, a large number of items is classified on a homogeneity principle 
(location). This is also true to acertain extent for faults in an electronic equipment, fora fault is 
located to acertain level, say to within a component or module. 

These observations suggest very strongly that the human partition process of placing objects 
inlogical categories somehow leads to a lognormal distribution of the number of objects inacategory. 
The precise manner in which this takes place and the necessary and sufficient conditions for a logical 
process to produce a lognormal distribution have not been discovered. However, these considerations 
indicate that it is not surprising that maintenance downtimes for electronic systems are lognormally 
distributed, since location and correction ofa fault in such an equipment is to a large extent a logical 
process of placing things incategories. The categories are functional assemblies or components, 
and failures are categorized according to the failed component. 

Practical considerations indicate that the distribution of downtimes should be asymmetrical. 
Asarule, most of the maintenance personnel are familiar with the operation and repair of the most 
frequently occurring malfunctions in systems that have been inservice anumber of years. At 
relatively infrequent intervals, however, anunfamiliar trouble occurs or an inexperienced techni- 
cian is called onto repairthe system. Suchrepairs require a longer timethan most repairs. This 
type of maintenance situation exhibits a pattern of repair times whichcluster about a relatively low 
value, but in whichthe distribution shows atailtothe right. Thus it is unimoded and skewed to the 
right. One member of aclass of suchcurves, the lognormal distribution, describes this situation 
very well over the observed range of variate (i.e. , downtimes). 

It is to be noted that other distributions of this class (e.g., Pearson Type III or Chi Square) 
may describe the observed data equally well. However, these distributions do not arise from the 
partitioning process which is inherent in the maintenance activity. 

An interesting property of the downtime data was observed by Howard [7]. The lognormal 
distribution function fitting a particular sample of data was found by fitting a straight line to the plot 
against downtime on lognormal probability paper of the cumulative fraction of the downtimes not 
exceeding a given magnitude. These same data were also plotted on semilog paper anda curved line 
was obtained, showing that this function is not exponential. Were this curve a straight line, the func- 


tion would be exponential. By fitting a straight line tothe tail of this curve and removing from the 


sample the downtimes falling inthis exponentially fitted portion, a smaller sample was obtained (the 


remainder) for which the process was repeated. Inthis way, Howard [7] was able to approximate 

the lognormal function by a sum of exponential functions. Further, it was found that the distribution 
function remaining after each removal of an exponential function was lognormal, although it 
approached an exponential functionas more exponential distributions were removed from the sam- 
ple. Ingeneral, one would expect the sum of exponentials to give a Pearson Type III distribution [8] . 
This indicates that over the observed range of variable, the data can be equally well fitted by the 


Pearson III (i.e., Chi Square) distribution. 
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Howard [7] also noted that the downtimes which were grouped together by the above process 


fell more or less into definite physical categories (e.g., routine adjustments, normal replacements 
of frequently failing part, andthe like). The downtime distribution within each of these categories 
was, to avery good approximation, exponential, while the distribution obtained by placing them 
together was lognormal. 

There is a growing body of evidence which indicates that, in general, downtime distributions 
within the various categories are either lognormal or exponential; and the total downtime distribution 


is lognormal, provided the equipment is sufficiently complex. Otherwise it tends to be exponential. 


Experimental Data 

Observed downtime distributions for anumber of equipments have been reported. Howard 
|7| has plotted total downtime for the fire control system for two large surface-to-air-missiles and 
also total downtimes for a heavy surface radar. Horne |9| plotted active work times inthe aircraft 
and inthe repair shop for the AN/ ARC -27 communication sets of two different aircraft squadrons 
and also active work time inthe aircraft for the AN/APS-20E radar. Scott |3] plotted total downtime 
distributions for three surface shipboard equipments (AN/SRR-13A Radio Receiver, AN/SPA-8A 
Radar Repeater, AN/SPG-48FCS) and listed the distribution determining parameters for two 
other equipments (AN/URT-2 Radio Transmitter, TED-4 Radio Transmitter). Scott also showed 
the distribution of part procurement time, preparation time, fault locationtime, andfixtime. In 
these data the fault location time, fixtime, andtotal downtime always showed a lognormal distribu- 
tion, but the preparation time and part procurement time were not asconsistent. Table 1 lists the 
dispersion parameter o,the median downtime m, andthe average downtime for all the equipments 
for which data were available. 

It may be seen in Table 1 that the data fall into two groups according to whether © has alowor 
high value. The first group includes the three distributions of repair times obtained on board air- 
craft. The second group includes all the time distributions obtained on board ships and on the ground. 
The average O for the first group is 0. 60 while that for the second group is 1. 40. 

The difference between these two groups lies inthe typeof maintenance performed. Airborne 
equipment like the APS-20E and ARC -27 are primarily maintained in operating condition by means 
of module replacement, whereas ground equipment maintenance may also include the repair of faulty 
components and part replacements, rather than just module replacements. Also, thepriority of 
ground maintenance is not usually as great as is that of airborne maintenance, for there may be spare 
systems available on the ground for immediate operation. For relatively simple equipments like 
radios this isusuallythecase. The situation is different, however, for acomplex system likea 
surface-to-air missile fire control system. 

It is to be noticed that the TED-4transmitter has ao = 0.69 and m =0. 561 whilethe URT-2 
transmitter hasao =1.4landm = 1.21. These are similar equipments eventhough the 0's andm's 
are quite different. It turns out that the TED-4 is manually tuned and a relatively simple equipment 
while the URT-2 is pushbutton tuned and relatively complex. The latter vives a great deal of mainte- 


nance trouble. 
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The SPA-8A is a highly modularized system in which a failed tube, for example, is replaced 


by replacing the module containing the tube. This would lead one to expect that it should be in the 


upper group in Table 1. However, the equipment has an ineffective fault location system (locating 


the failed module is difficult). The maximum maintainability permitted by its modular construction 


is thus not obtained. 


TABLE 1 


Experimental Data on Lognormal Distribution 
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Heuristic Considerations 


The relatively small variation in co among the diverse electronic equipments (Table 1) canbe 


understood when it is noted that electronic equipments all use similar types of components and that 


the manner in which these components are connected and arranged to perform a function is generally 


similar, regardless ofthe function. Further, a complex equipment can be viewed as being a logically 


connected set of simple equipments or functional groups. The designers of such equipments use 


very similar engineering techniques and principles in designing the fault locationsystem. Usually, 


the initial fault recognition includes an identification of the failed functional group. Asa result, 


greater system complexity due to an increase in the number, but not complexity, of functional groups 


comprising the system does not generally increase the time to locate and correct a fault, but it does 
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increase the failure rate of the equipment. Ontheother hand, greater system complexity due to an 
increase in the internal complexity of the distinct functional groups can appreciably increase the 
time to locate and correct a fault as well as increase the failure rate. 

This can also be viewed from the human factors aspect. In general, military maintenance 
crews will exhibit similar characteristics in that each is composed of people with similar training in 
similar environments and required to do similar kinds of tasks. In addition, the individual members 
of the maintenance crews are subjected to a renewal process inthe natural personnel turnover inthe 
Armed Services. Thus, ina statistical sense after a period of time, one would expect the various 
crews who deal with similar types of equipments to yield similar maintenance time distributions. 
However, the median time to perform maintenance is not expected to be a fixed quantity. These dif- 
ferences inthe median values, m, arise inthat some equipments are inherently more difficult to 
maintainthanothers, eventhough, onthe average, similar techniques of application and tasks are 
required. Thus the maintenance personnel's qualitative performance willbe similar, but the actual 
times required to do these tasks will vary between the equipments. 


Another factor is personnel motivation or incentive. In an aircraft, there may be more 


urgency to effect maintenance. Thus, one would expect the medians to be different from those 


obtained on aircraft maintenance action as compared with those same equipments when they are 
having ground shop maintenance work done on them. In the case of the ARC-27 data, the dis- 
persion parameter, 0 , as well as the median, m, are different for aircraft and ground repair 
shop maintenance. 

The line of reasoning just considered leads to the conclusion that variations in o and m 
among electronic equipments indicate variations in approach to the design of failure detection 
systems, personnel skill levels, motivation, troubleshooting routines, and component replace- 
ment techniques. Thus, variations in 0 and m among various electronic systems do not usually 
indicate variations in type or complexity of the comprising equipments, except as these factors 
affect the engineer's ability to design good fault location systems and to arrange components 
physically so that they are easy to replace. 

The dispersion parameter, o, and the median m, of the lognormal distribution function 
describing downtimes of an electronic system are indicators which are sensitive to the main- 
tainability of the system. They are directly determined by the interaction between the main- 
tenance characteristics of the system and the maintenance skill and motivation of the mainte- 
nance crew in the given environment. The parameter, o, indicates the relative proportion of 
short, medium, and long downtimes (i.e. , the variability in maintenance performance), while 
the parameter, m, indicates the general level or magnitude of downtimes experienced (i.e. , 


the level of maintenance performance). 


PREDICTING DOWNTIMES 
The preceding data and discussion indicate a way to predict downtimes for an electronics 
system which has but one mode of operation. Although there are reasons for believing that 


multimoded system downtimes also follow the same distribution function, the present discussion 
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is restricted to unimoded systems. A unimoded system is either nonfailed in its entirety, or 
else it is failed and out of operation. The method for predicting the downtime distribution for 
such a system rests on the fact that the form of the distribution function is known to be closely 
approximated by a lognormal function. This function is completely determined by the two inde- 
pendent parameters o and m. These parameters were observed to be strongly correlated with 
maintenance factors which can be estimated for a given equipment. Knowledge of the way in 
which o and m vary with the factors which determine them is as yet fragmentary. This knowl- 
edge can be increased, nevertheless, by analyses of downtime data and from experiments 

using equipment and personnel under controlled conditions designed to provide data on specific 
equipment. 

Examination of the data shown in Table 1 and the curves of Figure 3 shows that 0 = 1.4 
is a value for the dispersion parameter which is approximately correct for a wide range of 
electronic systems. This value is applicable to a system when its maintenance includes fault 
location, equipment adjustments, and spares replacement at all equipment levels from indi- 
vidual parts to full assemblies of parts, components, subassemblies, and replaceable modules. 

Another wide range of electronics systems is approximately described by 9 = 0.6. This 
second class of systems is repaired primarily by module replacement, and fault location to the 
module level is a relatively simple process. Further, administrative and spares procurement 
delays are virtually eliminated in this type of maintenance as, for example, in inflight mainte- 
nance of airbrone systems during the time the equipment is required to operate. 

Estimation of the appropriate median downtime for a particular system is more difficult, 
as inspection of Table 1 shows. However, in the authors' opinion it would be relatively easy 
for experienced maintenance personnel to estimate this parameter approximately for a specific 
system, provided the expected operational maintenance environmental conditions are quite 
specifically described. The variation of this parameter with the various factors which deter- 


mine it can be discovered in any particular case through appropriate experimental measurements. 


POSSIBLE APPLICATIONS 

The observations of the preceding sections show that it is possible to study maintaina- 
bility in a quantitative way. An approach has been indicated. The approach places the meas- 
urement and prediction of equipment downtime at the center of the study of maintainability and 
relates both system parameters and operational parameters to this maintainability measure. 
This measure in turn is given operational significance through its effect on system availability 
for satisfactory functional performance. Downtime thus provides a means of deriving quantita- 
tive functional relationships between system characteristics, the operational environment, and 


Significant operational effects. A way by which this might be accomplished is now outlined. 


Maintainability and Availability 
The user of a system is interested in its operational availability (i.e., in the fraction of 
total demand time that the system is capable of performing required functions satisfactorily). 


This important operational measure of support effectiveness depends on the nature of user 
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demands, on system characteristics, and on the capabilities of the maintenance and logistic port 
support organizations. These combine to yield system operational availability in the way shown are 
in Figure 4. . dete 
Figure 4 shows that from the user's point of view, maintainability has two major ele- tain: 
ments: maintenance and logistics. Maintenance refers to fault location and removal activities of th 
whereas logistics refers to the supplying of material, particularly spare items required to opel 
support the maintenance activity. The time required for each of these two activities determines cept 
the time to return the system to service once it has failed. This is the duration of a failure the | 
which, together with the frequency of failure, determines the total downtime (and hence uptime) 
experienced by the system during a period of operational use. man 
mal 
expt 
cont 
AVAILABILITY Of p 
(UPTIME) ; 
mal 
[ ‘ and 
(T me TOF AILURE ) Cie 10 RETURN 
| [1 
tat vere | | MivetO.SE 
Figure 4 - Relation of maintainability [2] 
elements to availability 
[3] 
Maintainability Parameters 
The time to repair a failure (excluding supply time) depends strongly on the operational | 4| 
conditions at the place of repair (such as the number and capabilities of the maintenance per- 
sonnel); the tools, test equipment and manuals available to them and the nature of the mainte- 
nance facilities. The time to get a spare depends on the distance of the supply stores from the [ 5 | 
place of repair, on the administrative procedures required to obtain a spare, and on the nature 
of the logistics communication and transportation facilities. In addition to these operational [ 6 ] 
factors, repair and supply times are also strongly influenced by system characteristics such 
as built-in fault detection and location capabilities, and equipment packaging and installation [7 
for removal and placement of failed items. 
Parameter Investigation 
The system designer can influence system characteristics affecting support effective- [8 , 
ness whereas the system user can influence operational conditions important for system sup- 
[9 
The point of view in Figure 4 is developed in greater detail in[10]. See also [1] fora 
closely related treatment of this subject. 
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port. In addition to support parameters which can be controlled, there are many others which 
are uncontrollable and frequently random in nature. This constellation of support parameters 
determine the operational availability of the system by way of operational reliability and main- 
tainability. The approach to maintainability measurement described in the previous sections 
ofthis paper provides a means of investigating quantitatively the functional relationship between 
operational availability and the many parameters which determine it. Figure 4 shows a con- 
ceptual framework in which this might be done which relates the user parameters to those of 
the designer. 

One avenue of investigation which appears rewarding is a study to determine the precise 
manner in which the two parameters of the lognormal distribution ( oand m) are related to the 
maintenance factors which determine them. The result would permit accurate prediction of the 
expected downtime distribution for a new system and would provide a means of measuring and 
controlling the operational maintenance performance of an equipment and its maintenance crew. 
Of particular importance in such a study is the discovery and evaluation of the precise effect on 
maintenance performance of various human factors such as training, experience, motivation, 


and operational environment. 
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The problem has been considered by Isbell of determining the 
path that minimizes the maximum distance along the path to a line in 
the same plane whose distance from the starting point of the search is 
known, while its direction is unknown. We consider in this article 
the analogous problem for the search for a circle of known radius, 
of known distance from a starting point in its plane. It is further 
shown that Isbell's solution is a limiting case of the problem posed as 
the radius of the circle tends to infinity. The problem appears to be 
of some practical significance, since it is equivalent to that of search- 
ing for an object a given distance away which will be spotted when we 
get sufficiently close—that is, within a specific radius. 








INTRODUC TION 

Suppose that it is desired to trace a path 
from a point A (Figure 1) in such a way that it 
minimizes the maximum distance travelled along 
the path until a coplanar circle, with center °o and 
radius s, is reached, given only that the distance 
from A to Oo is r+S. 

A heuristic argument will be presented 
justifying the solution that is given below, although 
a rigorous proof could be demonstrated, which 
however would be lengthy and tedious. We shall 


present a couple of examples for different values 


following properties. 





*Manuscript received August 9, 1961. 





of r/s, and shall also demonstrate that as s> « 
the solution tends to Isbell's | 1| optimal path for Figure 1 
hitting a straight line a fixed distance r away. 


THE SOLUTION 
The path comprises the straight lines AP, PQ, the arc QTR of the circle fo} ; 
A, radius r, and the straight line RS, where the points P, Q, R, S and the lines satisfy the 





center 


Let OA be an arbitrary straight line of length r + s andlet O be the 


circle with center O and radius s. Then Pison ©; PQ isa tangent to QO; also, P andQ 
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are such that the tangent XPY to © at P has the property that the angles XPA, QPY are equal; 
finally R and S are such that the tangent at R to the circle 9' passes through S, which is on 
© and through O. 

Firstly, it should be noted that, if 7 = OAO is the angle that the line joining A to the 
center of the circle it is desired to reach subtends to AO, then (a) the path APQ covers the 
possibilities 0< 7 < OAQ, (b) the arc QTR covers the cases OAQ <n < OAR, where OAR 
m and (c) RS accounts for OAR <n < 27. 

Let us first justify the path RS. After leaving ©', it is necessary to hit © in 
order to cover the remaining possibilities for 7 <~ 27. It then follows that the path of short- 


est distance that leaves © ata point R and hits © atS is such that RS is the tangent to 
©' that passes through O. This is because (a) if we leave ©' before the prescribed R, 
then we must remain outside it in order not to miss possible positions of the circle we wish to 
find that would be missed if we went inside ©’ , and then with this restriction the path is 
obviously greater than the one prescribed, (b) if we leave the circle after R all possible paths 
are greater than RS, and (c) if we leave at R, RS is the shortest path to O 

Now consider the path APQ. Let us first show that, of all paths hitting © at P' and 
then proceeding to Q' on the circle O|, the one prescribed is best. Firstly, to minimize 
AP'Q', given Q', AP' and P'Q' must be straight lines. Now the ellipse through P' with foci 
A and Q' is such that S = AK + KQ' is a constant s' for K on the ellipse, S <s' for K within 
the ellipse, and S > s' for K outside it. Hence it follows that the ellipse with foci A and Q’ 
that touches © minimizes AP’ + P'Q’ for fixed Q’. If this ellipse touches at P' = P”, the 
well-known elliptic property that the angles the tangent subtends to the foci are equal implies 


tA tA 


that the angles X P A, Y P Q are equal, where X Y is the tangent at P commonto O_- and 


the ellipse. The optimal path also has the property that P Q must be a tangent to the circle 
°' Hence, the only path having these properties is that for which P" is P and Q' is Q. 
This completes the heuristic proof. It should be noted that, in effect, we have re- 


stricted ourselves to convex paths. 


EQUATIONS FOR COMPUTATION OF SOLUTIONS 


From Figure 1 it may be seen that the following equations are satisfied by the solution. 


a. 5 
(1) y Y =? 


(2) a+B=>-0, 


(3) AP=rsecy, 
and, from the triangle AOP, 





(4) _ Sin (8-y)_ __sin 6 
A r+S rsee } 
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Manipulation of these equations to eliminate a, y and 6, finally gives the equation for: 
(5) =(1 - sin 28)= 4 sin® B cos 8B (2 cos 8 -1), 


from which 8 may be easily computed for particular values of r and s by successive approx- 


imation. 


LIMITING CASES 
1. s/r—0 
As s/r— 0, it is obvious that the solution should tend to a path in which we proceed 
along the radius to the circle ©' until the perimeter is reached, and then go round the 
circle. Equations (1) to (4) confirm this, in thata, y and 90. It might incidentally be 
, noted that B—> 7/4. 
2. Ss 
As s—~, the problem tends to Isbell's problem of searching for a straight line. It 
should first be noted that @ satisfies the equation 


r 
(6) cos ¢ =——- 


Then, ¢—»7/2 as s—», and from Eq. (5), 2 cos 8-1— 0 or B-—+7/3, since the other solu- 
tions 8 = 0, 7/2, to Eq. (5) obviously do not give meaningful paths. The values (8, 0) = (7/3, 


1/2) correspond to Isbell's solution. 


-" TWO EXAMPLES 
Before computing the solution for two particular values of r/s, we should first give the 


equation for the maximum distance D travelled along the path; that is, the distance APQTRS: 


(7) D=r[secy + tany + 27-a@ -y -9 + tang]--s. 
1. r/s=1 
Solving the Eq. (1) to (7), we obtain 
9 = 60°, 
mn. B = 56° 36’, 
y = 3 af, 
a = 15° 59', 
and D = 6.8007s. 
2. r/s=3 
@ = 41°25, 
B = 53°50, 
y = 17°40, 
a = 8°29', 
D 21. 0721s. 
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It can incidentally be seen that the distances D are a little smaller than r + 27r, 7.2832s and 
21.8450s respectively, the distances corresponding to paths proceeding radially to the circle 


© ' and then around its perimeter. 
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MATERIAL LOGISTIC SUPPORT OF WEAPON SYSTEMS” 


D. W. Whelan, SC, USN 


BACKGROUND 

In recent years the experts on defense posture have spoken with ever increasing frequency on 
the subject of weaponsystems. Neophytes on the subject have been at once charmed and confused by 
theterm. It would appear that some of the experts participate in the confusion if the vagueness with 
which they sometimes describe a weapon system is any indication of the exactness of the definition. 

The evolution of Naval warfare has been from massive fire power toward sophisticated 
weapons which have the capability singly or in small salvos to seek out and destroy the target 
by increasing the kill probability of each warhead. The evolution is marked by the transition 
from thousands of rounds of heavy ammunition to hundreds of surface launched guided mis- 
siles; from millions of rounds of aircraft ammunition and tons of bombs to thousands of air 
launched missiles. 

The WorldWar Il concept of gaining the objective through saturation by volume fire- 
power has yielded to highly selective placement of highly accurate missiles. This is the era 
of measuring distances in miles and speed in mach numbers. This is the era of measuring in 
seconds the time to acquire targets, select warheads, firing order, firing, and reloading to 
repeat the cycle. Accomplishing this feat presumes high-speed electronic computation of the 
fire-control problem and presupposes the earlier imput of a myriad of instructions to the 
computer. 

These are the eclements of the weapon system. CNO defines a Navy Weapon System as 
the combination of a weapon or multiple of weapons and the equipment employed to bring the 
destruction of the weapon against the enemy. Navy weapon systems may be shipborne by 
surface ships or submarines or they may be airborne. By definition, within a Navy weapon 
system the weapon itself is associated with the detection, location and identification devices 
necessary for the employment of the weapon; the delivery units employed to launch the weapon 


and the control units used to exercise control of the weapon after it has been launched. 


DEFINING WEAPONS SYSTEMS 

It is vital to any consideration of the support of the weapon system of tomorrow's 
Navy to define definitely and precisely the equipments and the personnel who make up the 
weapon system. It is not necessary to contrive a single definition of a weapon system into 
which all weapons fit nicely. It appears best for the purposes of the Navy, to define each 
weapon system by listing all the related equipments which are required for the full and proper 
functioning of the weapon. Any definition in these terms requires some arbitrary decisions 


LD 
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as to the extent of the weapon system in order to preclude the possibility of the definition be- 





coming so monstrous that it defies handling by a closely enough knit organization. For ex- 




































v 
ample, a ship's electric, pneumatic, and hydraulic systems may support a number of weapons a 
as well as administrative and housekeeping functions in the ship. The service-type facilities cor 
are better considered systems in themselves and are better supported as such. But, again, the 
each weapon system must be clearly defined before any specific consideration can be given to we: 
the material support of the weapon system. 
WEAPON SYSTEMS SUPPORT ORGANIZATIONS ig 
When the decision was made to develop an intermediate-range ballistic missile (IRBM) “% 
capability for the Navy on a compressed time schedule, the task was assigned to the Bureau = 
of Ordnance Special Projects Office. A number of facts and proposed goals precipitated this “a 
decision. The benefit of such an organization lies in the fact that the organization has a sin- - 
gle purpose and the authority, staff, and funds for its realization. The task itself is ambi- 
tious enough to stagger the layman's imagination but the organization which was created for ~ 
the purpose of accomplishing it provides for the most effective coordination of effort for its - 
realization. - 
Other missile systems such as Regulus and Terrier have been developed and made - 
operational within the regular basic organizations of the Bureau of Ships and the Bureau of - 
Ordnance. That these have been noteworthy accomplishments, goes without saying. But that sil 
there appears to be a more efficient method for the accomplishment of desired goals in new ” 
weapon systems, is evidenced by the recommendation of the committee on organization of the - 
Navy that BUORD and BUAER be consolidated in a Bureau of Naval Weapons. The organi- 
zational pattern of the Bureau of Naval Weapons considers under the Assistant Chief for Pro- i 
gram Management the vertical organization requirements of the weapon systems. It follows an 
the old organization axiom that where two or more people or groups of people, at the same ion 
level, labor on separate subtasks which are contributing to the whole task, there must be a ™ 
common supervisor or supervisory group that provides the central direction for established ~~ 
projects. THI 
Just as the term weapon system has been the victim of several definitions, the con- 
cepts of logistic support organizations on a weapon system basis have been described and em- rath 
ployed in several ways. coul 
It has, for example, been the practice of the U. S. Air Force-~ when a missile system then 
graduates from the Research and Development stage to production—to assign logistic respon- nanc 
sibility to one of the several air material areas. These air material areas are designated vari 
as lead air material areas (LAMA's) for the weapon systems assigned to them for program test 
support. These activities become, in effect, a combined technical bureau and SDCP for their succ 
weapon systems. It is inherent in the concept of vertically integrated systems of logistic may 
support for each weapon system that there will be built in plural managers of many common tend 
items of supply and a certain amount of material cross hauling. Such a logistic system has 
built into it duplication among weapon systems on a functional assignment basis and material 
item basis. 
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In the Navy we consider this a disadvantage when it is contrasted with the compara- 



























tively clean assignments for material management among inventory managers by cognizance 
. assignment and separation of functions between technical bureaus and program supply demand 
control points. The apparent disadvantage to the Air Force, however, must be evaluated in 
the light of the offsetting advantage which centralizes logistic support responsibility for each 
weapon system in the hands of a clearly identified logistic manager (LAMA). 

The debate between vertical weapon system support management and an integrated 
supply system may be illustrated by the example of a recent provisioning of an FBM equip- 
ment. The provisioning decisions developed a list of repair parts in which some $900, 000 
would be assigned for inventory management of ordnance parts and about $50,000 for common 
electronic and general stores material. No actual tally was kept on the number of man hours 
required to accomplish the ordnance parts procurement and supply support required. How- 
ever, it is estimated conservatively that the effort in the provisioning activity per item to ob- 
tain supply support was twice that required to procure and establish each ordnance item. 

This observation is not intended as criticism of any inventory control point. It is intended 
only to illustrate the point that total adherence to the integrated supply system concept can 
divert the lion's share of the effort toward the decidedly less valuable portion of the inventory 
investment. There are, of course, long range, offsetting advantages in favor of segregating 
inventory management on a commodity basis. The relative value of these advantages at the 
outset of weapon system support of vital programs being worked against compressed time 
schedules are deserving of continuing consideration. 

The U. S. Army system of support for weapon systems more closely parallels that of 
the Navy than does the Air Force. The Army Missile Command, consisting of the Ballistic 
Missile Agency and the Rocket and Guided Missile Agency has what is comparable to the 
technical bureau and program support supply demand control point responsibilities. Supply 
support of common items is negotiated with the several Army Technical Services (Signal, 


Engineer, and the like). 


THE NEED FOR WEAPON SYSTEM MAINTENANCE POLICY 

The beginnings of material logistic support requirements of weapons systemis come 
n= rather carly in the development of weapons systems. Since no one has ever guaranteed or 
could guarantee that the equipments of a weapons system will be failure-proot throughout, 
m then the design criteria for the equipments mus! include some provisions for their mairt< 
nance. The level of maintenance may be a functicn of the availability of technical skills in 
various operating echelons. The levels of maintenance may be a function of capacities of 
test equipment and other physical facilities to isolate ‘he failure and to give an indication o} 
ir Successful repair. These things are indications of the possible range of repair parts which 
may be required to maintain the equipments of Naval Weapon Systems on combatant ships, 


tenders and repair ships, in Naval industrial activities »nd at manufacturer's plants. 
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It is necessary, therefore, in the design of the several equipments which constitute a 
weapon system to give consideration to the maintainability of the equipment. When the con- 
sideration of maintenance is careful and complete, the complete, clear maintenance policy is 


easily evolved and the selection of the range of repair parts is greatly facilitated. 


PROVISIONING DOCUMENTATION AND STANDARD PARTS 

There are several objectives in the preparation of weapon systems which, when pur- 
sued independently are at cross purposes with each other and inventory policies which have 
been established in the Navy Department and in the Department of Defense. These are: 

1. The enthusiastic proponents of complete provisioning documentation have tended 
in some services to deluge the equipment manufacturers with such ponderous repair parts 
paperwork that the manufacturers are rightfully asking the service customers to decide which 
they want more—paperwork or hardware. Obviously, in selecting repair parts for new equip- 
ment it is always desirable to have a few more facts before making a repair-part decision. 
Just as obviously, there is a point beyond which the pursuit of additional facts (or documenta- 
tion) is not justified for small incremental increase in decision effectiveness which is 
attained. 

Great strides have been made in recent years toward providing the necessary docu- 
ments for the orderly, timely, and complete provisioning of equipment repair parts. More, 
however, can and must be done to keep repair parts provisioning specifications in step with 
the rapidly changing technology of weapon system. Military specification P-21078A is one 
example of a provisioning documentation effort for a maximum of useful provisioning informa- 
tion at acceptable cost. 

2. A second perplexing problem lies in the fact that equipment designers and pro- 
ducers are under great pressure in these times to deliver workable products on a compressed 
time schedule so as to meet early installation dates in new construction and major conversion 
ships. The equipment designers are required by specifications to use "standard" parts 
wherever possible. The advantages of standardization are well known for the benefits to 
inventory managers through the reductions in inventory range. The designers’ cry is that they 
are compelled to design around existing components. Their position is that they cannot 
design tomorrow's weapon systems with yesterday's components. Obviously, these diverging 


view points can be reconciled only by compromise. 


PROGRAM SUPPORT BY A NAVY INVENTORY MANAGER 

The logistic support of Naval Weapons System makes mandatory a strong and effective 
system of program support and supply support supply demand control points. Program sup- 
port of all the principle equipments of a weapon system in a single SDCP is the only sound 
approach to the support of weapons systems from an integrated Navy Supply System. 


It must be recognized at the outset that weapon system managers are nct generally 


supply oriented and Navy Supply Systems conscious. This is not because the Navy Supply 


System is a thing of mystery and comprehensible only to people who draw their pay from an 


activity under the management control of BUSANDA. It is because the weapon system 
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managers devote their attention to the development and production of weapon systems which 
meet the operational requirements of CNO. This does and should require the full attention 

of the weapon systems manager. He generally has neither the time nor inclination to become 
intimately familiar with the mechanics of the integrated Navy Supply System. It matters little 
to him that there are such things as cognizance symbols, federal catalog requirements, and 
mandates against dual stocking of material in the Navy Supply System. He is, in a word, 
impatient with this relative trivia. 

This condition is either unalterable or so nearly so that to do other than live with it 
would require such a prodigious supply orientation program that it is unlikely to be accom- 
plished. Under the circumstances the wiser course is to concentrate upon developing a 
general confidence among weapon system managers for the Navy Supply System. This con- 
fidence is most easily realized through dealings with one responsive program support inventory 
manager of the Navy Supply System for each weapon system. 

It is then the special responsibility of this office (in this case, the program support 
SDCP) to receive, utilize ani disseminate all the necessary logistic, maintenance, and 
technical information for all the equipments of the weapon system to the rest of the supply 
system. 

How the program support supply demand control point accomplishes these difficult 
tasks is a matter in which the supply system will be forever trying to improve itself. Safe- 
guarding against dual stocking by good material identification, cataloging and supply support, 
maximum utilization of assets, and financing the inventory are all vexing problems in that 
their solutions take time. Time is a commodity that is always in short supply in trying to 
assure adequate range depth and quality of repair parts for a weapon system against a com- 
pressed schedule for an operational capability. 

It is to be recognized that, by the time the customer begins requisitioning repair 
parts for a weapon system the preponderance of requirement-determination decisions have 
already been made. In order to insure the availability of repair parts by the time the equip- 
ment is operational, initial delivery of on board repair parts and supply-systems stock must 
be concurrent with the equipment delivery. When this requirement is coupled with the re- 
quirement to buy repair parts for which there is no other requirement than the weapon sys- 
tem, it is obvious that repair parts must be ordered when the equipments go to production. 

At that time, there is no demand history from which to project range and depth requirements. 
Yet, it is at this time that all the requirements decisions must be made. There is no 


opportunity to test the validity of these decisions until the equipment goes into operation. 


WEAPON SYSTEM PROVISIONING 


Historically the initial range and depth of equipment repair parts has been established 


by a spares technician who tries to recall experience with similar parts and who supplements 


his experience with intuitive judgement as to what items and how many of each should be 


procured. He has been required to consider such factors as consumption rate and variability, 
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cost, weight and cube, and compensability and shortage cost. In truth, however, many re- 
pair part technicians do not even comprehend the meaning of these terms. 


Efforts are being made to set down an orderly pattern for repair parts decision mak- 


ing. The first step is a mathematical model which gives balanced attention to all the factors 


which influence depth of buy decisions in provisioning. This work is being conducted under 
Contract Nonr 2713(00). The ambitious goal of this program is a balanced load of repair 
parts which will fit within the available stowage cube of the Naval ship. The decisions as to 
what and how much material is to be carried in a combatant ship cannot be made logically 
without consideration of the type and frequency of resupply of deployed vessels. Obviously 
the economics of inventory management must be considered at this point. 

One fact remains unalterable in the treatment of any weapon system. That fact is 
that the task of requirements determination of repair parts for new equipments is a most 
difficult one with which to cope. A sound approach to new provisioning techniques is in the 
making with the development of parameters and the determination of parameter values. Some 
of these parameter values are known now or can be assigned for weapon systems and equip- 
ments within weapon systems. Unit price, reorder cost, holding costs are relatively straight 
forward. 

Penalty cost or the cost of equipment downtime is more nebulous. The absolute value 
of a sustained operational capability of a given weapon system in the scheme of national de- 
fense is not easily determined. Weapon system managers have several methods at their dis- 
posal to determine the probable inherent reliability of a weapon system. Whatever this in- 
herent reliability may be, it must be assumed that it may be brought up to a figure approaching 
100-percent reliability by the provision of the proper depth in the necessary range of repair 
parts. In substance, the proper maximum amount to invest in repair parts is the difference 
between the value of the offensive or deterrent capability to our defense position and the raw 
costs of development, equipment production, and manpower and facilities required for the 
system. Obviously, the investment in repair parts approaches infinity as the reliability 
approaches 100 percent. While 100-percent reliability and on-line availability is desirable 
as agoal, any sober consideration of the situation leads to the conclusion that something 
less than 100 percent is all that is attainable within the funds which may possibly be made 
available. Thus, it appears that military strategists must make the initial determination of 
the total value of each weapon system. After the best efforts have been exerted by the design 
and production people to deliver the system which has the best attainable inherent reliability, 
the difference remaining between that total cost and the total value of the weapon system 
should be the maximum amount available to operate and maintain the system. It is realistic 
in this manner to assign a fixed maximum dollar figure for the repair parts investment. 

It remains, however, that an accurate determination of projected requirements must 
be made. The most essential ingredient in this determination is a forecast of the useful life 
of each repairable and replaceable part in the equipment. As the complexity of equipments 


increases, it becomes increasingly important to design equipments for easy maintenance by 
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the service personnel. In order to hold equipment downtime to a minimum, self-checking 


features must be incorporated in the equipment which will isolate the failure to a specific 


part or component. Then this component must be designed for universal interchangeability 


with the replacement part or component. The trend in weapon-system equipments design is 
toward the replacement of subassemblies instead of piece parts. 

This trend complicates the determination of the required depth of replaceable parts. 
As the size and complexity of parts increases, the cost increases. In the past overestimating 
requirements for individual electronic piece parts, for example, was not relatively signifi- 
cant because of the comparatively low cost of parts and the high probability of the excess 
being consumed in other equipment applications. In today's context, the repair part isa 
group of 50 to 100 or more electronic piece parts assembled in a plug-in printed circuit 
module which probably costs between several hundred and several thousand dollars. 

Relatively little has been done to develop techniques for estimating probable failure 
rates for modules which are in most cases single equipment application items. To over- 
estimate parts requirements in today's equipments is to generate large stocks of high value 
material which must be scrapped when the equipment is superceded. The need is indicated, 
therefore, to do an accurate job of estimating requirements by determining the probable 
frequency of component failures in the terms of the periods which the equipment may be 
operating. Failure rates stated in terms of hours of operation fit the requirements of equip- 
ments which operates regularly or for extended periods. Fire control, launchers and test 
equipment fall into this category. Failure rates for components of the missile itself are best 
expressed as the expected failures for a given number of tests whenever the missile operates 
only periodically for test and check-out prior to firing. 

In the case of both flight hardware (missile) and capital equipment failure data, these 
data must be estimated until there are empirical data generated to replace them. These 
estimated rates must be developed from the known functional behavior of the bits and pieces 
of which the modules are manufactured. 

Printed-circuit electronic modules, for example, may be made up of 50 to 100 or 
more resistors, capacitors, diodes, transistors, and coils in various combinations. It is 
necessary, at the outset of the equipment production to develop arithmetic sums of weighted 
failure rates of the individual pieces to estimate the probable failure rate for the whole 
module. 

This technique is obviously time consuming and the cost of this approach is prohibitive 
for the lower cost modules. Its value exists for higher cost modules, particularly where the 
state of the art indicates that there is a good spread in the useful life of the several com- 
ponents. 

The value of this approach was demonstrated dramatically for a small number of very 
high cost components in the initial provisioning of Talos missile plug-in units. The con- 
tractor's first estimate of repair parts cost for USS GALVESTON (CLG-3), was approximately 


4 million. Failure rates were subsequently developed for the three highest cost items. The 
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reduced requirements for the three items reduced the total cost of all Talos components 
down to $0.9 million. 

The rank and file of equipment specialists who make the initial requirements deter- 
mination are inclined to reject this technique in favor of their own intuitive judgment. The 
need is indicated for an education program which leads to every-day use of this tool for 
formalizing the thought processes which go into first, the selection of the item and, second, 
the procurement depth of the item. 

The need is also indicated for developing techniques for distinguishing between con- 
sumable and repairable repair parts. This need is also an outgrowth of the trend toward 
plug-in units of repair parts for new weapons. The trend toward modular construction is 
accompanied by the trend toward miniaturization of piece parts and modules. There isa 
logical inclination toward the consideration of modules as repairable because the failure is of 
one or more of the piece parts in the module. The inclination is that if the faulty part or 
parts are replaced, the module is repaired. 

In many instances this 1s true. However, it is also true that the cost-in time and 
effort to isolate the fault and correct it oftentimes approaches the replacement cost of the 
whole module. There are investigations underway to develope trade-off problem criteria for 
distinguishing consumable items from repairable items. One such study is being conducted 
by TEMPO for the Special Projects Office. 

The criteria for this type-decision-rule contain many variables. One of the more 
perplexing of these is occasioned by the trend toward miniaturization, particularly in elec- 
tronic equipments. The present state of the art includes a capability to compress 100 of 
today's electronic bits and pieces into a space of 1l-cubic inch. It is quite likely that in 10 
years miniaturization shall have been advanced to the point that 100, 000 of today's pieces will 
be compressed into a space of 1-cubic inch. 

When the art has advanced that far, one thing will be quite clear. That is that there 
will be no point in designating the miniature module as repairable, since the whole thing will 
be less large than a lady's watch. Until that time, there is a need for "fix or replace" deci- 
sion rules which consider space, both work and storage, training for the crew, and piece 
parts inventory on the one hand and the space to stow replacement modules and their higher 
cost on the other. 

We have barely scratched the surface in this all-important area of provisioning. The 
Navy Supply System, including the technical Bureaus as principal item inventory managers, 
need, very much, a set of decision rules which will assure a consistently correct approach 


in this area which draws down so many material dollars. 


MAINTAINING LOGISTIC SUPPORT 
Up to this point, the discussion in this paper has been directed toward first-stage 


logistic support, that of the logistic planning and determining material requirements. It is 


clear that much of this work is accomplished on the assumptions that a certain set of operat- 


ing conditions will obtain, that operating forces will utilize the weapons system furnished 
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them in accordance with instructions furnished them, and that the weapons system will func- 
tion in accordance with engineering estimates. 

At the present time, most of the "'fact"' input to tomorrow's weapons is based upon 
the projections of trends established in yesterday's weapons. It is too much to assume that 
the logistic package will be entirely suitable and that the weapon systems will operate without 
further guidance and tailoring where necessary. The need is well recognized for a feedback 
of information on operating performance of the equipment itself and the logistic support 
system. 

The first need for such information is to correct deficiencies in the weapon system. 

Some of the deficiencies are of such magnitude that they will require engineering 
changes to the equipment. More of them can be eliminated by adjustment in logistic support 
(e.g., adjustments in range and depth of repair parts) and increased tempo of resupply. 

The normal mechanics of the replenishment requisitioning process, given enough 
time, will restore the balance in the logistic support equation which may be in imbalance as 
the result of unpredictable factors in the equation. When the consumption rate of repair 
parts is higher than that which was anticipated, more frequent ordering of the repair parts 
temporarily relieves the situation. However, without a feedback of information on the 
operating performance of the equipment involved (failure reports) the imbalance is not fully 
satisfied. Abnormal consumption of repair parts presents a two pronged problem. 

First, abnormal consumption of repair parts is an effect and does not, in many cases, 
pinpoint the cause. The flat tire and the blown fuse are effects which presuppose an earlier 
cause of tube failure and an overloaded circuit, respectively. Similarly abnormal mechani- 
cal failure in weapon equipment presupposes faulty lubrication, strain, or some other cause; 
and abnormal electronic part failure may be caused by a fault elsewhere in the circuitry. 
Without investigating the causes of equipment failure thoroughly, the continuous replacement 
of failed parts is as unproductive as bailing out the boat without repairing the leak. Failure 
reporting is the first step in the cycle which is not complete until the reporting is followed 
by investigation, analysis and correction. However, failure reporting is the very important 
first step. In the Naval Service, the operating failures occur many miles at sea, a long 
distance from the laboratories and factories from which the equipment is produced and where 
the investigations and analysis of operating failures will lead to satisfactory correction of 
the faults. 


If the failure reporting system fails, and many have, it is usually because the work 


falls to operating personnel who are already overburdened with paper work and who usually 


see no tangible result of their reporting effort particularly ''on their watch." It follows, 
therefore, that failure reporting systems must be simple. Reporting formats must be 
concise enough to be completed quickly without being meaningless. The failure-reporting 
program should continue only long enough to serve its purpose. It should be limited in scope 
to the reporting of significant failures. The efforts of failure reporters should be rewarded 


with a timely feedback of failure analysis and corrective instructions. 
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The second problem is related to the first. The Navy logistic system is based on 
onboard allowances of repair parts, replenishment stocks in tenders and supply ships, back 
up stocks in the several echelons of the supply system ashore and a procurement program to 
replenish the pipeline. Abnormal consumption on board ship causes abnormal demands 
through the whole system. The first unsatisfactory condition is a rash of not-in-stock posi- 
tions in the logistic pipeline at all levels. More frequent ordering of the parts in question 
provides partial relief. However, when the onboard allowance and tender/supply ship stock 
levels are not adjusted, the second unsatisfactory condition obtains, that of inadequate endur- 
ance. The need, therefore, is obvious for a feedback of usage information to effect an up- 
ward adjustment of stock levels afloat to compensate for the increased consumption. 

In the case where there is a feedback of failure data and usage data to the equipment 
engineers and supply program managers, respectively and independently, a more complicated 
interaction may develop. At the time the equipment engineers are working toward more 
reliable components to reduce equipment downtime, the supply system is usually procuring 
more of the failing parts in response to demand data. 

Again, it is obvious that there must not only be a failure-reporting system and a usage 
reporting data system, but the two systems must be coordinated. The supply program 
manager must know what the equipment engineer is planning and doing in order to avoid 
stocking the now obsolete part and to stock, as necessary, quantities at all levels of the 
superceding part. Where the correction of equipment deficiencies is formalized in an altera- 
tion program the supply manager must participate to the extent of deleting and disposing of 
superceded parts and providing the replacement part. There is no assurance of responsive 
support from the supply system if the alteration program is not formalized to the extent that 
the supply program manager participates in it. 

At this point it is desirable to review the process the supply system follows placing 
the initial repair-part support for equipment in ships and the means of sustaining support for 
equipments. The COSAL program includes an allowance parts list section (APL) in each 
COSAL. APL's for new equipments are usually developed through the provisioning process. 
Each ship is now, or will be in the near future, provided with a new COSAL incident to con- 
struction or major conversion. This COSAL will be revised at the time of periodic ship 
overhaul. 

The APL's for new equipments are, as has been discussed, the result of a series of 
best estimates of part failures and are the best estimates of what will be required to support 
the installed equipments. If the COSALS are to be kept in step with the actual requirements 
of the equipments they support, the APL's must be revised to reflect the repair parts usage 


experience which comes from a usage data program. This writer's knowledge of all inven- 


tory managers' usage data and APL maintenance programs is admittedly incomplete. How- 


ever, 40 months of more than casual inquiry in this area leads to one person's conclusion 
that while interested and alert to the problems, the inventory managers usually do not have 
the time or funds to do a comprehensive job of updating the range and depth of APLs for 
COSALS. 
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Excellent progress is being made in the area of supply availabilities at the time of 
ship overhaul. Armed with a new COSAL for the ship, the supply-availability team purges 
out repair parts which were applicable to equipment which has been removed. Stock records, 
bin tags, and locators are brought up to date to reflect stock number and cognizance changes. 
The big questions are several. Has the quality of the APL's in that COSAL been improved in 
range and depth to reflect the historical repair part demand data? Have quantities been in- 
creased where demands have been higher than the original COSAL estimate? Have quanti- 
ties been reduced where equipment part reliability is higher than originally anticipated? Are 
new items being added on the basis of demands which could not have been anticipated in the 
original COSAL? If the answers to these questions and a battery of others like them are 
negative, the overhauled ship goes back to sea with a repair part allowance which has been 
brightened up with new labels but without substantial qualitative improvement. 

The first emphasis on ship's allowances has been to ensure that all the repair parts 
called for in the allowance lists are on board. This has been done without substantial 
challenge to the validity of the lists themselves. Unquestionably, ships must carry neces- 
sary Spares. There are also two important reasons why ships should carry only those repair 
parts they will need and can use. The first is money. Ina climate of fixed budgets, invest- 
ments in excessive range and depth of unnecessary parts prevents the necessary investment 
in needed parts or limits the number of ships which can be so outfitted. The second reason 
isa space constraint. As ships and their weapon systems become more complex the number 
of required repair parts increases. As these new weapon systems are placed on more 
smaller ships (e.g., destroyer and submarine types), the task of finding places for these 
repair parts becomes more complex. The USS TRITON reportedly has repair parts stored 
in 900 locations in the ship. 

Weapon system reliability, dollar constraints, and space constraints are irrefutable 
arguments in favor of continued surveillance and modernization of Allowance Parts Lists. 

If we can accept as a fact that it is necessary to update APL's systematically based 
upon historical usage data, the problem is how and where to begin when it is well known that 
the inventory managers are already hard pressed to keep up with the job of writing new 
APL's for new equipment. It is manifestly impractical to expect that it is possible to update 
all existing APL's. The key to approaching the task lies in concentrating first on the high 
military worth equipments. 

Military essentiality studies have been instituted as a means of ensuring first that 
critical repair parts for critical equipments are always available. Military essentiality is 
subdivided into mission effect and compensability. Mission effect designations are the means 
of ranking the importance of equipments to the accomplishment of the ship's mission. It is 


submitted that, with certain modifications, it is entirely practical to use these mission 


effect designations as a means of determining the order in which weapon system APL’s should 


be updated. 
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NAVAL WEAPON SYSTEM LOGISTIC AGENCIES 

There is discussion earlier in this article of the general pattern of material logistic 
support of weapon systems. The most recent major departure from the general pattern of 
Naval weapon system support is in the Fleet Ballistic Missile System. The Chief of Naval 
Operations directed in March 1959 that CINCLANTFLT establish a Fleet Polaris Material 
Office. This office is charged with responsibilities for processing FBM requisitions, exer- 
cising issue control over critical FBM Material, recommending changes in allowance lists 
and maintaining liaison with supply activities and inventory managers. 

It should be noted that the PMO is a Fleet activity. Its job is to ensure that the fleet's 
requirements for the sustained support of the FBM weapons system are relayed to, under- 
stood by and satisfied by the shore establishment. The PMO is not charged with the func- 
tions of supporting the FBM weapons system. These are the functions of the Shore Estab- 
lishment. In a sense, the PMO is the Fleet voice in announcing the what, when and where of 


required support. 


The value of the PMO rests in the fact that its attention is directed entirely toward 


the FBM weapon system and toward the entire FBM weapon system. It is intended that the 
PMO will be expert in knowing all FBM requirements. One of its most effective tools for 
doing its job is a thorough understanding of the support contributions of all the inventory 
managers. It will receive the fleet support requirements from the operating FBM units and 
pass them in supply language to the highest levels necessary in each inventory manager's 
distributive system to satisfy these requirements. 

As the logistic agent ashore for the FBM fleet, the PMO must ensure that the logistic 
requirements of that fleet are understood by the shore establishment and satisfied as stated. 
It should not be inferred from the foregoing stated objectives of PMO that its establishment 
is a denial of the integrated Navy supply system. The establishment of PMO is recognition 
of several significant facts: first, that within weeks of the delivery of operational missiles 
the first SSBN will be deployed on operational patrol with its strategic target assignment; 
second, that once operational, each fleet unit will be deployed outside CONUS continuously 
until overhaul; and third, that the strategic value of the FBM weapon system varies directly 
with the capability to keep the weapon system "up" or operational at all times. 

The FBM weapon system is a most ambitious undertaking with the promise of being 
one of the most effective deterrent weapons 1n the United States military arsenal. The PMO 
is one of the devices which will contribute to the logistic assurance of this capability. The 
value and urgency of the FBM weapon system to this country's military posture justifies the 
establishment of the PMO. 

In accepting the concept of the PMO for the FBM weapon system and assuming the 
substantial contribution which is expected of it in the support of FBM ships, shouldn't the 
concept be extended to other Naval weapon systems? It is submitted that the answer is 
probably negative. The FBM program presents a pattern of considerations which are 


singularly adapted to the PMO concept (i. e., identical capabilities in all ships, single and 
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identical assignments, and similarity of installed equipments). New construction surface 


ships are being configured with several combinations of surface-to-air, surface-to-surface 
and surface-to-underwater weapon systems. Terrier, Talos, Tartar, ASROC and others are 
separate weapon systems which will be arranged in various combinations in CG, DDG, DLG, 
and CV type ships. The individual differences in ships themselves, the schedule for their 
completion and fleet assignments present a condition which is not conducive to improving 
fleet support by a PMO-type logistic agent. 

Surface ships which are commingled in various combinations for various periods in 
striking forces are probably susceptible of most effective operational support through the 
existing logistic channels. 

However, many of the support problems in initial outfitting and reoutfitting are pres- 
ent in surface ships which are being met for SSBNs through the PMO. The Ordnance Supply 
Office has recognized—for purposes of initial outfitting—that each ship, with its combination 
of armament, is a weapon system. OSO has been, for over a year, assigning hull coordi- 
nators (for such types as DDG, DLG, and CG) to see each ship through initial outfitting. The 
Aviation Supply Office, similarly, has provided coordination of initial support of new aircraft 
types. 

Many of the elements of weapon system support such as integrating individual pro- 
visionings of equipments which comprise a weapon system, COSAL development, material 
distribution and stock positioning, and filling initial allowances require a certain amount of 
coordination as a weapon system in order to be handled effectively. Usage data collection 
and analysis for the purpose of maintaining APL's for COSALS cannot be effective except on 
a weapon system basis. 

It is submitted that the program support Supply Demand Control Points are best 
suited to this type of work and can accomplish it with the least amount of modification of 
organizations and make ready time. What would be required is a positive communication 
channel to and from the fleets to bring the inventory managers closer to the operators ona 
weapon-system basis. An expansion of the hull coordinator and aircraft support liaison con- 
cepts of OSO and ASO, respectively, appears to be in order. It does not necessarily follow 
that it is necessary to reorganize personnel and add new organizational units for this pur- 
pose. It would seem more desirable to program as many as possible of the functions of 
sustained weapon system support on electronic data processing equipment. For example, 
mathematical decision rules are now being applied to weapon system provisioning. A by- 
product of this exercise is the equipment APL for the fleet COSALS. It follows that usage 
data and failure data (actual) could be fed into the computer to redetermine the previously 
estimated onboard range and depth of repair parts. 

A program of sustaining the support of equipments on a weapon system basis in the 
program SDCP's would effectively ''close the loop" to the operating forces for weapon sys- 
tems. It would provide the means for feeding back the data to evaluate and improve initial 


provisioning decisions. It would be the means of insuring in the supply availability program 
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that ships would receive true onboard allowance overhauls as well as equipment overhauls. 
A demonstrated capability inthe SDCP's to perform these functions would do more than any- 
thing else to preclude proposals for the establishment of special logistic offices to support 


new programs and weapon systems. 


SUMMARY OF CONCLUSIONS AND RECOMMENDATIONS 

Actual experience in several branches of the armed forces has demonstrated the need 
for concentrating the program support of each defined weapon system in a single organiza- 
tion. It is recommended that the Navy Supply System continue to focus upon the needs of the 
operating forces and the specific requirements of weapons system in making weapon system 
program support assignments to inventory managers. 

Clear, concise statements of weapon system operation and maintenance, when stated 
in a maintenance policy are the best guide lines for program support SDCP in picking up 
logistic support. It is concluded that whenever one is already not provided, the Inventory 
Manager should press the responsible technical bureau for maintenance objectives stated as 
a policy by the time of the first system equipment provisioning. 

There are a number of justifications to the fleets, technical bureaus, SDCP's and 
equipment contractors for the early adoption of computer processing of mathematical pro- 
visioning rules to raise, uniformly, the quality of provisioning decisions. It is concluded 
that BUSANDA should continue to press for the acceptance of mathematical provisioning 
decision rules by technical bureaus and program supply demand control points. Pockets of 
resistance to these rules are present. The resistance is usually born of unwillingness to 
try to understand the rules. An education program is necessary—but, as time goes and as 
the program gains favor, the remaining resistance should be treated with calculated 
impatience. 

Directly related to the provisioning decision rule conclusion is another concerning the 
assignment of parameter values for making decision rules. Not only are parameter value 
assignments necessary for completing the mathematical calculation, they are invaluable as a 
technique for demanding that the equipment designer, builder, and provisioner alike sharpen 
their thinking on repair part decisions and thereby do a much better job. It is recommended 


that the technical SDC P's and bureaus who have provisioning responsibility give early and 


sustained emphasis for defining parameters and setting parameter values for provisioning. 


The value of failure and usage reporting is well recognized. On the other hand the 
inventory managers and others responsible for maintaining allowance lists have very little 
hope of maintaining, through usage data all APL's. 

It is recommended that the mission effect element of military essentiality of equip- 
ments be used as the criterion for preference ranking of equipments for the scheduled main- 
tenance of APL’s on these equipments. 

The establishment of the Polaris Material Office is justified on the basis of the un- 
usual urgency and importance of the Fleet Ballistic Missile Program. It does not 


‘ 
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necessarily set a precedent for all subsequent weapon systems. It does fill what could be a 
void in subsequent weapon systems unless an alternative is adopted. 

It is recommended that the program support supply demand control points become 
the material offices for new weapon systems by providing the means for "closing the com- 
munications loop" with the Fleet on a weapon system basis in order to insure that adequate 


support of weapon systems is provided through the life of these systems. 
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A NOTE ON PYRAMID BUILDING* t 
Zivia S. Wurtele 


It is of interest for planning purposes to determine the effect on the rate of balanced 
growth of adding to an economy industries whose outputs are used neither directly nor indi- 
rectly for the production of consumption goods. Examples of such industries may be found 
among defense industries; and a classic example is Keynes' pyramid building. 

Let us assume that we have a closed Leontief-type economy which is divided into n 
sectors and which is lagged so that the outputs of one period are the inputs of the next. LetA 
and B, ae Veale be n-order matrices of — coefficients a and capital coefficients 
bi . Let x bea column vector of the outputs x; (i=1, ... , n) of the sectors during the period 
t. We assume that (I-A) is a Leontief-type decomposable matrix and that the industries of the 
economy have been aggregated into sectors so as to achieve a matrix of technical coefficients 


of the form: 


where A. (e=l, ..., k) are square, indecomposable matrices of order ny with 


e=1 
where all entries above and to the right of A, are zero, and where the last sector of A is labor. 
Thus the sectors in the ” block Ay for e <k, are needed neither directly nor indirectly for 
consumption. Similarly, let B. be the submatrix of capital coefficients corresponding to the 
coefficients of Aw C= 4s. seg ee 
The equations of our system, lagged so that outputs of one period are inputs of the 


next, are: 
(1) x = as”* + B(x‘*! - x') 


Setting x! = (1+ a)’ x? in Eq. (1), we obtain: 


“Manuscript received July 1, 1960. 
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(2) ( I/a - (I-A)? (A+B)) J” 0%. 





By means of some fairly straightforward algebraic manipulation, it can be shown that Eq. (2) 


is of the form: 





(va-c ) 0 0 y 
_" I/d-C peor 0 ; 
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where the matrices C, = (I-A,)" (Ag + Be) are indecomposable and ve (e=1, ..., k) isa 
vector of the outputs x° i sp so of the sectors of the et? 
Ny+Not --- not! nj+ng+ ... Ne 


block. 





Let 1/d, be the maximal characteristic root of C, ; and assume that if m ¢ e, = + d.: 
We now make the additional assumption that the matrix of coefficients in Eq. (3) is not separa- 
ble. (The separable case may be studied by first dividing the matrix into nonseparable sub- 
matrices.) The elements of C, are positive. If Eq.(3) holds for nonnegative y’, then d must 


1 
equal qd. for some e such that 1<e<k. Suppose d= ¢.. Since 1/4; - C, + 0, we must 
have ret = 0. Since ret = 0 and |i, - Cy | + 0, we must have Yo = 0. In this way, we may 
prove that vy, = Oforh<j. Since c, is indecomposable, ; >0. Consider now the eth set of 


equations of Eq. (3) where e >j. These can be written in the form: 
(4) y? =(1/a,-c.)i2z 
e j e e’ 


where Ze is a nonnegative vector. A necessary and sufficient condition for the solution ve to 
be positive is that qd, ae . ™ ” 

Suppose that to an economy consisting of the sectors of the} through k™ groups, with 
balanced growth rate 4, » are added the sectors of the (j-1)** group. There are two cases to 
distinguish. If as <4, , uniform growth is possible with positive outputs for all sectors of 
the enlarged economy, and the new growth rate is qa - On the other hand, if dj_4 me, » uni- 
form expansion with positive outputs for all the sectors in the 9-1" through kth groups 
would be physically impossible. In evaluating a plan for enlarging the economy by the addi- 
tion of industries whose outputs are not used directly or indirectly for consumption, it is thus 
appropriate to determine whether it is compatible with balanced growth and positive outputs 


for all the sectors of the economy. 


1 The argument of this paragraph rests heavily upon the results in [1] and[ 4]. 
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AIRCRAFT COMPARTMENT DESIGN CRITERIA 
FOR THE ARMY DEPLOYMENT MISSION*T 
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ment of army equipment to areas of actual or potential warfare. The 


large number of vehicles required to support army forces makes 
Weight and cube, the traditional measures of airlift capacity, insuf- 


ficient for analysis of the amount of such cargo which can be carried 
by various aircraft. This paper uses a representative list of army 


An important mission of military cargo aircraft is the deploy- 
| 
| 

deployment cargo to investigate the capacity of alternative aircraft | 

designs. A loading model is developed to estimate the relationship } 
| between cargo compartment dimensions and loading efficiency. The 
relationship between compartment floor area and floor area of vehicles 
loaded is found to be remarkably stable over a large range of alterna- 
tive compartment sizes. A second model is used to minimize the 


distribution of total weight per aircraft when the aircraft are loaded 
in accordance with the criteria of the first model. The resultant dis- 
tributions, presented for various compartment sizes, suggest the 


relationship between cargo compartment and payload capacity which 
would be optimum for this particular mission. 
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INTRODUCTION 

Aircraft loading problems are intimately connected with design criteria. Transport 
design affects both the preferred loading technique and the maximum amount of materiel the 
uircraft can transport over any given route. One of the most persistent problems in air 
transportation has been how to determine the capability of cargo aircraft to perform various 
tasks. It has long been recognized that it often may be either impossible or impractical to 
load aircraft with the max.mum payload allowable over the chosen route; when this is so, some 
measure other than weight must be used in calculating airlift capability. Cargo cubage ap- 
pears to be an acceptable measure for the type of materiel airlifted by commercial carriers 
and for most equipment moved by the m litary :n peacetime; in both cases, most of the cargo 
consists of small units which lend themselves to flex.ble arrangement, so that the cargo com- 
partment space can be efficcently used in all three dimensions. Under these circumstances, 
cubage may be used as the sole measure of space requirements in airlift planning. 

The important relationship for peacetime cargo 1s between its weight and its cube. If 
an aircraft's allowable payload over a route exceeds the weight of the cargo which will fit into 


the cargo compartment, cubage is the constraint on airlift capability. 








Manus« ript rec eived Cctuber 3, 1960. 


This study was performed while the author was a member of the Logistics Department of the 
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If the allowable payload is less than the plane could carry, weight is the con- 
straint. Considerable research has been devoted to this important problem area and the 
implications for airlift planning and aircraft design have been examined in some detail.! 

Unfortunately, the procedures for studying peacetime cargo airlift do not appear ap- 
plicable to army deployment, which involves a different kind of cargo. The bulk of the cargo 
consists of vehicles which cannot easily be arranged to fill the cargo compartment. The 
cargo cannot be stacked to take full advantage of compartment height; most of them are longer 
than the compartment width, so they cannot be turned sideways; and it is impractical to turn 
them on their sides. Cubage alone, then, is not an adequate measure in this case. Detailed 
investigation of the space requirements of army deployment cargo is necessary in order to 
obtain a more accurate measure. Such a measure can then be compared with cargo weight 
and fruitful design criteria for this mission derived. 

This article describes a method which can be used to attain these goals. To aid the expo- 
sition a particular application of the method will be described. We will examine the problems of 
loading a particular list of equipment: that required by a hypothetical planning force designed by 
the Office of the Deputy Chief of Staff for Logistics of the U.S. Army.* The conclusions reported 
here will apply to deployments similar to those of the force studied, of course. Moreover, although 
it cannot be said for certain that other deployments will have similar requirements, it is extremely 
likely that they will, since most of the vehicles studied are standard equipment for almost every 
sizeable army organization and make up a relatively constant percentage of any major force. The 


analysis should thus apply to a number of alternate mixes. 
This article considers the following problem. There is a specified list of army deploy- 


ment cargo to be airlifted. All the items can be placed in aircraft with cargo compartments 
at least 110 inches high, 120 inches wide, and 400 inches long. A number of aircraft designs 
are being considered for this task. Each has a cargo compartment at least 110 inches high but 
widths vary from 120 to 150 inches while compartment lengths vary from 400 to 1200 inches. 
Which design is best for this mission? What is the desirable relationship between compart- 
ment dimensions and payload capacity? 

We will approach the answer to these questions in two stages. First, a loading model 
will be developed to estimate the relationships between cargo compartment dimensions and 
loading efficiency. Compartment width will be examined first, followed by compartment 
length and cubage. A second model will then be developed to minimize the distribution of total 
weight per aircraft when the aircraft are loaded in accordance with the criteria of the first 
model. The resultant distributions, presented for various compartment sizes, suggest the 
relationship between cargo compartment and payload capacity which would be optimum for the 


specified mission. 


_ for example, R. E. Bickner, Cargo Density and Air Transportation, The RAND Corpo- 
ration, Research Memorandum RM-1380 (ASTIA Document No. AD 85398), May, 1955: and 
idem, Cargo Density and Airlift, The RAND Corporation, Research Memorandum RM-1853 
(ASTIA Document No. AD 112417), January 14, 1957. 








*A copy of this list of equipment was made available to the RAND Corporation in September 1959. 







































AIRCRAFT COMPARTMENT DESIGN CRITERIA 
COMPARTMENT WIDTH 


Compartment width claims our first attention; this section attempts to specify the most 





efficient method of using it in loading vehicles. The problems of compartment length are dis- 
regarded here, to be taken up in the following section. 


For expositional simplicity, assume it 1s possible to build an aircraft of any length but 


PEO with a specified width; given this constraint, we wish to determine the shortest possible air- 
craft which will just hold all the vehicles to be transported. The compartment is wide enough 

or to hold some vehicles side by side, but no more than two, since the narrowest vehicle is more 

” than one third the width of the compartment (this applies to all compartment widths studied, 

led since the narrowest vehicle is over 50 inches wide). 

’ The technique is simple analytically and could be applied on the flight line. The vehi- 

t cles are first lined up with their left sides against the left wall of a very long cargo compart- 
ment; the widest vehicles go in first. Graphically, this is merely the cumulative column 

aia co length as a function of vehicle width. In order to account for the space required for hold-down 

ems of 


cables or Straps, as well as problems of maneuvering the vehicles into aircraft, the vehicle dimen- 


gned by sions have been increased. Both the vehicle widths and lengths reported here include this adjust- 


eported ment. Figure 1 shows the function for the vehicles in this study. The horizontal axis gives the 


» although column length which can be made up of vehicles whose widths exceed the corresponding value on 


extremely § the vertical axis. Column lengths are expressed in terms of the total column length (the length of 

















edad. the vehicle column when no doubling of vehicles is permitted). 
ce. The Next, the last—and narrowest—vehicle pulls out of the line and drives forward, with its 
right side against the right compartment wall. The next-to-last vehicle follows it. The 
eploy- 
<i process continues until one of the vehicles can drive no farther because the space is too narrow 
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between the vehicle to its left and the compartment wall to its right. It is easy to perform this 
process graphically by using a transparent overlay of the original distribution: the overlay is 
rotated 180 , with the original horizontal axis placed on the horizontal line corresponding to 
the compartment width; the overlay is then moved from right to left until the two curves touch. 
Figure 2 shows the result of applying the techmque to a compartment width of 150 inches. The 
total length of the new column is about 0.72 of that of the original. The minimum column length 
has been attained because the vehicle represented by the shaded area can move no farther for- 


ward (to the left in Figure 2). 
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Figure 2 - Aircraft loading: compartment width 150 inches 


An intuitive argument for the use of this technique is relatively simple. The vehicle 
represented by the shaded area cannot move to the left of the line EF since the original vehi 
cles (Groups I, II, and IV) do not leave enough space. But the vehicles in Groups VI and VII 
are at least as wide as the vehicle in the shaded area since we have arranged them in order | 
decreasing width. Thus, they cannot be placed in the area to the left of EF. Therefore, the 
only vehicles which can possibly fit in this area are those in Groups III and V, which have 
already been placed there by the technique. Since no further rearrangement of vehicles can 


thus decrease the length of the column, the minimum column length has been attained. 
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By repeating this procedure for each of the compartment widths to be evaluated, we 
1S can find the corresponding minimum possible column length. This can be done in a matter of 
oO minutes once the transparent overlay of the initial distribution is made. In Figure 3 the solid 


ich. line indicates the minimum column length associated with each of the aircraft compartment 
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The horizontal axis in Figure 3 represents compartment width, while the vertical axis 
represents the total length of cargo compartment(s) required for the entire group of vehicles. 
Any rectangular hyperbola of the form xy = K represents all combinations of compartment 
width and column length which have the same total floor area. Three such functions are shown 
by the dotted curves in Figure 3. (Note that since the axes do not start at zero the hyperbolae 
appear almost linear.) We know the total floor area of the vehicles themselves, so any total 
compartment floor area can be related to the vehicle floor area by means of a load factor. 
Thus, a load factor of 76 percent means that the vehicles will cover 76 percent of the total 
floor area; stated in another way, the total floor area equals 1/0. 76 times the required vehicle 
floor area. 
For the vehicles included in this study, Figure 3 shows that the simple rectangular 
hyperbola provides an excellent equation for estimating the relationship between aircraft width 
and column length. The function everywhere lies between the curve representing a 76-percent 
load factor and that representing an 83-percent factor, and is well approximated by that 
representing a 79-percent factor. Within the range considered, we may conclude that air- 
craft width is irrelevant: 79 percent of the compartment floor area is covered with vehicles 
regardless of compartment width. 

There is no a priori reason that a constant floor-area load-factor should adequately 
summarize the effect of varying compartment width. Its use will simplify the remainder of 
the analysis, but the analysis could also be performed with some less convenient estimating- 
equation for this relationship. Other lists of equipment might yield a curve which could be 
better approximated in some other manner. But the list under study here was selected be- 
cause it was thought to be typical of army cargo to be deployed in contingency missions, and 
it is likely that many alternative lists would exhibit the same relationship between aircraft 
width and column length. 


COMPARTMENT LENGTH 

In the preceding section we minimized the length of the total column of vehicles to be 
airlifted, subject to various constraints on the compartment width. We concluded that if 
compartment length could be perfectly utilized, 79 percent of the total floor area could be 
covered with vehicles. This is an impossible ideal; to assess the real possibilities we must 
now examine the effect of compartment length on aircraft loading. 

The same physical fact prevents perfect efficiency in using either compartment length 
or width: vehicles cannot be chopped into convenient sizes. But compartment length is more 
likely to be used efficiently than is width, because compartment lengths are greater in com- 
parison with the average vehicle length than compartment widths are in comparison with the 
average vehicle width. Moreover, we would expect longer compartments by nature to be more 
efficient than shorter ones. Obviously, the approximate magnitudes of these relationships 
need to be determined. 

There are no simple rules for the optimum use of compartment length; and in assessing 
the capability of various aircraft, no rule should be considered at all if it cannot be applied 


with some speed on the flight line. A further complication arises when we attempt to reconcile 
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a technique for obtaining optimum use of compartment length with the technique described in 
the preceding section for optimum use of compartment width. We shall initially disregard 
width (as we did length, in the preceding section) and return to it later. 

Since there is no simple optimum rule for length, we will have to look for some simple 
and fairly efficient working rule. Tentatively, we will use the following technique. 

We are given a group of vehicles to be loaded in single file, none side by side. The 
aircraft have a given compartment length. The goal is to load the vehicles but use as few air- 
craft as possible. 

The vehicles are first lined up in order of decreasing length. The longest vehicle goes 
on the first aircraft, then the second longest, and so on until the remaining space is too short 
for the next vehicle in line. At this point, the longest vehicle which will fit is loaded. 

The general rule is to load no more than one vehicle of any particular length on any 
single aircraft. The rule must be broken toward the end of the process, however, since one 
each of the few remaining vehicle lengths cannot fill the aircraft. In such a case it is permis- 
sible to load more than one vehicle of the shortest available length in a given aircraft. 

When the vehicles are all loaded, we count the aircraft used; this quantity, times the 
compartment floor-length of the individual aircraft, gives the total compartment floor-length 
required. The total length of the vehicle column, divided by this total compartment floor- 
length, gives a loading-efficiency ratio with respect to compartment length. 

The rule suggested has virtues other than simplicity. It imposes at least a small 
degree of combat loading; that is, it makes the loss of any one aircraft somewhat less crucial 
since the great majority of the aircraft contain only one vehicle of any given model. Further, 
it will tend to reduce the a.stribution of weight per aircraft. $ Finally, it may prove fairly 
efficient. 

To get a first approximation of the efficiency of this rule, it was tested on the full list 
of vehicles, assuming that none would be placed side by side. The longer compartment was 
more efficient, as we expected. Cargo compartments of 500 and 1000 inches were evaluated; 
the efficiency of the shorter compartment was 95.5 percent, of the longer, 98.1 percent. The 
size of these percentages indicates that the rule is surprisingly efficient and quite satisfactory 
for our purposes. 

When we turn to a consideration of double columns of vehicles, we once more encounter 
the problem of compartment width, whose optimization was discussed in the preceding section. 
Doubling also compounds the length problem, however, since the end of a vehicle .n one 


column may come in the middle of a vehicle in the other column. If the resulting ioss of 





._— are two constraints on aircraft capacity: the vehicle capacity (determined by the above 
analysis) and the payload capacity. The payload per aircraft loaded according to vehicle 
capacity is the ratio of these two dimensions for the equipment to be moved. The more this 
distribution clusters around its mean, the less inefficiency will arise because of the presence 

of the two constraints; and hence the smaller the number of aircraft of any given design which 

will be required to move the equipment. 








































































388 W. F. SHARPE 
length-efficiency is substantial, it may even be wise to consider less efficient methods of using 
compartment width; but as a first cut at this problem, we will see how efficiently length can 

be used if we retain the doubling techniques described in the preceding section. 

The method adopted can best be seen with the aid of Figure 2, which shows the arrange- 
ment of vehicles in an aircraft 150 inches wide. In any aircraft, there will be some wider 
vehicles which will have no others alongside them (Group I, in this case). At the end of the 
column, there will also be a few pairs of vehicles with the same or nearly the same width 
(Groups VI and VII); we can regard these pairs as single vehicles and mix them at will with the 
wider vehicles (Group 1). We can now apply our length-utilization rule to this composite group 
of vehicles to determine the number of aircraft required to hold them. 

Next, we take the remaining vehicles and divide them into groups according to the 
following criterion: the vehicles can be rearranged within a group in any manner, without 
violating the width constraint. Groups II, Il, IV, and V in Figure 2 are examples of such 
groups. Next, we apply the length-utilization rule to each such group, obtaining the number 
of aircraft required. The aircraft required for two groups which go alongside each other is 
merely the larger of the aircraft required for the individual groups. ‘Thus, if Group II requires 
more aircraft than does Group III, the aircraft required for Group II will suffice to carry both 
groups. This technique allows us to determine a number of aircraft which can contain all the 
vehicles. We may express the efficiency of the technique by dividing the length of the initial 
(doubled) column by the total compartment length of the aircraft required. The results for 


the four extreme compartment sizes are as follows: 














Compartment Compartment Width 
> 120 in. | 150 in. 
(®) (%) 
500 91.8 89.4 
1000 95.4 91.2 














Again the efficiency is surprisingly high. The narrower compartment uses length more effi- 
ciently since relatively few vehicles have been doubled, while the wider compartment pays for 
its efficient use of width by achieving a slightly less efficient use of length. The advantage of 
longer compartments is evident for each compartment width. 

These cifferences are not insignificant; they suggest that a really sophisticated analysis 
might profit from a study which applied several possible techniques to this problem, evaluating 
the efficiency of such techniques for a large number of possible compartment sizes. The re- 
mainder of this study, however, adopts the very simple assumption that length can be used 
with 91.5-percent efficiency regardless of compartment size. This assumption has the virtue 
of simplicity and should not prove far wrong for most compartment dimensions under con- 


sideration, since compartment length and width are usually correlated. The efficiency of the 
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short, narrow compartment is 91.8 percent, while that of the long, wide compartment is 91.2 
percent. If a compartment is to be considered which sharply diverges from these proportions, 
the final values obtained in the analysis should be adjusted appropriately. 

We may now summarize the analysis of this section and the preceding section. In the 
preceding section, we found that vehicles would cover 79 percent of the floor area in a perfect 
utilization of compartment length; if the policies suggested in this section are adopted, how- 
ever, only 91.5-percent efficiency appears obtainable in terms of length utilization. In that 
event, only about 72 percent of the floor area will actually be covered—91.5 percent of 79 
percent. Stated another way, the minimum number of aircraft required to hold the vehicles 
can be computed by dividing the vehicle floor area by 0.72 to determine the aircraft floor area 
required. This figure is then divided by the floor area of the compartment size in question, 
to determine the required number of aircraft. And this technique applies to the full range of 


aircraft under consideration. 


COMPARTMENT CUBE 

Our analysis thus far has been concerned only with the vehicles to be airlifted. The 
study must now encompass the remaining materiel, which accounts for one third of the total 
weight to be moved. This section deals with the requirements for cubage in which to place 
this equipment; the next section deals with its weight. 

Having determined the minimum number of aircraft required to move the vehicles, we 
must now find out if enough cubage is available in those aircraft to carry the remaining equip- 
ment. Two sources of space are available for this purpose. First, there is the space not 
covered by vehicles. This amounts to 28 percent (100 - 72) of the total aircraft floor area, or 
approximately 40 percent of the vehicle floor area (28/72). The second source is the platform 
area of many of the cargo vehicles. This area is roughly 30 percent of the total floor area of 
the vehicles. The total area available for the remaining cargo is thus equal to approximately 
70 percent of the vehicle floor area. 

The cubage of such equipment in the lists used in this study was approximately 1.4 
times the floor area of the vehicles. (It is considerably denser.) If it were spread evenly 
over the available area the average stack height would thus be 2 feet (1.4/0.7). Since the 
minimum aircraft height considered is over 9 feet from the floor, and between 5 and 6 feet 
from the average truck bed, there is considerable latitude in loading this equipment. For all 


practical purposes, cubage per se can be disregarded in the analysis which follows. 


AIRCRAFT PAYLOAD 

We have determined the minimum number of aircraft required to hold the vehicles to 
be moved, and we have found that this number can easily hold the remaining cargo. The next 
step is to examine the distribution of weight per aircraft which would result if aircraft were 
loaded in the way the preceding Sections have described. By comparing this distribution with 
current designs, we can determine which constraint is more important for such aircraft— 
weight or floor area; more important, we can suggest desirable design criteria for future 
aircraft. 
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We will begin with the distribution of weights which results from loading the vehicles, 
temporarily leaving the remaining items out of the analysis. First, we must determine each 
vehicle's weight per square foot of floor area. For convenience, we will use the term Area 
Load Factor (ALF) to indicate the ratio of weight to floor area. We then calculate the distribu- 
tion of these ALF's, using floor area as the frequency measure.” In the particular case de- 
scribed in this study the mean ALF was 47.6 lb/ rt? and the standard deviation of a fitted 
normal curve was 33.4 lb/ft. 


The average floor area of the vehicles studied was 108 rt, and we have determined 


uo pe se o> 


that 72 percent of the floor area of any given aircraft can be covered with vehicles. There- 
fore, the average number of vehicles per aircraft (N) is a simple function of the aircraft floor 


} 
area in square feet (A): 


0.72 A. 
108 
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If a large total number of vehicles is to be loaded, and randomly with regard to weight 
per square foot of floor area, it is a simple matter to estimate the distribution of the average 
ALF's of the individual aircraft. Each aircraft represents a sample of size (N) drawn from 
the population of mean 47.6 lb/ it? and standard deviation 33.4 lb/ it. The average ALF (per 
square foot of vehicle floor area) in a given aircraft corresponds to the sample mean, and the 
distribution of such average ALF's is merely the distribution of the sample means. This 
distribution has simple properties. The mean is, of course, the mean of the population, or 
47.6 lb/ rt. The standard deviation is simply the standard deviation of the population divided 
by VN (i-e., 33.4/VN).” 


cr 


This is the proper measure for the special case where floor-area utilization is constant. In be 
the general case the particular measure which best represents the amount of an aircraft ite 
utilized by each vehicle should be chosen, 


4 ; ae ; , it j 
In the general case N would have to be determined for each individual compartm«t:.t size in 
accordance with the results of the analysis of length and width utilization. mc 


Of course this will give a somewhat optimistic estimate of the actual standard deviation, be- 


cause N is merely the average number of vehicles loaded in each aircraft; since the actual * 
number varies around N, the standard deviation of the resulting distribution will exceed the cre 
value calculated by this method, In fact, however, the error resulting from the assumption ; 
that N vehicles are loaded in each aircraft appears to be small. The actual weights loaded in wil 
the smallest aircraft under consideration--120 inches wide by 500 inches long—were computed tio 
and compared with the distribution predicted by the use of this formula. Approximately 90 

percent of the aircraft contained loads within the range which the formula predicted would con- cra 
tain 95 percent of the cases; approximately 95 percent of the aircraft loads were within the 

range which the formula predicted would include 99 percent of the cases, While the error in the 
this case was not insignificant, it appears that using the formula will not seriously bias the ; 

study, for two reasons, First, no attempt was made in the example to adjust loads so as to an 


reduce the range of vehicle weights; such a readjustment would have made it possible to ob- 
tain a distribution much closer to that predicted by the formula. A second reason for re- 

garding this error as acceptable is that the particular aircraft chosen for the test is the rep 
smallest considered, having a floor area of only 420 ft®; the error for larger aircraft is 
likely to be considerably less important. 
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Thus, we can easily determine the distribution of average weight per square foot of 
vehicle floor area once we know the floor area of the individual aircraft. The next step is to 
consider the effect of loading the remaining cargo. 

Figure 4 shows the cumulative distribution of weight per square foot of vehicle floor 
area for an aircraft of 450 tt? floor area. Area is expressed in terms of total vehicle floor 
area, the ALF is measured in lb/ tt? of vehicle floor area; the cumulative distribution is shown 
by the curve OPNR. Now we must add to this weight another group of items which weigh half 
as much as the vehicles themselves. Since the preceding section demonstrated that there can 
be considerable latitude in loading these items, they will be added here so as to yield a desir- 


able weight distribution. 
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Figure 4 - Cumulative distribution of area load factors for 
aircraft of 450 ft2 floor area, loaded with vehicles 


What weight distribution would be the best? If the aircraft is so designed that all air- 
craft loaded with vehicles already contain more weight than the allowable payload, weight will 
be the only effective constraint on aircraft loading; the manner in which we load the remaining 
items is irrelevant, since they will eventually have to be reloaded in other aircraft. However, 
itis much more likely that only a portion of the aircraft loaded with vehicles already contain 
more than the allowable payload weight. In these aircraft, weight will be the effective con- 
straint on cargo carried, but vehicle loading will be the initial constraint in the remaining air- 
craft. Now if additional weight is added to the former aircraft, the number of aircraft required 
will likewise increase—weight being the constraining dimension. However, if we load the addi- 
tional weight on the latter group of aircraft (those with a low initial ALF) the number of air- 
craft required will not be increased—up to a point, at least—since the initial loads are below 
the allowable payload. Thus the preferred loading policy for this remaining equipment is quite 
simple: each pound should be loaded in the aircraft with the (then) smallest average ALF. 

This rule can be easily performed graphically. Since the horizontal axis in Figure 4 
represents ALF's measured in Ib/tt?, and the vertical axis represents floor area measured in 
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rt”, an area of any given size on this graph represents some total weight. Since our task is to 
distribute a given total weight, we need only apply the specified rule until the required area is 
obtained graphically. The rule that each pound must be loaded in the aircraft with the smallest 
average ALF can be followed by merely moving a vertical line from left to right. The area to 
the left of such a line and under the curve will represent the weight added; the line is merely 
moved to the right until this area equals the total weight to be distributed. In Figure 4 the 
final position is shown by the line MN, while the shaded area OPNMO is equal to the weight of 
the additional cargo. The final distribution of ALF's is now given by the curve OMNR. 

Thus far, we have determined the distribution of weight per square foot of vehicle floor 
area. We must now convert this into the weight per square foot of aircraft floor area. This 
is a simple operation since the vehicle floor area equals 72 percent of the aircraft floor area. 
By merely multiplying the ALF's already obtained by 0.72, we obtain the corresponding ALF's 
stated in terms of aircraft floor area. 
Figure 5 represents the final distribution of weights per square foot of aircraft floor 
area for a cargo compartment of 450 rt? floor area. Consider an aircraft which can carry a 
payload over the relevant range equal to 52 ib/tt? of floor area. From Figure 5 it is evident 
that 90 percent of the aircraft loaded up to the floor-area constraint will not require the full 
payload available; for these aircraft, the limiting factor is their floor area. To determine the 
number of aircraft required to carry this equipment, the total vehicle floor area is divided by 
0.72 and then multiplied by 0.9; the resulting figure is then divided by the floor area of the 
individual aircraft. Now the remaining 10 percent of the aircraft (as initially loaded) contains 
cargo which weighs more than the allowable payload; thus some of the equipment must be off- 
loaded and placed on additional aircry ie Weight will be the limiting factor in loading this 
equipment. The total weight is shown by the shaded area in Figure 5; the total number of air- 
craft required to carry these items can be determined by dividing the total weight by the 
relevant payload of the aircraft. This procedure, although somewhat cumbersome, can pro- 
vide a rather accurate estimate of the number of aircraft of a particular design required to 
move such a group of equipment. 

Let us now summarize the ALF relationships for aircraft of various floor areas. In 
Figure 6, the minimum ALF in 1b/tt? of aircraft floor area is shown by the curve labeled 
"0 percent," for aircraft floor areas from 400 to 1100 tt”. It has a very slight slope and is 
asymptotic to 51.4 lb/ft”. The other three curves show the values below which lie the average 
ALF's of 95.0, 97.5, and 99.5 percent, respectively, of the aircraft for each compartment 
size. The equations for these functions can be readily determined from the confidence limits 
of the normal distribution and the relationship among the aircraft floor area, the number of 
vehicles per aircraft, and the standard deviation of the relevant ALF distribution described 
earlier in this section. 


But it cannot be placed in any of the aircraft which are floor-area-limited since this equip- 
ment also requires floor area. 
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Figure 5 - Cumulative distribution of area load factors for aircraft of 
450 ft2 floor area, loaded with vehicles and other ca rgo 
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Figure 6 - Confidence limits of area load factors for 
various aircraft compartment sizes 


One of the aims of this paper is to arrive at relationships which are simple but none- 
theless take into account the more important factors to be considered in evaluating the cap- 
ability of various aircraft in performing the army deployment mission. In keeping with the 
quest for simplicity, it seems reasonable to suggest that, at least for the larger aircraft, no 
great harm would be done if the distribution in Figure 6 were represented by a single line at, 
say, 51 tb/tt?, ' If this simplification is made, a straightforward rule can be adopted. Any 
aircraft which is designed for a maximum payload (over the relevant range) of less than 51 


Ib/tt? of its compartment floor area will be weight-limited in performing this mission; and the 


7 
The slight loss of accuracy from such an approximation should be well worth the resultant 
gain in simplicity. 
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total number of aircraft required can be determined by dividing the total weight to be carried 
by the aircraft payload. On the other hand, any aircraft which has a payload greater than 51 
lb/ft? of floor area will be floor-area-limited in performing this mission, and the total number 
of aircraft required can be determined by first dividing the total floor area of the vehicles to be 
moved by 0.72, and then dividing the result by the aircraft floor area. 


CONCLUSIONS 

As it is formulated here, the design criterion which indicates an aircraft's efficiency 
in performing the army deployment mission is fairly obvious. The most efficient design would 
provide for a payload capacity (over the relevant range) of approximately 51 lb/ft” of floor 
area. If the payload were greater than that, the payload limitation could be reduced (thus 
reducing the total cost of purchase and operation of the aircraft) without decreasing the amount 


of cargo which the aircraft could in fact carry. On the other hand, if the payload were less 


than 51 lb/ rt’, the compartment floor area could be reduced (again reducing the purchase and 


operating cost) without decreasing the amount of cargo which the aircraft could actually carry. 
Thus, the most efficient aircraft for this mission would have payloads (over the required 
range) equal to about 51 1b/ tt? of compartment floor area. 

' It is important to remember that the army deployment mission is not likely to be the 
sole mission for any aircraft. Thus, the criterion obtained from the application of the method 
described here will be but one of a number of criteria to be considered in determining aircraft 
design. Indeed, there is no a priori reason why aircraft should be designed to accommodate 
existing army equipment designs and mixes. It may well be that the redesign of army force 
equipment, taking into account the constraints on airlift of present equipment, may be the 
better approach. The proper policy will depend on the extent to which changes for the purpose 
of this deployment will harm the efficiency of the vehicles or the aircraft in meeting their 
other requirements, and the importance of the deployment mission relative to those other 
requirements. 
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—— 

Problems of inventory control for items with extremely low 
demand (say less than one unit per month average demand) have re- 
ceived relatively little attention in the literature. However, for mili- 
tary supply systems it is well known that an extren*:ly large propor- 
tion of the total number of items fall into this category. The problem 
of designing a control system for these low demand items is therefore 
one of the most critical inventory management problems faced by the 
National Military Establishment. 


In the following, a single echelon, multidepot supply system is 
studied for low demand items having a stationary Poisson probability 
distribution for demand. Instantaneous information concerning inven- 
tory levels is assumed to be available. Procurement lead time is 
assumed to be constant as well as the time required for either of two 
available modes of redistribution. Items are ordered one at a time 
and decision rules are developed for allocation of new procurement, 
redistribution of stocks among the depots, and for determining system 
and depot stockage objectives in order to minimize the expected costs 
resulting from system and depot stockouts, cost of redistributing 
stocks among the depots, and costs of transportation from the source. 
The model is an extension of previous models in that redistribution 
costs and depot stockout costs are considered in the determination of 
the stockage objective for the system as a whole. When more than a 
single unit is on order, the allocation of a unit ready to be delivered is 
determined by the solution of a dynamic programming problem. When 
only a single unit is on order (and this is the one ready to be delivered), 
the optimal allocation procedure is reduced to allocating the unit to the 
depot which has the greatest probability of using it in a time period 
T + (1/A) where T is the procurement lead time and A the system de- 
mand rate for the item. 











INTRODUCTION 
A supply system consisting of N depots and a central control point will be studied. It 


is assumed that the system operates using transactions reporting so that at any instant of time, 





* Manuscript received August 26, 1960. 


tResearch supported by the Advanced Logistics Research Division, Bureau of Supplies and 
Accounts, U. S. Navy, and United Research, Inc. 


G. HADLEY AND T. M. WHITIN 


the central control point knows the state of each of the N depots. The central control point has 
the authority to allocate procurement, to make redistribution decisions, and to determine how 
much stock will be carried at each of the depots. It will be supposed that the low demand item 
has characteristics which make it desirable to procure it one unit at a time. Demands will be 
assumed to occur one at a time and at random at each of the depots, the demand at one depot 
being independent of that at the other depots. The random variable u representing the number 
of units demanded at any depot j in a time interval t will be assumed to be Poisson distributed. 
The average demand rate a, at depot j is assumed to be invariant with time. Furthermore any 
units sent to a depot will ultimately be demanded. The procurement lead time will be imagined 
to be a constant T. 

For this system a number of problems face the decision maker at the central control 
point. First there is the question of how many units should be stocked at each of the depots. 
For some of the depots it may be best not to stock any at all. Then there is the problem of 
selecting the depot to which each new unit procured will be allocated. As will be seen later, it 
will not always be optimal to allocate the unit to the depot where a demand initiated the order 
for the unit. Finally, there is the problem of deciding when to redistribute and what mode of 
transportation should be used. Going along with this is the similar problem of deciding upon 
the mode of transportation to be used in shipping from the source. A model will be developed 
herein which helps to answer these questions, when two modes of transportation referred to 


as fast and slow are allowed. 


SYSTEM AND DEPOT STOCKAGE OBJECTIVES 

When units are ordered one at a time as demanded, the amount of stock on hand plus 
on order minus backorders for the system remains constant throughout time. Denote by M 
the stock position of the system (on hand plus on order minus obligations). The optimal value 
of M is found by balancing carrying costs against the costs of stockout and redistribution. 
Average procurement costs (for a given mode of transportation from the source) are fixed in 
this model, since units are ordered one at a time, and hence do not enter into the 
optimization. 

It will be supposed that the cost of carrying one unit in inventory for one year will be 
IC where C is the unit cost of the item. Furthermore it will be assumed that there exist two 
kinds of stockout costs, depot stockout costs and system stockout costs. The depot stockout 
cost will be denoted by of at depot j and is the cost per unit per unit time out of stock (cost of 
i unit year of shortage). The system stockout cost per unit year 7 will be assumed to be 
incurred in addition to the depot stockout costs if the entire system has no stock on hand. Let 
p(x) be the probability that the system procurement lead time demand will be x. Take A to be 


the average system rate of demand for the item; then 


N 
= Dd, 
J 


j=! 
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where as indicated previously ., is the average rate of demand at depot j. If the demands at 
the depots are independent and Poisson distributed, then p(x) also has a Poisson distribution, 


j.€., 


x -AT 
AT) e ei 


p(x) = x! 


To introduce the approach, consider first the case where redistribution is free and 
instantaneous. Then there can only be system stockout costs. Furthermore, the question of 
how much should be stocked at each depot is irrelevant. Only the system stock position is of 
interest. It is found by balancing the holding costs and the system stockout costs. Recall that 
the system level M remains constant in time. Also within a leadtime all M units over and 
above the backorders will become available to meet demands. This means that the probability 
that w units are on hand in the system at any instant of time is simply p(M - w),i.e., is the 


probability that exactly M - w units are demanded in a procurement lead time. Thus, the 


expected on-hand inventory at any time is 


M 
>; (M - x) p(x), 
X=O 


and the expected number of stockouts (backorders) is 


[-.¢) 


Z (x - M) p(x). 


x=M+1 


The expected cost per year K(M) is, then, 


oO 


M 
(1) K(M) = IC >: (M - x) p(x) + 7 S (x - M) p(x). 


x=O x=M+ 1 


IfiM minimizes (1), it is necessary that AK(M) = K(M+1) - K(M) 2 0, and AK(M-1)< 0. These 


relations show that the optimal M is the smallest M satisfying 


co 


p | IC 
(2) p(x) <————.- 
x=M¢+ 1 +s 


This result has been obtained previously |1]. The M obtained in this way will be called the 
"free and instantaneous" M since it is optimal when redistribution is free and instantaneous. 
Now consider the case where redistribution is neither free nor instantaneous. In order 
to include the costs of redistribution and depot stockouts in the total cost expression it is 
hecessary to say something about redistribution. The purpose of redistribution is to prevent 


or eliminate stockouts in the period between allocations. For the purpose of setting the system 
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and depot stock positions it will be assumed that a redistribution is never made to a depot 
until there is a backorder. The assumption is reasonable because for the low demand items 
under consideration, the inventory levels at most depots will frequently be either 0 or 1 and 
hence normally a redistribution will not be made until there is a backorder. This does not 
imply that the redistribution rules will never indicate the need for a redistribution until there 
is a backorder. It only implies that the probability of this event is small enough so that it 
does not have an appreciable influence on the system stock position. It will also be assumed 
that when there is a stockout a redistribution will be made and that the stockout lasts for the 
full redistribution lead time, i.e., no allocation arrives during the redistribution lead time. 
With these assumptions, the cost of a stockout at depot j (to be added to any system stockout 


costs) is 
(3) 


where v is the cost of redistribution by the appropriate mode of transportation (the selection 


of this mode will be discussed later) and T,, is the redistribution lead time. Both v and T 


will be average values since they really mane on which depot redistributes to depot j. . 
Since the purpose of redistribution is to protect against stockouts in the time period 
between allocations, it will be necessary to compute the probability that there will be u de- 
mands at depot j in the time period between allocations. Because lead times are fixed the 
time period between allocations will have the same distribution as the time period between 


; -A : 
demands for the system, that is, Ae . when demands are assumed to be Poisson. Then the 


marginal distribution for demand p;(4) at depot j in the time period between allocations is, 


for Poisson demands, 


. u 


A(A.t) 
io =f ang oe +h * 2 +}. gta 


ae «a 


Oo 


which is the geometric distribution. Equation (4) follows immediately on integration by parts. 
Denote by « (é) the probability that there are ~ units on hand at depot j immediately 
after an allocation. The expected number of stockouts at depot j in the time period between 


allocations, if units are on hand at the beginning of the period, is ' 


oo 


(5) 2, 88, (u). 
- u=¢+1 


Note that the probability of more than one demand at depot j per period will normally be small 
so that the expected number of stockouts will be essentially the same as the probability of one 


backorder, i.e., 
oo 


x (u - 2) B (u) 2 Big +1). 


u=é+1 
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Thus, on averaging over bis the expected cost of stockouts and redistribution per year is 


N eo 
(6) aD, D, w- Ob al), 


j=l €— u=f+1 


Since on the average there are A periods per year. 

It appears difficult to compute accurately ai(e ), Since it depends on the allocation 
policy, the redistribution policy and the joint distribution of the stocks at all N depots. We 
shall not attempt to compute the dj(é). Instead we proceed as follows: Let ., be the safety 
stock maintained at depot j. If p is the average procurement lead time demand, then p is the 


average amount on order in the system, and hence 


N 
(7) M=p+) S,. 


j=l 
Now M must be an integer. It will also be assumed that the s; are integers also, since S. will 
normally be on hand at depot j.. This assumption is not completely consistent with Eq. (7) 
since Pp may not be an integer. In order to make M an integer when the s; are integers, p 
will be replaced by the smallest integer greater than p. 
The s; will be 0 or 1 for most depots for this low demand item. To simplify Eq. (6) it 
will be assumed that sufficient accuracy will be obtained by assuming S is on hand at depot j 


after an allocation. Then, the average redistribution cost per year becomes 
N 20 

8 r 1D. (u - s,) Btw. 

(8) j , p,( ) 
j=l u=S}+ 1 


The total yearly expected system costs, exclusive of procurement costs, are then 


y 


N 
(9) TC = K(M) + 7 
j=l )u-Si+1 


u - S.)p.(u). 

( i) b,(u) 

In order for any given set of S to minimize Eq. (9), it is necessary that 
arc(s,)= TC (s, + 1)- TC(s,) > 0, atc(s, ‘ 1)< ee or 


Thus, the smallest Sj » )=1, .-. , nis chosen so that 


ie) 


(10) aTc(s,)= Ic - (IC + 7) 2. p(x) - Am, Zz. p, “2 oO. 
x=M+1 u=Sj+1 
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In general, each of the inequalities (Eq. (10)) depends on all of the S}- However, if the system 
stockout term is negligible, _ the inequalities become uncoupled and the jth inequality depends 
only on S}- In this case, the smallest Sj is chosen so that 


(11) 2. pu) < Jk» i=, 

u=5,+ 1 Tj 
It would seem that normally the system will be decoupled because enough additional safety stocks 
will be needed to avoid depot and redistribution costs so that the system stockout costs will be 
negligible. The difference between the M computed from Eq. (11) and the "free and instan- 
taneous" M will be called the ''augmented inventory" if this quantity is positive. The augmented 


inventory is carried to avoid redistribution and depot stockout costs. 


As the redistribution lead time increases, 7 j increases, and hence the amount of 
augmented inventory required will increase. This observation points out the fact that the above 
system will not be suitable in certain cases. For example, if the redistribution lead time were 
longer than the procurement lead time, it would not make sense to hold safety stocks to avoid 
redistribution. Redistribution would never be used; units would always be ordered from the 
source. Before accepting the above model, then, it should be compared with a completely 


decentralized model in which there is no redistribution. In such a model, a stockage objective 


M; (on hand plus on order minus backorders) is determined by minimizing the costs of carrying 


inventory and stockout costs for each depot individually. Thus the optimal M; is the smallest 


M, such that 
(12) 


where P; (u) is the probability that u units will be demanded at depot j in the procurement 
lead time. In the event that the expected costs for this decentralized model are less than those 


for the centralized system, then the decentralized model should be used. 


If the procurement lead time is long, compared to the redistribution lead time, much 
larger safety stocks will be required for the decentralized model than for the centralized model 
which allows for redistribution. It might be noted that even without introducing redistribution, 
the decentralized model could be improved upon by using an allocation rule which allocated 
units to depots which needed them most at the time the unit is ready to be shipped, rather than 


to the depots which ordered them. 


2 ° 

The reader may question whether or not there should be a system stockout term in addition to 
the depot stockout costs. It was included for generality. It might arise, in practice, as the 
cost of expedited procurement, for example. 
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Selection of Mode of Transportation from the Source 
It will be observed that the mode of transportation used to ship units from the source to 

the depots influences p(u) but not p,(u). Thus, if the system is decoupled, the s; can be 

determined independently of the mode of transportation used from the source. The value of 

the mean lead time demand and hence M will, however, depend on the mode of transportation 

used. The costs of depot stockouts and redistribution are then independent of the mode of 

transportation used from the source; only K(M) depends on this mode of transportation. If 

Cg is the cost of shipping one unit from the source to depot j by mode r the optimal mode 


of transportation from the source is the one which minimizes 


N 


r 
P Cis d; + K(M), 
j= 


where M is the one appropriate to the mode of transportation under consideration. In the 
event that augmented inventory is needed for both modes of transportation, the slow mode * 


from the source will normally be preferred‘. 


It is not hard to see that as demand increases, other things being equal, it does not pay to 
utilize the fast mode of transport from the source. The reason for this is that stockout costs 
and inventory carrying charges change (approximately) proportionately with the square root of 


lead time demand while transportation costs increase proportionately with demand. 


Also other things being equal, it pays to use the rapid mode of transportation for high cost 
rather than low cost items since the potential savings in carrying charges are larger for high 
cost items. Similarly, as 7 increases it pays to utilize the fast mode for sufficiently high 


values of 7. 


Allocation Rules 

The decision as to where to allocate each unit ordered can be made at the time the unit is 
ordered or up until the time the unit is ready for shipment into the system. Clearly, it is best 
to delay allocation as long as possible. Deciding on the allocation at the time a unit is ordered 
may lead to a poor decision. It will here be assumed that the allocation decision is made at 


the time the unit is ready to be shipped from the source. 


ee 


3 
The terminology fast and slow for the modes of transportation does not necessarily refer to 
the mode of transportation used, For example, the same carrier might be used in both cases 


@ \ut the fast mode implies that expediting of the redistribution occurred. In any event the fast 
g™ode costs more and takes less time than the slow mode. 


qi 
@ The possibility of simultaneously using two modes in shipping from the source is not considered, 


4e,, the possibility of shipping to some depots by the fast mode and to others by the slow mode. 
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The allocation decision will be made by minimizing the expected costs. The planning 
horizon over which costs are to be minimized will now be discussed. Clearly, the best allo- 
cation will depend on the time until the next allocation. For example, because of the nature of 
the Poisson probabilities, if the next unit will not be available for a long time in the future, 
then the unit should perhaps be allocated to one of the higher demand depots where it will almost 
certainly be needed. However, if the next unit to be allocated will be available very soon 
it might be better to send the unit to one of the lower demand depots. The best allocation also 
depends on the number of units on order since this will determine the number to become 


available in a procurement lead time. 


Any units on order at the time the decision is to be made will arrive within a time inter- 
val T (the procurement lead time. Furthermore, any units not on order cannot arrive in time T. 
In fact, the time before the next demand is on the average 1/A; if a unit is demanded and 
ordered at this time it will arrive at time T + (1/A). Thus, on the average, there will be no 
additional allocations from the time the last unit on order arrives until time ti 1/A. Alloca- 
tions are made to protect against stockouts until the next allocation. On the average the units 
on order will be the only units available to the system until T+ 1/A. The planning horizon 
will therefore be chosen to be T + 1/A. If it is chosen to be longer than this, units not on 
order must be taken into account. If it is chosen to be shorter, not all available information 
will be used. In reality, of course, all future arrivals have an influence on the present 
decision. However, what happens after T + 1/A is not known and should have little influence 
on the allocation of the present unit. 

A dynamic programming model will be developed to obtain the optimal allocation. This 


model takes into account the times of arrival of all units on order. Suppose that at the time 


the allocation decision is to be made there are w units on order including the one to be 


allocated. These units will be denoted by 1, 2, ... , w, and the times of allocation represented’ 


by tj, tg, --- » ty, where t; = 0 is the time at which the present allocation is to be made. Let 


— . -_ 1 
(13) _ =1, 6 » W-1;TWHT+Z-t,. 


Note that T, 2 1/A. The time between ty, and t,, 1 will be referred to as period n. 


Denote by fy (x) the probability that exactly x units will be demanded at depot j in 


period n. Assume that the inventory on hand at depot j at time ty is L,. 


For simplicity, the time required to ship the unit from the source is neglected and it is taken 
to be on hand immediately after the allocation decision is made. These times are typically one 
third or less of the redistribution times since the redistribution times include the administrative 
lead times as well. 
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At the time the allocation is to be made, there may be one depot with a backorder which 
has no redistribution on the way to meet it. There can be one backorder because if a demand 


which generates a backorder occurs within a time T,, Of the next allocation no redistribution 


M 
ismade. In this case, the unit is automatically allocated to the depot with the backorder. There 
cannot be more than one backorder, since by assumption any others will be met by redistributions 


(even if they occurred within a time T,, of the allocation under consideration). 


M 


Consider now the case where there are no backorders at the time ti In setting up the 
dynamic programming problem it will be assumed that redistributions are made only when there 
isa backorder. It will also be assumed that all backorders are met by redistribution, so that 
at times to, mente ty there are no backorders which do not have redistributions on the way 
tomeet them. The only consequence of this assumption (which considerably simplifies the 
formulation) is that expected stockout costs are slightly overstated since any backorder which 
is met by allocation will avoid the redistribution cost and will also incur stockout costs for a 
time somewhat less than Tu" This error will be very small provided that Tu is not a large 
fraction of the average time between procurements. 

In the dynamic programming problem the costs of stockout and redistribution are mini- 
mized. Carrying costs need not be included since each unit will be carried somewhere in the sys- 


tem. Transportation costs are not included either since on the average they are incurred anyway. 


A to be the number of units allocated to depot j at time tt Since only one unit is 


to be allocated at a given time it must be true that 


Take y 


N 
- - 
(14) y,= Oorland )" yi =1. 
j=1 


Since the demands at the various depots are assumed to be independent, the probability that 


k= (x > as x) is a vector of the demands at the depots in period n is 


“ j 
(15) f(x) = * fi (*]- 
Ifthe stock at depot j is af at time t,, before the allocation of unit n, then the stock at 
depot j is 
4 ee = 
(16) fo 7 & + 7h, POOF E,, * O 


because backorders are eliminated by redistribution. 


Define T. (4 on ate é,) to be the minimum cost for periods n through w when an 


@ °ptimal policy is followed and when the stock at depot j at the beginning of period n (before 


the allocation) is i- Then the recurrence relations have the form 
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_ min 
T (4, eee ty) 


~ 


(17) 2 (x) Tao [fa + yi -x. | . 


; : = i ca >e 


: j .xJ 
0 if $+ Y¥), x S< 0, 


Also 


(18) Ty(Sy> ++ > by) =H 


~ 


The minimum cost is then T(Ly coe 9 Ly) The values yi tell what the optimal allocation 
should be. 


The computational procedure is to compute first Ln fy ar and tabulate it for 


so § 
all the possible fp oe oN’ Then, using ar compute . AE and oe. Se it for all 
possible fi nher td EN? etc., until T, Li Saat og Ly is found along with the optimal vy: If 
the - or u, could vary over a wide range, it would be quite impossible to solve a dynamic 
programming problem of this sort even with the largest computer available today. It is only 
because the assumption is being made that demands are very low that anything can be done 
For almost all the depots, the range of variation of §. will be 0 or 1, although it is possible to 
allow a small number of depots to have Si go upto2or 3. Furthermore, for almost all the 
depots x; =Oorl1. The probability that x; > 1 will usually be negligible. In addition 2 
will normally take on only the values 0, 1, 2. The probability of a system demand of more 
than 2 before the next allocation will be small since on the average, the time between alloca- 


tions is the average time between system demands. 


With the above simplifications, it seems quite possible to solve a two stage problem with 
about 10 depots on a high speed computer in something like 30 seconds or less. It seems 


unlikely that it would be worthwhile in practice to ever try to include more than two stages. 
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This should give all the accuracy needed. Note that in Eq. (17) with the simplifications suggested 


>. - 2(zx, = 0) (zx, = 1) + (Ex, = 2). 


all x 


For example, u(x, = 1) is the sum over all x; such that one and only one x; is unity and all 
other x; vanish. Furthermore in 2x; = 2, all possibilities need not be considered, because the 
probability of two demands is negligible at most depots. There can still exist quite a few 
possibilities for the set of fj since if there were 10 depots and each depot could have qj Oor 1 
there would be 2/0 ~ 1000 possibilities for (£1, -.. , y)- 


In the event that the unit to be allocated is the only one on order, a tremendous simplifica- 
tion can be introduced into the dynamic programming problem since it reduces to a one stage 
problem. The single stage extends over the entire planning horizon T + 1/A. Denote 


vd 4 Wd (u) by Pj *(u). Then the problem reduces to that of minimizing 


N 
A 
1 = : . - L. - y.)p.*fu. 
(19) Zz Z *, pie (¥ L, v5); (,) , 
where 


yj = Oorland Zy;=1. 


* 


Suppose Yip +++ » Yy Constitute an optimal set of y,, that is, a set of y; which 
minimize Eq. (19). Let z* be the optimal value of z. Assume that ve = 1 and all other vj =0. 


If the unit were allocated to depot i # k, then the value of z for this allocation must satisfy 
z-z*>0O. However, 


co 


iin » 2 


KK) -"™ p; (uj) > 0. 


A 
(20) Z-2*=7 i" p 
k uL= L,+1 i uj=Ly+l 


The optimal allocation of the unit is then to the depoi with the largest value of 


x 


* 2a pt Ww. 


u=Ly+1 


In the event that the m are all equal, the result is that the unit is allocated to the depot which 
has the greatest probability of using it in the time period T + 1/A, which depot may or may not 
be the depot whose demand originated the order for the unit to be allocated. 

For many low cost items, it will not be worthwhile to use the dynamic programming 
model to make the allocation. For some cases, it will be sufficient simply to allocate the unit to 


the depot which has the greatest probability of using it inthetime period until the next allocation. 
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In the initial provisioning of the depots, each depot will receive the safety stock S; computed 
using the rules derived earlier. If the amount which is normally on order is allocated to the 
depots, the simplest thing to do (which should also be quite satisfactory) is to allocate the units 
one at a time to the depot with the greatest probability of using it in the time T + 1/A. 










Redistribution 
It remains to consider the problem of redistribution. Since the on hand inventory levels 











































will generally be very low (0 or 1 usually), it will be desirable to consider whether or not there b 
should be a redistribution each time there is a demand in the system. The purpose of a re- b 
distribution is to protect against or reduce stockout costs in the time interval until the next 0 
allocation or in the time interval T + 1/A if nothing is on order. 0 
While it is permissible (provided that the redistribution time is substantially less than the Ww 
average time between allocations) in setting the stockage objectives at the depots and in making ni 
allocations to assume that a backorder will be satisfied by a redistribution and that the cost of S| 
a redistribution to j will be ° ’ these assumptions cannot generally be made when a redistri- @ 
bution is being considered. Recall that a redistribution is to be considered each time there is C 
a demand in the system. If the redistribution time is greater than the time until the next by 
allocation, no redistribution is considered. In such a case, if the demand caused a backorder, 
the next unit to arrive will automatically be allocated to the depot with the backorder, provided x 
that another backorder does not occur before the allocation time. di 
Consider now the case where a second backorder occurs in a time period Tu before a If 
unit is to be allocated. At the point in time when the second backorder occurs, the next unit to ar 
be allocated is tentatively allocated to one or the other of the backorders by minimizing stockout Ca 
and transportation costs under the assumption that a redistribution will be made to the other mi 
depot. After the allocation is tentatively made, a redistribution is considered to the remaining mi 
depot with a backorder. The time period to be used here is the time until the second unit or tin 
order is ready to be allocated or a time t + T + 1/A if only one unit is an order, where t < Tu of 
is the time until the unit on order is ready to be delivered. Of course, if a third backorder a 
occurs in the time period T)y before the arrivai of the first unit on order, the tentative alloca- 
tion is reexamined, and a redistribution is then considered to the depot not to receive the 
allocation. 
When a redistribution is considered, the above discussion has shown that if t = T + Tu 
is the time until the next allocation, then 7 >0. The next allocation here means the next one 
that is not tentatively already allocated; it may in certain cases be the time of arrival of the a 
second unit on order. The only depot considered as a possible receiver of a redistribution is, "In 
of course, the depot where the demand occurred. There are two cases to be examined. The be 
first is that where the demand caused a backorder at the depot j where it occurred. The “ 
other is that where it did not cause a backorder. 
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In this case the mode of transportation for redistribution can be selected at the outset. It is 
the mode which minimizes the sum of transportation and stockout costs i.e., minimizes 
A 


F =v+t uF yy It will almost always be true that the fast mode of transportation should be 


used in this an The depot i which is considering a redistribution to j must have at least 

one unit on hand; let L; be the number of units on hand at i. It will be assumed in making the 
computations that no units will later be redistributed to i unless there is a backorder Ifa 
backorder occurs at i before the next allocation, it will be assumed that a redistribution will 
be made by the mode of transportation decided upon above at a cost a if the redistribution 
occurs at a time at least T,, before the arrival of the next allocation. If the time of the back- 


M 
order is less than T,, before the arrival of the next allocation it will be assumed that the unit 


will be allocated to im the stockout costs will be at where t < Ty 38 the time until the 

next allocation. When more than a single backorder occurs at i before the next allocation it is 
supposed that all after the first will be met by redistribution regardless of when they occur, 
except in the case where more than L demands occur in the time period of length 7. In this 
case, the first demand, if any, occuring within a time t < Tu of the next allocation can be met 


by allocation rather than redistribution. 


The expected costs of stockout and redistribution for depots i, j will now be computed 
both for the case where i makes a redistribution to j and the case where i does not make a re- 
distribution to j. The expected costs at the other depots have no influence on the decision.® 
If i does not ship to j, the cost at j will be m (7+ Tw) since there is a backorder at j. Costs 
arising from additional backorders at j occuring later have no influence on the decision be- 
cause the redistribution does not modify then. At depot i stockout costs will be incurred if 
mM The cost of the first backorder (de- 
mand L. + 1) depends on whether it occurs at a time t > Tu before the next allocation or a 
time t < T,, before the next allocation. It will be a; if t > T,, and a,tift <T,,- The cost 
of any backorders after the first (demands L + 2, etc.) will be a except as noted previously, 


more than L units are demanded in the time 7+ T 


in the case where more than Li units are demanded in the time period of length rT. 


The probability that demand L + 1 occurs in the interval 0 to 7 is 


_— 


bin actuality, depots i, j cannot be completely isolated from the rest of the system. There 
are second order interactions, For example, if a backorder occurs at i ina time t < Thy 
before the arrival of the next allocation, then whether or not this backorder is met by alloca- 
tion or redistribution depends on whether or not some other depot has also had a backorder in 
the time period of length T), before the next allocation. These second order effects are 
ignored here since they have only a very slight influence on the expected costs. 









The most important case is that where the demand at depot j generates a backorder there. 
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where f  (u) is the Poisson probability that exactly u units are demanded at depot i in a time 
period of length 7. The probability that demand L, + 1 occurs between t and t+ dt, 
0 <t < Tu (t being time until the next allocation) is Ai (Lp dt. Thus the expected costs 
at depots i, j until the next allocation if no redistribution is made are’ 7 

x M 


bS f(u) +A, 7, \ t ft, (Ly) at 
0 


D = 1,(1+Ty4)+ % pe (o-L-1) +74, ; ' u=L,+1 
1 


u=L,+2 
i 


x 


T 
. M : 
(21) am}. fw) | t f¢(0) dt, 
u=Li+1 0 





it 


m L,+1 " i i 
t f, (L.)dt => f (u) - ik f_(u) 
A. 
still i? U=L+2 TM u=L,+2 rn 


oH 


x ioe 


(22) i) a ee ba f(u) |. 
i u=L,+1 M u=L+1 


In particular, 


T 
‘i. j ’ A.t 1 A.T 
i 2 a oe at 
(23) 2 t f, (O)dt s , te re "o E e i M (1 + AT )| js 


If a redistribution is made to j, the expected costs are 


Z a 4 @ Tu 
+ a » i >: i i 
u=L.+1 u=L, 0 
i 
i iol (Mg 
(24) +A TIS i) \ t f, (0) dt. 


0 


' The term = 
pa (u - L. - 1) fl (u) 
u= ;t2 i T+T yy 


is an approximation for the sum of three terms, the first giving the expected backorders from 
demands Lj + 2, ..., in the period of length T , the second giving the expected backorders 
from demands L; + 2, ..., in the period of length T,, given that demand Lj + 1 occurs in this 
period, and the third giving the expected number of M ckorders above one in the period of 
length T,y given that more than Lj occur in period of length T , given that one occurs in the 
period of length TM. This should be quite satisfactory, since the probability of more than 
two demands ina time T + Ty will be small. 















INVENTORY CONTROL MODEL FOR LOW DEMAND ITEMS 


Expected costs will be reduced by making a redistribution if 4D = D* - D <0. However, 


i) 


ne 
AD=2.-1 t+ ra (uy + 1 (0,)+ 4 fl (L,+1) - #£(L,+1) 
ee oe iu=L,+1 T+T yy i 4 i 5 T+T yy i “A i 


2) 
oo 


T. p : 
i i > , ; ' 
-——-i- f (u) - i i 
+ ws i* 2 ™+Ty u= ;+2 f(y) 7" iin (L,) ~ (4) 


Ty 


i i 
(25) +72, (1, ) { t f, (0) dt. 
0 
For each depot i with L >0, SD is computed. If every AD> 0, no redistribution should be 
made. When one or more4 D<0, a redistribution should be made, and the depot to make the 


redistribution is the one with the most negative 4D. 


The endeffect terms in Eq. (24) which arise because a backorder may be met by allocation 
rather than redistribution will yield expected costs significantly different from those obtained 
if it is assumed that all backorders are met by redistribution only when L, = 1 (the shipping 
depot has only a single unit) and the time at which the backorder occurs is close to the time 
Tu before the next allocation. These differences can, in general, be important however, since 


for the case under consideration it will usually be true that L =1. 


It remains to consider the case where the demand at depot j does not generate a back- 
order. Thus, after the demand, L, 20. This case will be studied in two parts: (a) that where 
L, = 0 after the demand and (b) that where L, > 0 after the demand. Part (a) will be examined 
first. Not only must it be decided whether or not to make a redistribution, but in addition the 
mode of transportation must be selected if a redistribution is to be made. By using the fast 
mode of transportation, the expected stockout costs in the time until the redistribution arrives 


can be reduced; however the fast mode is more expensive. 


The expected costs at depot j and the prospective shipping depot i will now be found 
for the case where i does not make a redistribution. As above, the time until the next allo- 
cation will be written t = 7 + Tw For depot j (which now has no units on hand), the stock- 
out costs incurred if there is more than a single demand in the time period t will not be 
considered since a redistribution could not influence them. It will be assumed that if there 
is a backorder at a time t before the next allocation Tu <t<?T+ Tw it will be met by re- 
distribution at a cost *, and if O0<t< Tw it will be met by allocation at a cost 7.t. For 
from depot i, it will be assumed that any backorder will be met by redistribution ata cost t The end 


effect where a backorder can be met at i by allocation will be ignored, since two demands must 
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occur to generate such a backorder, thus making its effect small. Therefore, the expected 


costs for i, j until the next allocation are: 






" . ta 
(26) D=a uL.) fo (ue Ccthiebin I t £,,,(0) dt. 
i u=L,+1 ( i T+T yy j u=1 T j j/0 T+t 






The expected costs if i does make a redistribution will be found next. For depot j, it 
follows that stockout costs are incurred only if there is a demand before the arrival of the 



















redistribution, since stockouts from more than one demand need not be included. Now, the 
end effect will be included for depot i rather than j since typically, i will have only a single 
unit on hand before it makes the redistribution. The expected costs at i, j if a redistribution 
is made are then 


eo a | Ty 
D* = T, Ze alot .. ty” + 7, iL, f(u) +A, 7, J, t f7,4(L,-1)at 
m Ts 
(27) +A, 2 Ao] rf t £,(0) dt + 7, 


where sf is the expected stockout cost at j before the arrival of the redistribution. 


yl 
i M he 1 “ATi 
7 = - = ——_— - 
" af 0 a, (Ti t) e dt ™ (Tr x)( e ) 
-A,T 
j M. 
(28) + ™ Tu e 


If AD = D* - D<0O, expected costs could be reduced by making the redistribution. The AD 
values are computed for each mode of transportation, and if all AD > 0 no redistribution 
should be made; if one or more AD <0, then a redistribution should be made by the depot and 
mode of transportation which minimize AD. 


Finally, part (b) will be studied. Inthiscase, thereare L > 1 units onhandat j after 
the demand. Only on very rare occasions woulda redistribution be made to j inthiscase. It 
would be done, only if j was one of the relatively high demand depots andthe time until the 
arrival of the next allocation was considerably greater than the average time between allocations. 
Here the end effect can be ignored and the cost of a backorder at j can be taken to be a. Now, 
i A 


a +) Ph j 
(29) AD=viT, f (u) - J wFL + fT 


i u=L,; T+ Tu 


(u). 
M 


If AD > Ofor all i # j, a redistribution should not be made. If AD <0 for one or more i, a 
redistribution should be made by the depot which minimizes AD. Here, it will almost always 


be true that if a redistribution should be made, the slow mode of redistribution should be used. 
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This completes the discussion of redistribution. A high-speed computer could go 
through this redistribution routine each time an item is demanded in a couple of seconds or 
less. However, the details could in many cases be carried out quite well by hand without 
going through all the details developed herein (just described), since often it will be fairly 
obvious which depot should redistribute to the needy depot, if no redistribution is made until 
there is a backorder. 


A Numerical Example - A system consisting of three depots will be used for the 





example. The following data are assumed: 


IC = $200 

m = $100,000 

uF = $2000 for j= 1, 2, 3. 

N =3 

A =3, Ay = 2, Ao = 0.6, Ag = 0.4 per year. 


Procurement Lead Time = 0.25 year 
Redistribution Lead Time = 0.02 year 
p =0.75, (round to 1). 


Under the assumption that system stockout costs are negligible, the "uncoupled" 
model can be used. Inequalities (Eq. (11)) with the Dj (u) of Eq. (4) are used to determine the 
safety stocks. The Dj (u) are tabulated in Table 1. 


























TABLE 1 
Depot 1 Depot 2 Depot 3 
p,(0) 0.6 0.833 0.883 
B() 0.24 0.1388 0.1039 
B,(2) 0.096 0.0231 0.0122 
p,(3) 0.00384 0.00386 0.00143 
B,(4) 0.00154 0.00064 0.00017 
Inequalities (Eq. (11)) result in optimal values’, Ss =e Sy =1; S3 =1. For p = 0.75 


(rounded to 1), Eq. (7) gives a system stockage objective, (on hand plus on order minus back- 
orders) of 5. Note that the uncoupling assumption was quite realistic since the probability of 


system backorders is only about 0.0001 per procurement lead time. 


7 ee : 

It might be noted that the same safety stocks would be obtained if pj (u) was taken to be Pois- 
son distributed with mean Aj/r, i.e., if the demand in the meantime between allocations was 
used. 
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To illustrate the nature of the dynamic programing allocation rule, it is assumed that 
there are two units on order, and that at the time the first unit is ready to be allocated (t, =0), 
the on hand levels are aT =1, bo = 2, §3 = 0. The second unit will be ready for allocation 
after 0.3 years elapse, (ty = 0.3). The lengths of the two periods under consideration are 
therefore %* 0.3, Ty = 0.28 (See Eq. (13).). 


At the start of the last period, the possible values of g= (é 1? 59> $9) are given in 
Table 2 along with the expected number of stockouts in the last period. The optimal allocation 
8 
of the unit allocated at to = 0.3, and the expected costs in the last period. 


TABLE 2 





1 0 





0 1 
































Optimal de- 
pot (for al- 
location at 
ty = 0.3) 1 3 1 1 1 2 1 1 1 1 2 1 3 2 1 


Expected 
Stockouts /0.025/0.025/0.145/0.25/0.035/0.155/0.155/0.305/0.41/0.14/0.26/ 0.26) 0 .035)0.14/0.13 

































































Expected 
Cost 50 50} 290) 500 70} 310} 310} 610) 820) 280) 520); 520) 70) 280) 260 
TABLE 3 
| 
(a) £& (0,2,0)}(1,1,0) (0,1,1) (,1,1)|(1,0,1) (1,1,0) | (0,2,1))(1,1,1) (1,2,0)|(1,2,1) | 





(b) Period 1 
Demand] (1,0,1)|(0,1,1) |(1,1,0) (2,0,0)|(0,2,0) |(0,0,2) | (1,0,0)|(0,1,0) (0,0,1)/(0,0,0) 


(c) Proba- Sir 
bility 0 .0294/ 0.00875) 0.0495) 0.0745) 0.00487 0.00277 0.246 |0.073 | 0.049 |0.41 


(d) No. first 
‘ all 

period 
Redis- rol 
tribu- 
tions 0 0 0 1 0 1 0 0 0 0 


(e) Cost of ~ 
first the 
period 
redis- 
tribu- is ¢ 
tions 0 0 0 2000 (0 2000 0 0 0 0 





















































To simplify computations some possibleé with very low probabilities have been neglected here. 
Specifically, cases with system demands in a single time period of more than two units are 
ignored. 
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TABLE 3 (Cont'd) 
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tions 


(f ) Expected 
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ty = 0 to depot 3 





0 .0294/0 00875) 0.0495) 0.0745 


2000 


1280 


minimum} 
| cost in 
period 2 
|(g) Costs x 
| Prob. 


67.0 


1 io «= 0S 2 
t ae pee 
| | | 
| 
| 
2000 |0 0  (/|0 4000 
t =- tT T = T —_———_+—— —— 
| | | | 
ee 
280 i {500 1500 | 280 
19.90 |13.85 |b7.25 2.43 [11.85 





1 


|0.00487| 0.00277 





+ 


| 
| 
| 





0.246 





10.073 |0.049 |0.41_ | 
ss } Recast Seve 














(0 0 1 0 

ae 

| | | 

0 0 2000 (0 

| i tn 
260 1140 |140 |50 

$$} —______} 

64.0 |10.20 |105.0 |20.50 











It remains to calculate the optimal allocation at t, = 0, 


given that the initial € is (1,2,0). 


the last row gives total expected cost for the two periods of $322.07. 


Since depot 2 is clearly overstocked relative to the other two depots, for simplicity only allo- 
cation to depots 1 and 3 is considered here. Table 3 gives (a) the values of £ at tg = 0.3 if 
allocation of the initial unit (t = 0) is to depot 3 when the period one demand is that given in 
The table also presents: (c) probability of occurrence; (d) the number of distributions 
in period 1; (e) the cost of period 1 redistributions; (f) the costs of period 2 for the specified 


$ vector; (g) the cost for periods 1 and 2 weighted by the appropriate probability. The sum of 


Table 4 provides the same data as Table 3 for the case where the unit allocated at t, =0 
is allocated to depot 1, i.e., the inventory after the first allocation is (2,2,0). The expected 


total cost is $351.98, about $29 more than that resulting from allocating the unit available at 


Hence, optimal policy is to allocate the first unit to depot 3. 


G. HADLEY AND T. M. WHITIN 
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ON MINIMAX SURVEILLANCE SCHEDULES*t 


C. Derman 


Columbia University 


INTRODUCTION 


Consider, as a unit, an operating system, a component of such a system, or a lot of 
material maintained in storage for use at some future unknown time. Suppose the unit under 
consideration deteriorates and that after a random period of time the unit is no longer capable 
of functioning properly. We shallrefer to this length of time as the life of the unit. In many 
cases there are estimable costs associated with the continuance in service (storage) of a unit 
which has failed. If the unit is a production system, the costs are associated with the amount 
of defective product produced; if it is material in storage (e.g., certain kinds of missile fuels), 
the costs are derived from considering the various implications of using, unknowingly, the 
unserviceable material; ... andthe like. Thus, in such cases, it is usual to inspect the 
unit from time to time in order to determine its condition. However, since the inspections 


are often costly, some consideration should be given to the scheduling of the inspections. 


We treat this surveillance-scheduling problem assuming that the unit remains in serv- 
ice (storage) until failure is detected by an inspection or it reaches some stipulated time T in 


service. We assume that a failure is detected by an inspection with probability p (p > 0). 


We do not assume any knowledge concerning the distribution function of the life of the 
unit. Our aim is to derive the schedule which minimizes the maximum possible expected 


cost associated with the inspections and failure of the unit. 


Of course, if the distributiou is known, then our minimax criterion is not meaningful. In 
fact, our particular formulation of the problem! was suggested by the one given by Barlow and 
Hunter [1] . In their formulation the distribution function F is assumed known and p= 1, 
and they derive the schedule which minimizes the expected cost associated with F. Roeloffs 
(in current work to appear in his doctoral dissertation) has obtained results relating to optimal 
schedules when partial knowledge concerning F is available. 


THE FORMAL PROBLEM 
Let Y, with unknown distribution function F, denote the life of the unit. By a surveil- 


lance schedule consisting of n observations we mean a set of n numbers x, (i = 1,005 BR) with 


‘Manuscript received November 18, 1960. 
‘Research sponsored by the Office of Naval Research, 


lour interest in this problem area originally grew out of conversations with Dorothy M. Gilford 
of ONR and representatives of the Bureau of Weapons. 
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o< xy a xy <"T; x denotes the time at which the ith inspection of the unit is made; T 


denotes the maximum length of time a unit is kept in service (storage). For a given 



























x= (x,, ees x) there is a function G, (Y) which is the cost associated with the unit having 


life Y. The expected cost for a given distribution function F is 


E,G 





G,(t)dF(t) ; 








Since F is unknown, we are concerned with 
M, = > E,G, 


where the supremum is taken over all distribution functions of nonnegative random variables. 


The problem is to find that x = x* such that 






w= min Hy = Hy, . 


THE MINIMAX SURVEILLANCE SCHEDULE 


We assume that the cost associated with the unit is 


cN + vdZ, 


where c > 0 is the cost of an inspection, N is the number of inspections on the unit, v > o 
denotes the cost per unit time attributable to the unit during the interval of time between its 
failure and its detection, and Z denotes the length of time between failure and detection. On 
noting that, for any given Y, the cost is dependent on the number of subsequent inspections 
needed to detect the failure and on taking the conditional expectation of the cost with respect to 


this aspect of the existing randomness, we have that G,(Y) is of the form, 


GAY) = gx, ee9X5 Y) 


= p(k+1l)c + pv(x_41-Y) 


n-l 
\ . 
+ ~ i (1-p)'*{ (ie tpe+vex,, -¥)} 
i=k+ 1 J 
+ (1-p)""*{ ne+(T-¥) | x. <¥ < x -k=0 nel 
i k ee) a at ? 
= ne + v(T-Y), x, <Y<T, 


(1) =. ney pee. & 
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where Xo = 0. 




















It is easily seen that for any distribution function F of a positive random variable, 


. n Xi 1 2 


2. | g(X,,--.+,X,,t)dF(t) + ne { dF(t) 
k=0 *k (sk mw!) T 


' 


E,G,(Y) 


IA 


n 
2, E(Xj soe Xy Xt 0) (F(x,, 1) - F(x) + ne(1-F(T)) 


1A 


max (g cracergh ka \ 
(80,x’ “is ) 


les. 
= Hy» 
where Bix = &(X;> rey X aX + 0), 2. ae 
and Xn+1 = T. 
Hence, our problem is to find that schedule x = x* such that 
— ye = min Ky = mijn max (fo, x? Selene En, x): 
ts 
On The result of this article is as follows: 
ons The minimax schedule x* is given by 
ct to a 
—. Cc ( a{ (n+1)p+ 2] _ 
(2) = IP lope i *iv\~ mer 7 Ot) , i=o, 1, »n, 


where n is the largest integer such that 


(3) ep-n” + cp(2-p)n + 2(c - pvT) < o> 


Proof: Let N be such that N(N-1) < 2T v/cp and n(n-1) > 2T v/cp, if n>N. Denote by 


Vn the class of schedules of n inspections, and let 


N 
-1, .. “Usci , 
i=O 


Define the schedule x= x(n) as that schedule xeC,, satisfying 


(4) Fox(n) ~ Six(n) = *** = Snx(n) = Mx(ny? 2 SN. 
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Let C* denote the class of schedules x(O), ..., x(N). 


(a) We show that if xeC. for n > N, then a, 2 Hy (N)’ Therefore if x* exists, it is in 


The proof will proceed as follows. 


Cy: (b) We show that x* exists and is a member of C*. (c) We show that Eq. (3) allows 


us to identify the member of C* which is x*. 


Proof of (a). In order to exhibit x(n) explicitly, note that for any schedule x «€ Cy on rear- 


ranging terms, 


n-1 - 
(5) Sx = ¢ ‘2 cektcn)Keatt-n? 
i=k 
if 
“ss >? aan, k=0, ...,n-1, 
4 
nx = cn+v4,,- 

















where 4; = x1,4-X,, i=0,1,...,n. From (5) we obtain, 


(6) (1-p)g., + vO)_1 + ckp = Bc-1) x’ | oo) ae 


On solving (4) in (6), we get x;(i=1,...,n) given by (2), and 


vT + > {(n+1)p+2 } 


Px(n) = — 





a ‘fo | ae. B 
np + 1 


Because of the constraints 4,70 (i=0,1,...,n), x(n) is not defined if n>N; for, if so, 4,40: 


Now, let x be any schedule in Cc. where n>N. Trivially, My >cn. Hence, 


7 vT + x {(N+1)p+2} 
Mx P(N) 2 2 ee ee 
Np + 1 
vT + {(N+1)p+ 2} 
> ¢e(N+1) - = 





Np + l 


P N(N+1) - vT+c 
Np * 3 
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cp 2Tv 





— -w+ec 
> 2 & 
Np + 1 
= 
Np+1 
20. 
Proof of (b). Fixn<N. If n=o, X(o) is the only schedule. For n> 0, suppose x is a schedule 
of n inspections such that max ax 5. conele Enk ) < Mein)’ Let 


n 
5. : 2 Ekx * 
k=o 


Our supposition implies that Ss: < S(n) . From the last equation of (5) and those of (6), on 


adding, we have for all x eC 


nc 
p S. + (1-p)g5x = vT + “9 (p(n+ 1) + 2). 


However, Since the right side of the above equation is constant, S; < S(n) implies g.5 7 Box(n)? 
a contradiction. Therefore x(n) is a minimax schedule in C.: This, together with (a) and 


the finiteness of C* proves (b). 


Proof of (c). Notice that 





(7) ; _ cp’n? + (2cp-cp’)n + 2Ac-pvT) <x 
x(n) ~ Yx(n-1) ,n<N. 
| awe 2(np+1) ( (n-1)p+ 1) 
The sign of Eq. (7) is the sign of its numerator; therefore Hy (,) 18 a nonincreasiny function 


of n as long as the numerator of Eq. (7) is nonpositive. The numerator, however, is an 
increasing function of n and nonnegative for n=N. Therefore, the n such that 
x* = x(n) is given by Eq. (3) Q.E.D. 
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PROBLEMS 


It has been suggested that this journal might serve as a medium for publishing problems 
in the area of Logistics -- with the idea that other persons might be interested in those prob- 
lems and might submit comments. Readers are invited to submit brief statements on applied 
and theoretical problems in Logistics. Address letters to Managing Editor, Naval Research 


Logistics Quarterly, Office of Naval Research, Washington 25, D. C. 


ON PROFESSOR SCHELLING'S STRATEGY OF CONFLICT 





John G. Cross 


Princeton University 


Professor Schelling's Strategy of Conflict[ 1] can be divided conceptually into three 





parts: an extensive verbal analysis of the nature of Threats, Commitments, and Promises 

and their uses in bargaining; a theory which may help one to understand the forces underlying 
the outcomes of bargaining situations; and a series of criticisms of Game Theory regarding 

its analysis of human behavior in conflict situations. In general, the terminology and utility 
matrices common in Game Theory are used, and in fact one may consider the Strategy of Con- 
flict [1] to be an effort to bridge the gap between Game Theory analysis and actual develop- 
ments under conflict. 

The most fruitful part of the book is the excellent description of the use of threats, 
commitments, and promises, and the requirements which must be met before they can be put 
into operation successfully. Many provocative and entertaining examples are given, some of 
which shed considerable light on very ordinary experiences. The importance of commitments 
is emphasized, both as tactics in themselves (which can force opponents into taking the desired 
course of action), and as an important means of convincing an opponent that.one intends to 
carry out a (perhaps mutually painful) threat if necessary. 

Unfortunately, Schelling | 1] fails even to consider an important objection to the use of 


” , ] , . 
such ''strategies"’ as threats and commitments: that of unfavorable psychological reactions to 





h hesitate to use the word ''strategy"' for fear of causing confusion: Schelling's definition (Y of 
a strategy as something concerned with the ‘exploitation of potential force ( 1 » Pp. > 1 0) 

1s cuite different from the game theoretical definition of a player's preplay plan which speci- 
fies what choices he will make in every possible situation which may arise during the course 
ofa particular game (see [2] p. 79.) 
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them. The use of a threat, for instance, may very well antagonize an opponent (i.e., modify 
his utility values associated with the possible outcomes) so much that the threat does not 
work — to the disadvantage of both players, since the threat may have to be carried out. It 
would seem that such effects would have to be taken into account in any complete analysis of 
bargaining. 

Schelling turns next to a discussion of some factors which can influence the specific 
outcomes of bargaining, particularly in cases of "'tacit'’ bargaining in which the two players 
have no opportunity to communicate. He postulates that under such circumstances, the out- 
come is likely to be at some point of prominence (a "focal" point) of a mathematical or physical 
nature (such as a 50-50 split of a dollar, or a military division of land at some natural bound- 
ary). This analysis is primarily on an intuitive basis, in spite of an "informal" series of ex- 
periments which Schelling reports in order to substantiate his position. His analysis of why 
bargaining tends to settle at these focal points is rather sketchy — his conclusion that focal 
points provide obvious points at which a player can make a strong stand is accompanied by the 
remark that: ''This proposition may seem intuitively plausible; it does to the author, and in 
any event, some kind of explanation is needed for the tendency to settle at focal points."* 

In order to discuss Professor Schelling's criticisms of Game Theory, it is first neces- 
sary to distinguish between Game Theory proper and some of the more recent extensions to 
Bargaining Theory. The former, as originally presented by von Neumann and Morgenstern 
is strictly an analysis of the range of possible outcomes of conflict situations, based on a few 
specific axioms. As it turns out, these axioms are usually not sufficiently restrictive to yield 
one unambiguous outcome for thegame. For example, in the two-person, nonzero sum game, 
Game Theory concludes that "each player gets individually at least that amount which he can 
secure for himself, while the two get together precisely the maximum amount which they can 
secure together. we Thus, although we do know that the players will act in combination to ex- 
tract the greatest possible gain (from ''Nature''), we cannot say how this gain is divided between 
the two. The major contribution of Game Theory in this case is an elaborate framework into 
which new axioms can be introduced easily, to restrict further the number of outcomes. Bar- 
gaining theories consist of just such axioms, an example being Schelling's own assertion that 
focal points or prominences which are likely to have some subjective appeal to the players will 
affect, or even determine the outcome of the game. Unfortunately, most of these additional 
axioms involve much stronger assumptions about the behavior of the players, and therefore 
they are much more difficult to justify than those of Game Theory. Insofar as Schelling's 
criticism on page 119 [1| that Game Theory has been "pitched at too abstract a level of analy- 
sis" (referring to its inability to point to specific outcomes), we may reply that Game Theory 
never intended to give such results, its purpose being only to restrict the range of possible 


solutions as far as possible, using only its relatively weak assumptions, and leaving further 


ts ; a2 4 
restrictions to psychological analyses of bargaining. 


*See [1], ae ie 
>See [2], p. 555. 
“See [2] » p. 616. 
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A further issue stems from Schelling's assertions about Game Theory's "preoccupation 
with the zero sum game" and his many subsequent comparisons between nonzero sum games, 
bargaining, and the Game Theory analysis of two-person, zero sum games, particularly 
regarding the importance of communication between the two players: "There is probably no 
contrast more striking, in the comparison of the mixed motive [Schelling's term for a nonzero 
sum game] and the pure conflict (zero sum) game, then the significance of having one's own 
strategy found out and appreciated by the opponent. Hardly anything captures the spirit of the 
zero Sum game quite so much as the importance of "not being found out" and of employing a 
mode of decision that is proof against deductive anticipation by the other player."' 5 

It is difficult to understand how a theory which emphasizes the existence of incomplete 
and complete information, the "bluff," © and the coalitions between players which must occur 
before the game can be ''won,"' can be accused of basing its analysis on secrecy, or how a 
theory originally designed to treat economic problems "in which the sum of all payments, the 
total social product will in general not be zero and not even constant"? can be described as 
"preoccupied with the zero sum game." It is true that Game Theory bases its theoretical 
development on the two-person, zero sum game, and upon the mixed strategy, but not in the 
restricted sense in which Schelling has presented it. The "persons" in the Game Theory 
sense are merely groups of players who have combined, or who are classified together for the 
purposes of analysis (e.g., into the "winning" and the "losing" coalitions). These groups do 
use minimax strategies against one another when playing an n-person zero sum game. 
Furthermore, Game Theory is enabled to make use of its zero sum game analysis through the 
creation of a "fictitious player'' who takes a purely passive role in the game, paying out what- 
ever is demanded by the rest of the players according to the rules of the game — again only a 
mathematical device designed to simplify the analysis. In this manner, the two-person, non- 
zero sum game is treated as a three-person game in which the third player (''Nature") is 
exploited by the other two players in combination. The specific outcome, of course, depends 


tween 


ie (among other things) upon the rules of the game, communication, and bargaining between the 


7 two players. At no point in this analysis, however, do we have the two "real" players using 


- minimax strategies against one another. Professor Schelling has quite inappropriately used 


will the two-person, zero sum analysis as an example of the Game Theory approach to any two 
individuals in any game regardless of what other (perhaps "fictitious") players are also in it. 
Schelling goes on to point out that the threats, promises, and commitments which he 


has demonstrated to be so important in bargaining are not explicitly considered in Game 


Mee [1], p. 160. At this point Schelling quotes the Theory (op. cit., p. 147) in a footnote, 
quoting a passage emphasizing the avoidance of being found out. This remark, however, 
tefers only to certain zero sum games, and not to two players who have some interest in 
Joint operations, 


Schelling's concept [1] of a "bluff, '' a misrepresentation or distortion of one's own value sys- 
tem (see [1] pp. 23, 33, 36), differs somewhat from the Game Theory definition of an impedi- 
ment to the spread of information (concerning the state of affairs in a game in general or 
simply the strength of one's position) to an opponent. (See [2] » pp. 54, 188.) 


‘See [2], p. 46. 





424 NOTES 


Theory. This is because Game Theory generally treats the normalized form of the game, into 
which, as Schelling himself demonstrates on pages 150-156, these concepts can all be incorpo- 
rated. Schelling complains that these tactics cannot be studied in this form, which is quite 
true, but Game Theory was not designed to study such things in themselves, but only their 
effects upon the possible outcomes of the game. Game Theory begins with the assumption 

that all existing alternatives including threats, promises, and the like, are already known, 
and are included in the normalized game. 

Finally, Schelling criticizes Game Theory on the grounds that the normalization and 
randomization processes obliterate the "prominences" and "focal points" upon which his bar- 
gaining analysis depends.° For instance, he presents the matrix: 

Player B 
II Ill 


0,0 | 0,0 


ii | 0.0 | 10,10 | 0,0 
iii | 0,0 9,9 


— + 


L_0,0 


Player A 


iv 0,0 | 0,0 


i 


and states that if player A and player B are kept incommunicado, they will tend to choose iii 


III on the grounds that this payoff acts, through its un'queness, as a focal point and thus pro- 
vides a "clue to concert choices.'"' Presumably, many other such "clues" which are inherent 
in the structure of the game have been lost through normalization. In the first place, it is 
important to note that the figures given in the Game Matrices are utilities and not physical 
units. ? Therefore, in the many cases in which physical objects serve as the "prominences,"’ 
these utilities cannot be used as focal points, and thus the prominences are kept outside the 
matrix. There is nothing to prevent us, however, from considering the fact that the strate- 
gies (I, II, Ill, IV, i, ii, iii, iv) possess prominences external to the matrix. In fact, 
Schelling's theory of focal points in which each player has some reason (which may or may not 
be purely intuitive) to assume that the other will choose a particular strategy, can be inter- 
preted as the assignment of a probability distribution to each player's set of strategy choices by 
the other, who then acts accordingly. Hence the problem is reduced to one of decision making 
under risk — for each player. For instance, referring to the matrix above, player A may 
assume (subjectively) that player B will choose strategies I, II, III, IV in the ratio 1:1:2:1. 


The expected value of each of his strategies will be in the ratio: 


10 for strategy i 
10 for strategy ii 
2x9 = 18 for strategy iii 
10 for strategy iv 
Hence he will chose iii. Whether player B assumes the same, or goes further and assumes 


that A makes the above calculation will make no difference — the outcome will still be iii III. 


9v, 156. 


es in many cases assume money to be a direct measure of utility. (See | 2] pe 79.) 
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This interpretation also demonstrates how a focal point may not always be the solution. 


Consider the matrix of payoffs: 


Player B 

Il IV 
0,0 0,0 
3,3 0,0 
0,0 , 0,0 
“| 0,0 3,3 








Player A 





























If A again assumes that B will choose his strategies in the ratio 1:1:2:1, the expected values 
of his strategies will be in the ratio of: 

3 for strategy i 

3 for strategy ii 

2 for strategy iii 

3 for strategy iv 
and A will not take iii, preferring a purely random strategy among i, ii, iv. All this simply 
serves to demonstrate how it is trivially simple to consider focal points within the context of 
Game Theory. 

Schelling's assertion that Game Theory neglects the essential labels or prominences 
which make the above procedures possible evidently arises from his conclusion that Game and 
Bargaining Theories assume complete symmetry. It is true that certain forms of mathemati- 
cal symmetry are assumed, but certainly such symmetry as would affect Schelling's results 
(i.e., that both players have exactly the same role in the game) is treated as a special case. 
Clearly, if we did assume that both players have identical strategic opportunities, or that 
both players have identical personalities and physical characteristics, many of Schelling's 
"labels"’ would be lost. It should be pointed out here, though, that these labels must be con- 
sidered to be a part of the normalized game structure in the form of additional strategies, 
and that therefore the idealized "player'' A or B does not correspond precisely to the player 
George Washington Roosevelt (who counts among his ancestors Capt. Miles Standish, Benjamin 
Franklin, and Henry VIII of England), who is running for State Senator against Joseph Glotz 
from East Howesville. Thus, only an assumption of identical strategic opportunities for the 
players A and B would preclude Schelling's analysis, and it is certain that in general, Game 


Theory makes no such assumption 
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NEWS AND MEMORANDA 


Readers are invited to submit to the Managing Fditor items of 


general interest in the field of logistics 


SYMPOSIUM ON OPTICAL CHARACTER RECOGNITION 


A symposium on Optical Character Recognition will be held on 15, 16, and 17 January 
1962. It will be jointly sponsored by the Information Systems Branch of the Office of Naval 
Research and the Research Information Center of the National Bureau of Standards. The Sym- 
posium will be held in the Department of the Interior Auditorium on C Street, between 18th and 
19th Street, N.W., Washington, D. C. 

The objective of the Symposium will be to bring together equipment specialists, scientists, 
and potential users who are interested in exploiting automatic character recognition techniques 
for a variety of applications. Major purposes are to explore, for mutual benefit, the present 
state of the art, operational requirements and the implications of current research. It is 
planned that the meeting will be composed of several sessions to cover, respectively, presen- 
tations on available equipment, directions in current research, representative requirements 
of potential users, and a forecast of future progress. The program for these sessions will 
consist of invited papers from individuals and organizations active in research and development 
efforts in the field, and from potential users. Where feasible, demonstrations of optical 
character recognition systems in operation in the Washington area will be arranged. 


Attendance is open to all interested technical and management personnel. Further in- 


formation anda preliminary Symposium program, whenavailable, may be obtained by coatacting: 


Miss Josephine Leno 
Code 430A 

Office of Naval Research 
Washington 25, D. C. 
(Phone OXford 6-6213) 


DONALD K. POLLOCK, ONR 
BERNARD RADACK, BUSANDA 
MARY ELIZABETH STEVENS, NBS 


Symposium Co-Chairmen 





NEVS AND MEMORANDA 


SYMPOSIUM ON REDUNDANCY TECHNIQUES FOR COMPUTING SYSTEMS 


A symposium on Redundancy Techniques for Computing Systems will be held on 6 and 7 
February 1962; sponsored by the Information Systems Branch, Office of Naval Research. This 
Symposium will be held in the Department of the Interior Auditorium, C Street between 18th 
and 19th Street, N. W., Washington, D. C. 


The objective of the Symposium is to focus attention toward new ideas, research and 


developments which may lead to the introduction of redundancy techniques into forthcoming 


computing systems. It is apparent that many information processing systems of the future 
will be of such size that conventional methods of fabrication and of maintenance will be com- 
pletely infeasible. Computers for space applications have reliability requirements which aré 
presently not achievable without use of redundancy. Thus, some form of logical redundancy 
must be introduced if practical systems are to be operated. Accordingly, it is planned to 
present in a single symposium a collection of related papers concerning suggested new tech- 
niques for the use of redundancy in computing systems. This should permit relative evalua- 
ticn of such techniques with respect to their probable future utilization in various aspects of 
the computing field. The program will consist of papers invited from many of the organizations 
engaged in appropriate research and development activities. 

Attendance is open to all interested technical personnel. Further information and a 


preliminary Symposium program, when available, may be obtained by contacting: 


Miss Josephine Leno 
Code 430A 
Office of Naval Research 
Washington 25, D. C. 
(Phone OXford 6-6213) 
RICHARD H. WILCOX 


Symposium Co-Chairman 
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MODERN COMPUTING METHODS, Second Edition (completely revised) Philosophical 
Library, NEW YORK, pp. vi + 170, $6. 

This is the second edition of the book initially written by L. Fox, E. T. Goodwin, J.G. 
L. Michel, F.W.J. Olver and J.H. Wilkinson. To write this edition, C.W. Clenshaw, D.W. 
Martin, and G. F. Miller have collaborated with the 2nd, 4th and 5th auihors of the first 
edition. 

The book contains 15 chapters together with a comprehensive bibliography. It is well 
written and contains illustrative examples. It is excellent reading for mathematicians and 
scientists who come in contact with computing. 

The first five chapters discuss problems associated with different, computationally 
useful methods of solution of linear equations and latent root problems of matrices. Elimina- 
tion methods e.g., Gaussian elimination and variants of it are discussed along with iterative 
processes, Jacobi's, Givens’ and Householder's methods, for finding latent roots, with an 
emphasis on computation. Error analyses (e.g., for the Gaussian elimination method and for 
triangular decomposition) are given. In chapter six, methods of finding zeros of polynomials 
are discussed (e.g., Newton's rule, Bairstow's process and other iterative methods). The 
remaining chapters are on finite differences, Chebyshev series ordinary and partial differen- 
tial equations. use of recurrence relations for evaluating limits, the evaluation of integrals and 
the tabulation of mathematical functions. 

The book treats the topics covered in an interesting fashion with greatest concern for 
computational aspects but often developing the required theory for background material of 
which of particular interest are the presentations on matrices and on hyperbolic partial differ- 
ential equations. 

The subject of constrained maxima and minima is much in demand and occupies the 
attention of many people. It is not discussed in the book nor does the book contain general 
philosophical background which points out the nature of classes of general problems and their 
associated computational subtleties. The readers of the Quarterly will find it a welcome addi- 
tion to a subject much in demand, particularly for its broad coverage of topics and treatment of 
computational problems in an appealing way. 


(Reviewed by Thomas L. Saaty) 


ELEMENTS OF MATHEMATICAL STATISTICS, by Howard W. Alexander. John Wiley 
& Sons, Inc., New York, 1961. xi + 367 pp. 


In spite of the publication of many books on statistics, both students and teachers of 


Statistics are having difficulty in finding a textbook which is appropriate to the mathematical 
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backgrounds of students. In this sense we appreciate the author's intention to ''make mathe- 
matical statistics really accessible to the college senior or junior whose background includes 
a typical calculus course." 

In this book, Professor Alexander takes up adequately almost all statistical methods— 
but he does not refer to the concept of decision theory—without showing any of the bias which is 
apparent in several books, according to the authors' predilections. I think this impartiality, 
is shown in the following chapters in Alexander's book, is an important qualification for a 
good textbook. 

Chapter 1 Discrete Distributions 
Continuous Distributions 


Sampling Distributions I 


Sampling Distributions II 


Regression and Correlation 


2 
3 
4 Statistical Inference 
5 
6 
7 


Analysis of Variance 
Appendix Matrices 
His introduction of the concept of probability is especially excellent. In teaching sta- 
tistics, we have difficulty at times in proceeding from the discrete distribution to the continuous 
distribution. The author of this book gives a very natural and understandable explanation of this 
point. 
Generally speaking, I think the author has succeeded in realizing his purpose. None- 


theless, I have some questions as to his presentation from the pedagogical point of view: 


1. When we consider the importance of the concept of the expected value in many 
applications of statistical methods, it seems to me necessary to give an empirical explana- 
tion of this concept, in relation to the concept of probability as it is given in this book. We 
find many people who are applying this concept superficially to situations where repetition 
of the trial is not possible. 

2. In testing hypotheses, the author does not distinguish clearly between simple hypoth- 
eses and composite hypotheses. If we were to follow the traditional approach to statistics—as 
the author has, as shown in his neglect of the explanation of the statistical decision function— 
it seems necessary to make clear the difference between these two kinds of hypotheses. 

3. Iam afraid that Alexander's explanation of the difference between the concepts of 
the confidence coefficient and probability will confuse students. 

4. It seems necessary to give an empirical explanation of the meaning of the unbiased- 
ness in the point estimation. This point is related to my question 1. 

5. When we consider the powerful applicability of the least squares method, it seems 
unreasonable to omit the proof of the unbiasedness of the least squares estimates without 
assuming normality, when the proof of the distributions of these estimates under the assumption 
of normality is given. 

Notwithstanding these questions, I think this is a good textbook on statistics. 


(Reviewed by Koichi Miyasawa) 
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